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Abstract 


An  experimental  research  program  based  on  laboratory  test  studies  and 
scaled  slope  model  tests  was  conducted  with  specimens  of  Monterey  No.  0 
sand.  The  principal  objective  of  the  research  was  to  study  the  effects  of 
interactive  coupling  during  combined  compression  (normal)  and  shear  loading 
on  the  response  of  sands  to  dynamic  loading.  The  research  program  included 
the  following  experimental  studies: 

1.  Resonamt  Column  tests  on  cylindrically  shaped  specimens  using 
longitudinal  and  torsional  excitation,  both  separately  and  in 
combination . 

2.  Thin-walled  Hollow  Cylinder  tests  using  longitudinal  and  tor¬ 
sional  cyclic  excitation,  both  separately  and  in  combination. 

3.  Large  scale  shaking  table  tests  on  slope  models,  using  horizon¬ 
tal  and  vertical  cyclic  excitation,  both  separately  and  in  com¬ 
bination. 

During  the  course  of  this  research  significant  effects  were  observed 
during  combined  compression  and  shear  cyclic  loading.  The  primary  ob¬ 
served  effect  of  combined  loading  was  the  more  rapid  degradation  of 
modulus  with  strain  than  would  otherwise  occur. 

Two  methods  were  developed  and  presented  for  calculating  the  degra¬ 
dation  of  compression  modulus  with  strain  under  combined  dynamic  loading 
conditions.  The  first  of  these,  called  the  Strain  Ratio  Method,  requires 
the  computation  of  either  the  instantaneous  or  an  overall  average  ratio 
of  shear  strain  amplitude  to  compression  (normal)  strain  amplitude.  The 
amount  of  additional  degradation  in  compression  modulus  due  to  interactive 


coupling  may  then  be  determined  by  reference  to  a  set  of  typical  strain 
ratio  curves  which  are  presented  in  Figure  7-7  on  page  138. 

The  second  method  developed  is  called  the  Octahedral  Shearing  Strain 
Method.  This  method  requires  the  computation  of  either  the  instantaneous 
or  an  overall  average  value  of  the  octahedral  shearing  strain  amplitude. 

This  may  be  calculated  from  the  strain  tensor  using  either  Equation  2-42 
on  page  18,  or  Equation  2-44  on  page  19.  Once  this  value  is  determined, 
the  total  degradation  in  compression  modulus  due  to  strain,  including  in¬ 
teractive  coupling  effects,  may  be  determined  by  reference  to  a  set  of 
typical  octahedral  shearing  strain  curves  such  as  those  presented  in 
Figure  7-6  on  page  136.  Alternatively,  a  specific  set  of  octahedral 
shearing  strain  curves  may  be  developed  for  any  sand  by  converting  the 
strain  amplitude  from  conventional  degradation  curves  to  octahedral 
shearing  strain  by  use  of  Equations  2-42  or  2-44  as  outlined  in  the  text. 

These  two  methods  provide  a  reasonable  estimate  of  the  effects  of 
interactive  coupling  on  the  degradation  of  modulus  with  strain.  Both 
methods  lose  accuracy  as  the  straining  progresses  from  elastic  to  plastic, 
and  as  specimens  approach  failure.  Nevertheless,  the  use  of  either  method 
under  conditions  of  combined  shear  and  compression  loading  represents  a 
significant  improvement  over  the  practice  of  neglecting  interactive  effects. 

A  series  of  large  scale  shaking  table  tests  were  conducted  upon  slope 
models  in  an  effort  to  determine  when  yielding  in  the  slopes  began.  The 
slope  specimens  were  subjected  to  horizontal,  and  combined  horizontal  and 
vertical  cyclic  excitation.  The  slopes  were  thoroughly  instrumented  to 


An  analysis  method  was  developed  and  presented  for  predicting  yield 
accelerations  in  granular  slopes  under  combined  vertical  and  horizontal 
loading  conditions.  When  the  shaking  table  test  results  were  compared 
with  the  predicted  values  of  yield  acceleration,  it  was  concluded  that  the 
predictions  were  at  least  as  accurate  as  the  empirical  measurements  could 
be  made.  These  results  were  consistent  with  the  conclusion  that  the  ef¬ 
fects  of  combined  vertical  and  horizontal  accelerations  on  yielding  of 
granular  slopes  may  be  evaluated  using  simple  superposition. 
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Introduction 


The  objective  of  this  research  was  to  study  the  interaction  effects 
of  combined  compression  and  shear  loading  on  the  response  of  sands  to 
dynamic  loading.  This  research  included  the  following  experimental  studies 

1.  Resonant  Column  tests  on  cylindrically  shaped  specimens  of  sand 
using  longitudinal  and  torsional  excitation,  separately  and  in 
combination  (Chapter  5) . 

2.  Thin-walled  Hollow  Cylinder  tests  on  sand  specimens  using  longi¬ 
tudinal  and  torsional  cyclic  excitation,  both  separately  and  in 
combination  (Chapter  6) . 

3.  Large  scale  shaking  table  tests  on  slope  models  constructed  with 
sand,  using  horizontal  and  vertical  cyclic  excitation,  both 
separately  and  in  combination  (Chapter  8) . 

A  review  of  available  theories  for  analyzing  soil  behavior  under 
dynamic  loading  conditions  was  made,  and  the  strain  tensor  (Chapter  2) 
and  state  of  stresses  (Chapter  3)  were  developed  for  the  specimens  under 
test  in  the  first  two  of  these  studies.  Formulae  for  calculating  moduli 
and  damping  factors  from  the  results  of  those  studies  were  developed  and 
presented  in  Chapter  4.  The  results  of  these  first  two  studies  were 
further  analyzed  and  combined  in  Chapter  7. 

The  details  of  the  large  scale  shaking  table  testing  program,  in¬ 
cluding  the  development  of  an  analysis  technique  and  evaluation  of  the 
test  results,  were  presented  in  Chapter  8. 

The  mechanical  and  testing  details  of  the  various  equipment  used  in 
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these  experimental  studies  are  presented  in  Appendix  A.  Derivations 
(Appendix  C) ,  computer  programs  (Appendix  D) ,  and  example  test  results 
(Appendix  B)  are  also  included. 

These  experimental  studies  were  performed  on  specimens  constructed 
of  Monterey  No.  0  sand,  a  uniformly  graded,  fine  grained  quartz  sand  pro¬ 
cessed  from  beach  sand.  A  gradation  analysis  of  the  sand  used  in  these 
studies  is  shown  in  Figure  1-1.  The  maximum  and  minimum  densities  of  this 
sand  were  found  to  be  as  follows: 


v  = 

d  max 

1.707  gm/cc. 

(1.1) 

and 

V  = 

d  min 

1.425  gm/cc. 

(1.2) 

Experimental  studies  were  performed  at  a  variety  of  densities  and  con¬ 
fining  pressures  so  that  the  effects  of  these  variables  could  be  evaluated 
as  well. 
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Chapter  2 
The  Strain  Tensor 

Introduction 

In  both  the  triaxial  resonant-column  and  the  thin-walled  hollow 
cylinder  testing  series  forces  and  displacements  are  measured  in  the 
vertical  and  torsional  directions.  However,  because  of  the  geometry  and 
complex  loading  of  these  specimens,  the  principal  stresses  and  strains 
are  not  directly  measurable.  It  is  the  purpose  of  this  chapter  to  de¬ 
scribe  the  strain  tensor  for  these  two  testing  series  so  that  the  com¬ 
plete  state  of  strains  may  be  determined  from  those  strains  which  are 
directly  measurable. 

During  the  two  testing  series,  the  maximum  strains  measured  were  in 
the  order  of  10  !%,  so  that  a  linear  distribution  of  strains  may  be  as¬ 
sumed.  An  elastic  analysis  of  the  stress-strain  response  of  these  soil 
specimens  may  be  performed  with  the  use  of  Hooke's  Law. 

Hooke ' s  Law 

In  its  simplest  form,  Hooke's  Law  for  a  homogeneous,  isotropic  body 
may  be  expressed  as  follows : 
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The  use  of  Hooke's  Law  for  modeling  the  stress-strain  behavior  of  a 
material  implies  several  assumptions.  The  material  is  assumed  to  behave 
linearly  and  elastically  during  loading,  so  that  the  strains  considered 
must  be  small  enough  so  that  changes  in  shape  and  size  are  negligible. 
Hooke's  Law  further  implies  that  superposition  of  loading  effects  is 
valid,  and  that  the  principal  stresses  and  strains  are  in  the  same 
directions . 

Hooke's  Law  in  Plane  Strain 

For  the  special  case  of  plane  strain  in  the  xy  plane,  Hooke's  Law 
reduces  to  the  following  expressions. 


Y 


yz 


0z  can  be  expressed  as: 


a 

z 


(2.2) 


(2.3) 


And: 


u  •  (a  +  a  ) 
x  y 
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|a-u2) 

'  !<U+U2) 


'  |(u+u2) 

|aV) 


2 ( 1+U) 
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Resonant-Column  Testing 

The  strain  tensor  for  the  triaxial  resonant-column  testing  series  may 
be  expressed  in  "cylindrical"  coordinates  in  the  following  form: 


£r  £r9  £rz' 


£  £8r  £8  £9z 


e  „  e  o  £ 

ZR  ZD  Z 


where  the  R,  9,  and  z  directions  are  as  shown  in  Figure  2-1 (c)  and  £0z  * 
Yq2»  £Rg  *  YRg.  ezR  “  YzR-  To  maintain  equilibrium,  we  can  see  from 


Figure  2-1  (b)  that  | €rQ |  =  |e0R|, 


i£zr'  '  and  IS91 


|£0z|;  and 


so  [e]  must  be  a  symmetrical  matrix  which  is  uniquely  defined  by  six 
strains . 

Strain  Equations 

If  the  displacements  u,  v,  and  w  are  defined  as  shown  in  Figure 
2-1 (c),  the  strains  can  be  defined  as  follows: 


u  dv 
r  +  d8 


(2.8) 


(b) 


FIGURE  2-1  FREE  80DY  DIAGRAM  AND  STATE  OF  STRAINS 
FOR  RESONANT  COLUMN  SPECIMENS 
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R0 

1.1  du 
2  r  d0 

dr  r 

6r 

z9 

=  £.  = 

+ 

2  'dz 

r  d0 

9z 

... 

h—  + 

2  az 

dw 

dr' 

Rz 

zR 

(2.9) 

(2.10) 

(2.11) 


For  this  testing  series,  the  following  boundary  conditions  exist: 

1.  w  =  a  •  z,  a  =  €  (2.12) 

z 

2.  v  *  r  •  z  •  A8q  (2.13) 

Equation  2.12,  states  that  the  vertical  strain  is  known.  This  is  true 
because  it  is  directly  calculated  in  all  of  the  tests.  The  assumption 
expressed  in  Equation  2.13  is  that  the  torsional  displacement  v  varies 
linearly  in  both  the  r  and  z  directions.  This  condition  is  illustrated 
in  Figure  2-1 (a). 


•  Strain  Tensor 


Applying  boundary  conditions  1  and  2  (Equations  2.12  and  2.13)  to 
the  six  strains  which  define  the  strain  tensor  (Equations  2.6  through 
2.11)  yields  the  following: 


£ 


(2.14) 


A  detailed  derivation  of  the  above  expression  appears  in  Appendix 

C-l. 


A 
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Hollow  Cylinder  Testing 


In  the  case  of  the  thin-walled  hollow  cylinder  testing  series  the 
strain  tensor  may  be  expressed  in  the  same  form  as  Equation  2.5,  but  the 
strains  are  now  as  shown  in  Figure  2-2. 


Strain  Equations 


As  with  the  resonant-column  testing  series,  [el  is  a  symmetric 
matrix  about  the  normal  strain  diagonal,  and  is  uniquely  defined  by  six 
strains.  From  Figures  2-2  and  2-3  those  strains  can  be  defined  as  follows 

dw 


dz 


(2.8) 


fl  ~  RAVG/r>  /  R2NEW  \ 

'6  _  L(1  '  WV  *  1  R20LD  ’  4 


(2.15) 


=  -e. 


R  <r<R 
1  AVG-  -  2 


“l^AVG 


.  _  R2NEW  ~  RAVG 

R  -R20LD  "  RAVG 


R]ir£R2 


(2.16) 


£  =  I(I  du  dv  _  v 

R8  2Kr  d6  dr  r 


(2.9) 


'z9 


1  .dv  dw 

2  dz  d0J 


(2.10) 


Rz 


1  du  dw 
2{dz  dr 


(2.11) 


w  (Z  DIRECTION) 


FIGURE  2-2  FREE  BODY  DIAGRAM  AND  STATE  OF  STRAINS 
FOR  HOLLOW  CYLINDER  SPECIMENS 
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FIGURE  2- 


(a)  TANGENTIAL  DISPLACEMENT,  v 


3 


O 

<c 


o~ 

on 


3  TANGENTIAL  (a)  AND  RADIAL  (b)  DISPLACEMENT  OF  HOLLOW 
CYLINDER  SPECIMENS  UNDER  LOAD  AS  A  FUNCTION  OF 
INITIAL  RADIUS 
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Tangential  Strain  Assumption 

Note  that  Equations  2.8  through  2.11  are  identical  to  the  expres¬ 
sions  developed  for  the  resonant-column  testing  series.  The  assumption 
implicit  in  Equations  2.15  and  2.16  and  in  Figure  2-3  is  that  the  thin- 
walled  hollow  cylinder  specimen  will  strain  in  such  a  way  that  the  aver¬ 
age  radius,  RAVG»  remains  unchanged  during  loading.  This  assumption  is 
equivalent  to  the  assumption  that  Eg  =  0.  IfyRg*Y,g*0»  such  as  during 
consolidation  or  vertical  loading  alone,  a  plane  strain  condition  exists. 

Because  of  the  geometry  and  the  end  platen  boundary  conditions  of 
the  thin-walled  hollow  cylinder  specimen,  the  assumption  of  £g  =  0  is 
reasonable.  AS  the  specimen  is  loaded  compressively  in  the  vertical  di¬ 
rection,  the  specimen  will  compress  vertically  and  bulge  radially.  This 
radial  straining  will  appear  as  a  bulging  on  both  the  inner  and  outer 
walls  of  the  cylinder.  At  some  radius  within  the  specimen  the  radial 
(and  tangential)  displacement  will  be  zero.  If  that  radius  is  the  aver¬ 
age  radius,  R  ,  then  the  net  tangential  strain  will  be  zero,  as  illus- 
AVG 

trated  in  Figure  2-3 (a).  If  the  no-strain  radius  is  some  other  value,  the 
net  tangential  strain  will  not  be  zero,  but  will  be  very  small  when  com¬ 
pared  with  the  net  radial  strain.  Marachi  et  al,  (1969)  have  shown  that 
under  this  type  of  net  strain  condition  the  stress-strain  response  is 
essentially  that  of  plane  strain. 

For  this  testing  series  the  following  boundary  conditions  exist: 

1.  w  ■  a  •  z,  a  =  e  (2.12) 

z 

2.  v  *  r  •  z  •  A6q 

Note  that  these  are  identical  to  the  boundary  conditions  for  the  triaxial 
resonant-column  testing  series. 
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Strain  Tensor 

Applying  these  boundary  conditions  to  the  six  strain  equations  de¬ 
fining  the  strain  tensor  (Equations  2.8  through  2.11,  2.15,  and  2.16) 
gives  the  following  expression  for  the  strain  tensor: 


1 W 


0 

0 


0 

0 

2RAVGA60 


2*AVGA60 


z  J 


(2.17) 


A  detailed  derivation  of  the  above  expression  appears  in  Appendix 
C-2.  It  was  assumed  that  e^  =  0  for  this  derivation. 

Incremental  Principal  Strains 

Incremental  principal  strains  may  be  evaluated  using  the  strain 
tensors  developed  in  Equations  2.14  and  2.17. 

Resonant  Column 

In  the  case  of  the  triaxial  resonant  column  specimens,  whose  strain 

tensor  is  given  in  Equation  2.14,  the  Mohr's  circle  in  strain  for  the  0z 

"plane"  is  shown  in  Figure  2-4 (a) .  The  condition  illustrated  in  this 

Ae 

figure  is  the  case  in  which  the  maximum  vertical  normal  strain,  — is 

ay9z 

occurring  simultaneously  with  the  average  shear  strain,  — - — .  As  re- 

4 

corded, AYgz  and  A£z  are  "peak-to-peak"  or  double  amplitude  values.  Thus 
the  single  amplitude  values  are  Ay  a  /2  and  Ae  /2.  Because  the  compres- 

02  Z 

sional  and  torsional  strains  were  out  of  phase  and  at  different  fre¬ 


quencies,  the  average  value  of  AYa  occurring  simultaneously  with  Ae 

oz  z 


FIGURE  2-4  MOHR'S  CIRCLE  IN  STRAIN  FOR  THE  "8 z  PLANE" 

FOR  RESONANT  COLUMN  (a)  AND  HOLLOW  CYLINDER  (b) 
SPECIMENS  UNDER  DYNAMIC  LOADING  CONDITIONS 
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was  approximated  as  AYq^/2.  Thus  AYgz/4  is  used  as  the  average  single 
amplitude  value.  Because  the  vertical  axis  of  Mohr's  circle  is  equal  to 
half  the  shear  strain  in  any  direction,  the  plotted  value  of  shear  strain 


on  the  0z  "plane"  is 


ay 


9z 


8 


From  Figure  2-4(a),  and  the  strain  tensor  (Equation  2.14)  the  prin¬ 
cipal  strains  may  be  expressed  as  follows: 


.  a  Ay8z 

A£j_  *  —  +  — g—  *  tanty 


(2.20) 


Ae, 


Ae 


Ae„  * 


=  -  u 


(2.21) 


Ae„  AYq 

Ae3  =  - (U  +  __  .  tanij^) 


(2.22) 


where. 


angle  of  orientation  of  the  principal  strains 


2  tan"1 


ay, 


0z 


Ae 

8  (~r~  -  Ae  ) 
2  ct 


(2.23) 


and. 


Ae 


ct 


(2.24) 


Hollow  Cylinder 


For  the  case  of  the  thin-walled  hollow  cylinder  testing  series, 
whose  strain  tensor  is  presented  in  Equation  2.17,  the  Mohr's  circle  of 
strain  for  the  9z  "plane"  is  shown  in  Figure  2-4 (b).  In  this  figure  the 
condition  of  maximum  vertical  normal  strain  and  simultaneous  maximum 
shear  strain  is  illustrated.  In  the  hollow  cylinder  series  the  measured 
value  of  AYq2  is  the  double  amplitude  value,  so  AYqz/2  is  the  single 
amplitude  value.  However,  the  compressional  and  torsional  strains  had 


the  same  frequency.  If  they  were  also  in  phase,  as  depicted  in  Figure 

2-4 (b) ,  then  the  maximum  value  of  AyQ  would  always  occur  when  Ac  is 

oz  z 

maximum.  Thus  it  is  not  necessary  to  divide  Ay0  by  an  additional  factor 

oz 

of  2  to  get  the  average  value  of  the  shear  strain  accompanying  the  maxi¬ 
mum  value  of  vertical  normal  strain,  as  was  done  for  the  resonant  column 
specimens. 

Because  the  vertical  axis  of  the  Mohr's  circle  in  strain  is  equal  to 

one-half  the  shear  strain  in  any  direction,  the  plotted  value  of  shear 

AYqz 

strain  in  the  9z  "plane"  in  Figure  2.4(b)  is  — - — .  From  the  figure  the 
following  equations  may  be  written: 


*£z  ,  AY6z 
2  4 


(2.25) 


(2.26) 


ip  =  ~  tan  1 


4(-r-  ^ 


(2.27) 


(2.28) 


(2.29) 


Non-Combined  Loading 


For  the  special  case  of  vertical  loading  alone.  Equations  2.20 


through  2.22  reduce  to  the  following: 


(2.30) 


and  Equations  2.25  through  2.28  become: 


As 


1 


(2.32) 


(2.33) 

(2.34) 


and 


i p  =  0 


(2.35) 


In  the  special  loading  case  of  torsional  loading  alone  Equations 
2.20  through  2.23  become: 


As. 


'0* 


(2.36) 


Ae„  =  0 


(2.37) 


ay 


As. 


9z 


(2.38) 


ip  =  45' 


(2.39) 


In  the  absence  of  vertical  loading  it  is  not  necessary  to  divide 

Aya  by  two  to  get  the  average  value  associated  with  As  .  Thus  the  de- 
oz  z 

nominator  in  Equations  2.36  and  2.38  is  4  rather  than  8. 

Although  torsional  loading  alone  was  not  performed  in  the  hollow 
cylinder  test  series,  the  applicable  equations  can  be  derived  from  Equa¬ 
tions  2.25  through  2.29. 

Because  Asz/2  =  0,  Equation  2.27  now  defines  Asz>  which  is  zero,  and 


where  is  the  second  strain  invarient,  and  may  be  written: 

I,  *  -<£  £  +  £  £  +  £  £  )  +  i-y  2  +  xY  2  +  tY  2 
2  yz  zx  xy  4  yz  4  zx  4  xy 


•(£9£z  +  £z£R  +  £9£R)  +  4*0z  +  Kr  +  bsl 


-(e ie2  +  £2£3  +  £1£3} 


(2.43) 


In  terms  of  the  principal  strains.  Equation  2.42  may  be  written  as  follows 


'OCT 


3  VU1  ‘  £3>*  *  <el 


ej  2  +  uZ  -  £,) 2 


(2.44) 


This  shearing  component  has  the  effect  of  causing  shearing  or  dis- 
torsional  deformation  of  a  material,  whereas  the  octahedral  normal  strain 
has  the  effect  of  causing  purely  normal  or  compressive  deformations. 

When  considering  the  process  of  degradation  of  modulus  with  in- 
creasing  strain  amplitude,  it  is  clearly  the  shearing  strain  component  of 
the  state  of  strain  which  causes  the  degradation.  The  normal  component 
of  the  strain  tends  to  increase  the  modulus  with  increasing  strain  ampli¬ 
tude. 

For  soils  assumed  to  be  isotropic  the  octahedral  shearing  strain, 

Y ,  would  appear  to  be  a  very  useful  tool  in  evaluating  the  effects  of 

OCT 

strain  on  modulus  for  these  testing  series  since  it  will  bring  into  con¬ 
sideration  both  the  effect  of  the  dynamic  shear  strains,  AYqz,  and  the 
shear  strain  component  of  the  anisotropic  dynamic  vertical  loading,  re¬ 
presented  by  the  diameter  of  the  Mohr's  circle  in  strain  in  Figure  2-4. 
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Chapter  3 

The  State  of  Stresses 

Introduction 

It  is  the  purpose  of  this  chapter  to  describe  the  complete  state  of 
stresses  of  the  various  soil  samples  used  in  both  the  triaxial  resonant- 
column  and  the  thin-walled  hollow  cylinder  testinq  series.  This  presen¬ 
tation  will  include  the  magnitude  and  direction  of  the  principal  stresses 
throughout  the  load  history  of  the  specimens. 

Assuming  that  the  material  is  behaving  both  linearly  and  elastically, 
the  principal  stresses  and  directions  may  be  easily  calculated  if  the 
values  of  all  stresses  on  three  mutually  perpendicular  planes  are  known. 

In  the  two  testing  series  conducted,  the  stresses  were  measured  in  dif¬ 
ferent  manners  and  the  specimens  were  subjected  to  different  pre-loading 
consolidation  stress  histories;  therefore,  the  state  of  stress  will  be 
discussed  separately  for  each  series. 

Resonant-Column  Testing 

Test  Procedure 

In  this  testing  series  cylindrically  shaped  specimens  of  sand  were 
isotropically  consolidated  to  one  of  the  following  confining  pressure: 

0.5  KSC,  2.0  KSC,  and  3.5  KSC  (1  KSC  *  1  kg/cm2  =  98.07  kN/m2) .  The  sam¬ 
ples  were  then  fixed  at  the  base  and  excited  from  the  cap  with  a  pure 
sinusoidal  force  of  very  low  amplitude.  The  frequency  of  this  excita¬ 
tion  was  adjusted  to  the  undamped  natural  frequency  (also  called  the 
resonant  frequency)  of  the  soil  column.  The  acceleration  response  of  the 


soil  column  was  measured  with  a  calibrated  accelerometer  mounted  in  the 
cap  of  the  specimen. 

The  resonant  frequency  was  obtained  with  the  use  of  am  oscilloscope 
by  displaying  a  Lissajous  figure  of  the  excitation  and  the  acceleration 
response.  When  the  two  axes  of  the  elliptical  Lissajous  figure  are  in 
the  same  directions  as  the  input  axes  (i.e.:  the  vertical  and  horizontal 
plates) ,  then  the  two  signals  are  at  the  same  frequency  but  90°  out  of 
phase.  For  the  resonamt-column  sample,  the  frequency  at  which  the  input 
excitation  and  the  acceleration  response  are  90°  out  of  phase  is  by  de¬ 
finition  the  undamped  natural  frequency. 

From  the  calibrated  acceleration  response  and  the  geometry  of  the 
specimen,  the  peak-to-peak  strain  response  may  be  easily  calculated. 

During  testing  of  a  specimen  the  excitation  may  be  axial  (vertical) , 
torsional,  or  both  simultaneously.  In  a  typical  test  the  specimen  will 
be  loaded  with  a  particular  excitation  for  less  than  two  minutes,  then 
the  excitation  amplitude  is  increased  and  the  process  repeated. 

Values  of  dynamic  unconstrained  compression  modulus  (also  called 
dynamic  Young's  modulus)  and  dynamic  shear  modulus  are  calculated  from  the 
resonant  frequency,  the  sample  geometry  and  weight  and  the  dynamic  re¬ 
sponse  characteristics  of  the  testing  apparatus.  A  more  detailed  dis¬ 
cussion  of  this  calculation  is  included  in  Appendix  C-3. 

Stress  Equations 

The  net,  peak-to-peak  vertical  and  torsional  stress  may  be  calcu¬ 
lated  by  multiplying  the  calculated  modulus  by  the  calculated  peak-to- 


peak  strain  as  follows: 
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and 


A  a 


E  •  A  £ 


At 


zfl 


G  •  Ay 


z9 


where : 

E 

G 

Ao 

2 

Ae 


At 


z9 


and  Ay 


z9 


unconstrained  compression  modulus, 
shear  modulus, 

peak-to-peak  vertical  stress, 
peak-to-peak  vertical  strain, 
peak-to-peak  torsional  stress, 
peak-to-peak  torsional  strain. 


Vertical  Loading  Alone 


(3.1) 

(3.2) 


The  typical  loading  sequence  for  the  case  of  vertical  loading  alone 
is  shown  in  Figure  3-1.  In  this  figure  it  can  be  seen  that  the  calculation 
of  the  magnitude  and  direction  of  the  principal  stresses  is  greatly  sim¬ 
plified.  The  principal  stress  directions  are  coincident  with  the  di¬ 
rections  of  the  cylindrical  coordinate  system  shown  in  Figure  2-1;  thus 
their  values  are  equal  to  the  values  of  0,0,  and  aa  as  a  function  of 

Z  R 

time ,  where , 

o  =  vertical  normal  stress, 

z 

o  =  radial  normal  stress, 

K 

=  tangential  normal  stress. 


The  mathematical  expression  for  the  vertical  principal 


is  as  follows: 


a  -  o 


Aa 


lc 


sin(u)  •  t) 
z 


where  a  =  vertical  principal  stress  during  consolidation 


(3.3) 
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and  =  2  •  tt  -  f  (3.4) 

with  f  as  the  vertical  resonant  frequency. 

If  a  free-body  is  isolated  as  was  done  in  Figure  2-1,  and  the  maxi¬ 
mum  vertical  stress  condition  is  plotted  upon  it,  the  illustration  in 
Figure  3-2 (a)  is  obtained.  Figure  3-2 (b)  shows  the  Mohr's  circle  plot 
for  this  loading  condition.  Note  from  this  figure  that: 

A  a 


a 


1 


=  a 

z 


(3.5) 


a 


2 


(3.6) 


where , 


0^  =  major  principal  stress, 

a2  =  intermediate  principal  stress, 
=  minor  principal  stress. 


Torsional  Loading  Alone 


Under  the  condition  of  torsional  stress  alone,  the  maximum  principal 
stress  Mohr's  circle  is  shown  in  Figure  3-3(a).  With  these  loading  con¬ 
ditions,  the  principal  stresses  in  the  Qz  "plane"  are  oriented  45°  from 
the  axes  in  cylindrical  coordinates.  These  principal  stress  values  may 
be  expressed  as  follows: 


a 


a 

a 


1 

2 

3 


lc 


'2c 


3c 


sin  (u)_ 
T 


sin (wT 


“] 

«] 


(3.7) 


(3.8) 

(3.9) 


where  ,  ^2c'  and  a3C  are  the  principal  stresses  during  consolidation. 
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(a) 


(b) 


FIGURE  3-3  MOHR'S  CIRCLE  IN  STRESS  (a)  AND  TWO-DIMENSIONAL 
FREE  BODY  DIAGRAM  (b)  FOR  RESONANT  COLUMN 
SPECIMENS  UNDER  CONDITION  OF  MAXIMUM 
INSTANTANEOUS  TORSIONAL  STRESS 
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2a  .  -  a, 

ct  i 


j  tan"1 


At 


0z 


sin (p  •  uj  •  t) 


Aa 

z  ■  , 

— r-  •  sin(a)  •  t) 
z  z 


(3.15) 

(3.16) 


(3.17) 


fAa  *] 

***  °ct  *  aR  +j^~4~  ’  sin(wz  •  t)J  (3.18) 

where  a  is  the  center  of  the  9z  Mohr's  circle, 
ct 

A  more  detailed  derivation  of  the  above  expressions  appears  in 
Appendix  C-4. 


Hollow  Cylinder  Testing 
Sample  Preparation 

these  testing  series  thin  walled  hollow  cylinder  specimens  of 
sand  were  anisotropically  consolidated  with  a  principal  stress  ratio  of 
0.54.  The  minor  principal  stress  during  consolidation  was  set  at  one  of 
the  following  values:  0.5  KSC,  2.0  KSC,  or  3.5  KSC. 

Consolidation  Stress  Equations 

As  discussed  in  Chapter  2,  the  load  response  of  the  hollow  cylinder 
specimen  is  closely  analogous  to  plane  strain  type  loading  in  the  zR 
plane  during  consolidation.  If  we  apply  the  expressions  for  Hooke's  Law 
in  the  special  case  of  plane  strain  loading  (Equations  2.2  through  2.4) 
to  the  present  example,  we  obtain  the  following  expression  for  the  in¬ 
termediate  prinaipal  stress  during  consolidation: 
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°2c 

— 

U‘°lc  *  °3=> 

(3.19) 

where: 

aic 

= 

major  principal  stress 

during 

consolidation 

a2c 
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intermediate  principal 

stress 

during  consolidation 

a3c 

= 

minor  principal  stress 

during 

consolidation 

and 

U 

— 

Poisson's  ratio 

Because  the  directions  of  the  principal  stresses  are  coincident  with  the 
directions  of  the  axes  in  cylindrical  coordinates,  and  because  the  princi¬ 
pal  stress  ratio  during  consolidation  is  known,  the  consolidation  stress 
state  for  the  hollow  cylinder  test  may  be  stated  as  follows: 


lc 
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1.85  •  a, 
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(3.20) 
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ae  = 

2.85  •  y  •  cr3c 

(3.21) 

3c 

°R 

(3.22) 

It  is  noteworthy  that,  should  the  Poisson's  ratio  be  less  than  0.35,  then 
the  tangential  principal  stress  would  actually  be  the  minor  principal 
stress  and  the  radial  stress  would  be  the  intermediate  principal  stress. 
This  consolidation  stress  state  is  illustrated  graphically  in  Figure  3-5. 

Test  Procedure 

Following  consolidation,  the  hollow  cylinder  specimens  were  excited 
at  a  frequency  of  approximately  0.33  hertz  with  a  moderate  amplitude 
sinusoidal  force.  The  loading  excitation  was  applied  in  the  vertical  di¬ 
rection,  the  torsional  direction,  or  in  both  the  vertical  and  torsional 
directions  simultaneously. 

During  the  excitation,  which  was  applied  for  10  to  20  cycles,  the 
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°2c  CTlc 
(c) 

FIGURE  3-5  FREE  BODY  DIAGRAM  (a)  AND  MOHR'S  CIRCLES  IN 
STRESS  (b  AND  c)  FOR  HOLLOW  CYLINDER  SPECIMENS 
UNDER  CONDITION  OF  ANISOTROPIC  CONSOLIDATION 


vertical  stress  and  strain  were  recorded  directly  on  an  XY  recorder,  and 
the  vertical  stress,  torsional  stress,  and  torsional  strain  were  plotted 
as  a  function  of  time  on  a  strip  chart  recorder.  Several  typical  load 
records  are  shown  in  Figure  3-6.  After  loading,  the  cyclic  stress  levels 
were  adjusted,  and  the  process  repeated. 

Vertical  Loading  Alone 

The  typical  loading  sequence  for  the  case  of  vertical  excitation 
alone  is  illustrated  in  Figure  3-7 (a).  In  this  figure  we  can  observe 
several  differences  from  the  case  of  the  resonant-column  testing  series 
that  was  shown  in  Figure  3-1.  The  peak  stresses  and  strains  occur  es¬ 
sentially  at  the  same  time,  but  the  zero  crossings  indicate  a  phase  lag 
of  strain  behind  stress.  This  is  further  illustrated  in  Figure  3-7 (b), 
as  a  hysteretical  stress-strain  response.  This  hysteresis  represents  am 
energy  loss  during  loading  and  is  a  measure  of  damping  during  the  loading 
sequence.  Although  some  non-elastic  response  occurs  in  all  loading,  it 
is  relatively  insignificant  at  the  very  low  levels  of  loading  encountered 
during  the  resonant-column  testing  series.  At  the  higher  strains  seen  in 
the  hollow  cylinder  testing  series,  however,  hysteretic  stress-strain  re¬ 
sponse  is  common.  This  damping  effect  will  be  discussed  in  more  detail 
in  Chapter  4. 

A  second  peculiarity  of  this  testing  series  is  the  anisotropic  con¬ 
solidation  which  the  samples  have  undergone.  This  has  the  primary  ef¬ 
fect  of  moving  the  dynamic  loading  effects  higher  up  the  stress-strain 
response  curve  to  a  region  where  higher  total  strains  (and  the  resultant 
non-linear,  inelastic  load  response)  are  seen.  In  this  region  of  the 
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material  response,  a  significant  level  of  load  resisting  stress  has  been 
mobilized  in  the  specimen. 

In  addition,  because  of  the  nearly  plane  strain  type  load  response 
of  the  hollow  cylinder  specimen,  the  tangential  normal  stress  varies  with 
the  vertical  normal  stress  and  the  Poisson's  ratio.  This  effect  is  il¬ 
lustrated  schematically  in  a  Mohr's  circle  diagram  in  the  9z  "plane"  in 
Figure  3-8. 


The  mathematical  expressions  for  the  principal  stresses  under  this 
loading  condition  are  as  follows: 
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where  A  is  a  constant.  These  equations  further  reduce 

to  the  following: 
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Torsional  Loading  Alone 

Under  the  condition  of  torsional  stress  alone,  the  maximum  principal 
stress  Mohr's  circle  in  the  0z  plane  is  shown  in  Figure  3-9(a).  Unlike 
the  case  with  the  resonant-column  testing  series,  the  principal  stresses 
in  the  "9z"  plane  are  not  oriented  at  an  angle  of  45°  from  the  cylindri¬ 
cal  coordinate  axes.  Indeed,  the  angle  of  orientation,  B,  varies  both 
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with  Poisson's  ratio  and  with  time.  This  orientation  effect  is  illustrated 
in  Figure  3-9  (b),  which  is  a  two-dimensional  free-body  diagram  in  the  Sz 
"plane". 

Under  these  conditions  of  loading,  the  state  of  stresses  is  defined 
by  the  following  equations: 
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Combined  Loading 


In  the  case  of  simultaneous  vertical  and  torsional  excitation,  the 
loading  sequence  is  shown  in  Figure  3-10.  The  greatest  complication  in 
this  stress  state  occurs  because,  although  the  two  excitations  are  at  the 
same  frequency,  they  may  be  at  different  phases.  This  condition  is  shown 
in  Figure  3-10 (a).  The  maximum  principal  stress  Mohr’s  circle  in  the  0z 
"plane"  is  shown  in  Figure  3-10 (b).  The  following  equations  define  the 
state  of  stresses  under  these  loading  cond  cions: 
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FIGURE  3-10  STRESS  TIME  HISTORIES  (a)  AND  MOHR'S  CIRCLE  IN 
STRESS  (b)  FOR  SIMULTANEOUS  DYNAMIC  LOADING  OF 
HOLLOW  CYLINDER  SPECIMENS  UNDER  CONDITION  OF 
MAXIMUM  INSTANTANEOUS  VERTICAL  AND  TORSIONAL  STRESS 
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a.  =  1.85  •  a,  (3-20) 

lc  3c 

and  a  is  the  phase  lag  of  torsional  loading  to  vertical  loading. 

A  more  detailed  derivation  of  these  equations  appears  in  Appendix 

C-S. 
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Chapter  4 

Soil  Moduli  and  Damping  Factors 

Introduction 

In  this  chapter  the  strain  tensor  developed  in  Chapter  2  and  the 
state  of  stresses  described  in  Chapter  3  will  be  further  developed  and 
combined  to  produce  a  complete  picture  of  the  stress-strain  behavior  of 
triaxial  resonant-column  and  thin-walled  hollow  cylinder  test  specimens 
under  dynamic  loading  conditions. 

Resonant  Column  Testing 

In  the  triaxial  resonant-column  testing  series,  the  moduli  are 
determined  directly  in  the  manner  described  in  Chapter  3.  Damping 
factors  are  also  obtained  during  this  testing  series  using  either  the 
amplitude  decay  technique  or  the  logarithmic  decrement  approach. 

Dynamic  Interaction  Theory 

When  evaluating  the  stress-strain  response  of  these  specimens 
under  combined  vertical  and  torsional  loading,  it  is  necessary  to 
consider  the  effects  of  the  reorientation  of  the  principal  stress 
directions  during  loading  to  evaluate  the  interaction. 

As  described  in  Chapter  3,  the  principal  stresses  and  their 
directions  will  change  throughout  each  combined  loading  test  because 
the  two  excitations  are  at  two  independent  frequencies.  As  a  test  is 
conducted,  any  interaction  effects  resulting  from  the  reorientation  of 
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the  principal  stress  directions  will  vary  continuously,  and  the 
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Lissajous  figure  observed  on  the  oscilloscope  will  be  an  integration 
of  the  combined  excitation  and  response  over  time.  The  modulus  and 
damping  factor  calculated  from  this  Lissajous  figure  will  be  "weighted 
average"  values,  averaging  any  interaction  effects. 

Because  it  is  being  assumed  that  the  specimen  is  behaving  as  a 
linear,  elastic,  isotropic  material,  the  direction  of  incremental 
strains  during  loading  will  coincide  with  the  direction  of  incremental 
stresses,  and  the  principal  strain  and  stress  directions  will  coincide. 
Consequently,  the  path  of  travel  of  the  dynamic  waves  will  be  distorted 
to  follow  the  direction  of  the  principal  stresses,  and  the  effective 
velocity  of  propagation  through  the  material  will  be  decreased.  This 
effect  will  evidence  as  a  decrease  in  the  resonant  frequency  of  the 
specimen,  and  a  decrease  in  the  effective  modulus  of  the  material.  This 
principle  is  illustrated  graphically  in  Figure  4-1. 

Data  Reduction 

A  computer  program  has  been  developed  which  reduces  the  raw  data 
from  the  resonant-column  testing  series  to  produce  values  of  the 
moduli  and  damping  factors  under  combined  vertical  and  torsional  dynamic 
loading  conditions ,  and  to  calculate  the  values  and  directions  of  the 
maximum  principal  stresses  during  loading.  This  computer  program, 
called  Program  RC,  is  included  as  Appendix  D-l. 

Example  results  of  this  testing  series,  as  reduced  by  Program  RC 
are  presented  in  Appendices  B-l  through  B-3,  and  summaries  of  the  re¬ 
sults  are  given  in  Chapter  5.  The  complete  results  are  available  in 
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FIGURE  4-1  EFFECT  OF  THE  DIRECTION  OF  THE  INCREMENTAL  STRESSES 
AND  STRAINS  ON  PROPAGATION  OP  DISPLACEMENTS  IN 
CYLINDRICAL  RESONANT  COLUMN  SAMPLES 


Hollow  Cylinder  Testin' 


In  the  thin-walled  hollow  cylinder  testing  series,  the  vertical 
and  torsional  stresses  and  strains  were  measured  directly  as  described 
in  Chapter  3.  Moduli  and  damping  factors  may  be  calculated  from  this 
data  provided  the  complete  state  of  stresses  and  strains  are  known. 

Consolidation  Procedure 

During  each  test,  a  test  specimen  was  first  isotropically 
consolidated  to  0.5  KSC  confining  pressure  by  intercell  vacuum.  The 
specimen  was  then  loaded  vertically  with  an  additional  compressive 
stress  to  provide  an  anisotropic  consolidation  ratio, 

This  occurred  when  the  total  vertical  compressive  stress  was  0.925  KSC. 
While  several  tests  were  conducted  upon  specimens  in  this  consolidation 
state,  a  number  of  tests  received  further  consolidation  preparation 
before  dynamic  testing  was  performed.  As  discussed  earlier,  tests 
were  performed  upon  samples  with  consolidation  lateral  confining 
pressures  of  0.5,  2.0,  and  3.5  KSC. 

For  those  tests  conducted  at  the  two  higher  confining  pressures, 
the  above  consolidation  process  was  continued  in  the  same  sequence  as 
described  above.  First  the  effective  cell  pressure  was  increased  from 
0.5  KSC  to  2.0  KSC,  then  the  additional  vertical  compressive  stress 
was  increased  from  0.425  KSC  to  the  value  causing  an  anisotropic 
consolidation  ratio  of  0.54.  This  additional  vertical  stress  value 
was  1.7  KSC,  resulting  in  a  total  vertical  stress  of  3.7  KSC. 


Similarly,  for  specimens  to  be  tested  at  a  lateral  confining 
pressure  of  3.5  KSC,  the  effective  cell  pressure  was  increased  to  that 
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value,  then  the  additional  vertical  stress  was  increased  from  1.7  KSC 
to  2.975  KSC  to  create  an  anisotropic  consolidation  ratio  of  0.54. 

Consolidation  Strains 

During  this  consolidation  process,  the  vertical  displacement 
proximeter  was  installed  immediately  following  the  initial  consolidation 
steps.  It  was  therefore  possible  to  record  the  vertical  strain,  e  , 

Z 

during  the  remaining  consolidation  steps  where  lateral  confining 
pressures  greater  than  0.5  KSC  were  used.  This  record  was  made  on 
the  XY  recorder  for  a  number  of  tests,  and  the  summary  of  these 
records  are  shown  in  Figures  4-2  through  4-7. 

In  Figures  4-2  and  4-3,  the  variation  in  vertical  strain,  E  ,  is 

z 

shown  for  the  increment  of  consolidation  loading  resulting  from  the 
change  in  intercell  vacuum  only.  This  change  in  stress  corresponds 
to  an  isotropic  increment  of  loading  applied  to  a  specimen  which  is 
already  at  some  anisotropic  state  of  stress.  Since  the  XY  recorder 
used  to  measure  this  strain  increment  was  plotting  it  against  the 
additional  vertical  stress  change,  which  was  zero  for  this  increment 
of  loading,  it  was  not  possible  to  show  the  true  shape  of  the  stress- 
strain  response.  The  stress-strain  response  is  therefore  shcvn  as  a 
dashed  line  between  the  known  values. 

Figure  4-4  provides  a  summary  of  the  data  presented  in  Figures 
4-2  and  4-3,  and  an  extrapolation  to  allow  for  predicting  the  total 
consolidation  strains.  The  predicted  isotropic  increments  of  vertical 
strain  shown  in  Figure  4-4  are  as  follows: 
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FIGURE  4-4  SUMMARY  OF  VARIATIONS  OF  ez  WITH  ISOTROPIC  INCREMENTS 
OF  LATERAL  CONSOLIDATION  LOADING,  o3;  AND  EXTRAPOLATION 
OF  INITIAL  INCREMENT  OF  STRAIN  FROM  0  TO  0.5  KSC 
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£Z(!J3c  =  0,5  KSC)  ~  °-10%  (4.1) 

£  (03c  =  2.0  KSC)  =  0.37%  (4.2) 

and  £z<a3c  =  3*5  KSC)  '  °-65%  (4.3) 

The  short  extrapolated  line  in  this  figure  is  an  estimate  of 
the  amount  of  vertical  straining  resulting  from  the  initial  application 
of  0.5  KSC  of  intercell  vacuum  upon  the  sample.  This  value  could  not 
be  directly  measured  as  the  sample  was  still  in  the  forming  mold  at 
the  time  of  application  of  the  intercell  vacuum. 

It  should  be  emphasized  that  the  stress-strain  response  approxi¬ 
mated  in  Figure  4-4  is  not  truly  representative  of  the  stress-strain 
response  of  the  samples  to  isotropic  loading.  Because  all  samples 
were  in  an  anisotropic  state  of  stress  before  the  various  isotropic 
loading  increments  were  applied,  the  lateral  stresses  were  increasing 
disproportionately  with  the  vertical  stress  during  the  loading  increment. 
The  primary  purpose  of  plotting  this  figure  was  to  allow  an  extra¬ 
polation  of  the  initial  vertical  strain  increment  value  and  the  estimation 
of  total  consolidation  strains. 

For  the  anisotropic  consolidation  increments  of  stress,  the 
vertical  compressive  stress  was  adjusted  as  described  earlier  to  provide 
an  anisotropic  consolidation  ratio  of  0.54.  Since  both  the  change  in 

vertical  stress,  Ac  ,  and  the  change  in  vertical  strain,  Ac  ,  were 
z  z 

plotted  on  the  XY  recorder  for  this  stress-strain  increment,  the  actual 
shape  of  the  stress-strain  responst  is  known. 

In  Figures  4-5  and  4-6  the  anisotropic  stress-strain  increments 
are  shown  combined  with  the  isotropic  stress  increments.  The  solid  lines 


FIGURE  4-6  VARIATION  OF  ez  WITH  BOTH  ISOTROPIC  AND  ANISOTROPIC 
INCREMENTS  OF  VERTICAL  CONSOLIDATION  LOADING,  oz , 


FROM  3.7  TO  6.475  KSC,  FOR  SEVERAL  TESTS 
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shown  in  these  curves  correspond  to  the  actual  anisotropic  stress- 
strain  response  of  the  soil  recorded  by  the  XY  recorder.  The  data 
on  these  two  figures  are  summarized  in  Figure  4-7,  along  with  an  extra¬ 
polation  to  estimate  the  total  consolidation  strain  values.  The 
extrapolated  line  on  Figure  4-7  includes  the  earlier  isotropic  increment 
extrapolation  from  Figure  4-4,  and  an  anisotropic  increment  extrapolation 
deduced  from  the  anisotropic  loading  increment  data  shown  in  the  figure. 
The  total  consolidation  vertical  strains  may  now  be  estimated  as  follows: 
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It  is  interesting  to  note  that  for  all  of  the  tests  shown,  which 
represent  a  wide  range  of  densities,  the  recorded  values  of  the  total 
vertical  strain  under  anisotropic  consolidation  conditions  fall  within 
a  few  percent  of  the  average  values  in  Equations  4.4  through  4.6. 

Consolidation  Modulus  Calculation 


From  the  strain  tensor  developed  in  Equation  2.17,  we  can  calculate 
the  volumetric  strain  as  follows: 
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and  from  Hooke's  Law  (Equation  2.4), 
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which  reduces  to: 
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Combining  Equations  4.7  and  4.10  gives: 
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which  reduces  to: 
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The  static  compression  modulus  E  was  calculated  from  Equation  4.12 
for  a  variety  of  values  of  Poisson's  ratio  and  the  results  are  shown 
in  Table  4-1.  The  similarity  in  the  calculated  moduli  for  a  specific 
Poisson's  ratio  is  a  measure  of  the  similarity  in  the  principal  stress 
ratio  during  consolidation,  and  of  linearity  of  the  stress-strain 
response. 


Tangent  and  Secant  Moduli 

During  dynamic  excitation,  the  dynamic  vertical  stress  and  strain 
will  be  superimposed  upon  the  consolidation  stress-strain  condition. 
Consequently,  it  is  possible  to  view  either  the  tangent  modulus  or 
the  secant  modulus  under  these  loading  conditions. 

Vertical  straining  during  dynamic  loading  will  be  in  accordance 
with  the  strain  tensor  developed  in  Equation  2.17  and  Hooke's  Law  as 
presented  in  Equation  2.4.  The  dynamic  tangent  compression  modulus 


56 


under  vertical  or  combined  vertical  and  torsional  dynamic  loading 
conditions  may  be  calculated  by  the  following  expression: 

A  a 

E  -  (1  -  U2)  (4.13) 

z 

A  derivation  of  this  expression  is  given  in  Appendix  C-6. 

Similarly,  Equation  4.12  may  be  rewritten  to  provide  a  means  for 
calculating  the  dynamic  secant  compression  modulus,  producing  the 
following: 


where  o  =  a,  ,  *  0,  ,  and  e  =  e,  . 

Z  lc  R  3c  z  lc 

The  secant  compression  modulus  may  be  simply  calculated  using 
Equation  4.14.  Because  this  modulus  is  an  overall  average  modulus  from 
the  preconsolidation  state,  the  values  will  not  vary  greatly  from  those 
values  in  Table  4-1  except  under  conditions  of  very  high  dynamic  straining. 

If  the  soil  sample  is  considered  representative  of  soil  conditions 
in  the  field,  with  anisotropic  consolidation,  then  the  compression 
modulus  of  interest  is  the  tangent  compression  modulus  determined  using 
Equation  4.13.  This  modulus  will  reflect  the  dynamic  stress-strain 
behavior  of  soil  which  has  reached  a  steady-state  consolidation  condition, 
and  whose  strain  will  vary  throughout  the  range  of  strains  imposed  in 


this  testing  t  ~ogram. 
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Effect  of  Torsional  Loading  on  Equations  for  Moduli 


Under  conditions  of  torsional  dynamic  loading  alone,  the  principal 
stress  and  strain  directions  will  change  during  loading  as  discussed  in 
Chapter  3.  The  Mohr's  circles  in  strain  and  stress  for  this  condition 
are  indicated  in  Figure  4-8,  with  the  reorientation  of  the  principal 
strains  in  the  6z  "plane"  indicated  by  the  angle  p.  If  the  material 
behaves  linearly  and  elastically,  then  ip  must  equal  3. 

Note  from  Figure  4-8  (a)  that  the  principal  strains  may  now  be 
written  as  follows : 
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where  £,  =  £  ,  £_ 

lc  z  2c 


£0  =  0,  and  £3c  =  eR. 


The  volumetric  strain  may  now  be  written: 
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Because  Equation  4,19  is  exactly  equal  to  Equation  4.7,  the 
expression  for  the  volumetric  strain  during  dynamic  loading  is  shown 
to  be  independent  of  reorientation  of  the  principal  stress  and  strain 
directions.  The  expressions  for  the  compression  modulus  (Equations 
4.13  and  4.14)  are  derived  from  the  expression  for  volumetric  strain. 


FIGURE  4-8  MOHR'S  CIRCLE  IN  STRAIN  (a)  AND  IN  STRESS  (b)  FOR 
HOLLOW  CYLINDER  TEST  SPECIMENS  UNDER  CONDITION  OF 
MAXIMUM  INSTANTANEOUS  TORSIONAL  STRESS  AND  STRAIN 


59 


thus  those  expressions  are  also  independent  of  reorientation  of  the 
principal  stress  and  strain  directions. 

The  dynamic  shear  modulus  under  torsional  loading  conditions  may 
be  calculated  by  the  following  expression: 


At 

Ay 


9z 

9z 


(4.20) 


Under  conditions  where  only  one  modulus  can  be  directly  calculated 
from  the  data  presented,  or  where  it  is  desired  to  predict  the  effective 
Poisson's  ratio  under  combined  vertical  and  torsional  loading  conditions, 
the  following  expression  may  be  used: 


E  =  2(1  +  p)  *  G 


(4.21) 


Data  Reduction 

A  computer  program  has  been  written  which  calculates  values  of  the 
tangent  and  secant  dynamic  compression  moduli,  the  dynamic  shear  modulus, 
and  the  magnitude  and  direction  of  the  maximum  principal  stresses  and 
strains  from  the  raw  test  data  as  a  function  of  the  Poisson's  ratio. 

This  computer  proaram  is  called  Program  HC,  and  is  included  as  Appendix 
D-2. 

Example  results  of  this  testing  series,  as  reduced  from  Program  HC, 
are  presented  in  Appendices  B-4  through  B-6,  and  complete  results  are 
summarized  in  Chapter  6.  The  complete  results  are  available  in  Griffin 


(1980)  . 
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Damping 

Introduction 

In  both  the  triaxial  resonant-column  and  the  thin-walled  hollow 
cylinder  testing  series  energy  is  dissipated  during  dynamic  loading  due 
to  damping.  Generally  speaking,  the  greater  the  strain  amplitude, 
the  greater  will  be  this  energy  loss.  With  closely  monitored  laboratory 
tests  such  as  these,  the  energy  loss  measured  is  due  solely  to  the 
inelastic  response  of  the  materials  under  test,  and  is  sometimes 
referred  to  as  hysteretic  damping  or  friction  damping.  If  the  material 
behaves  linearly  and  elastically,  the  energy  loss  and  thus  the 
damping,  will  be  zero.  As  long  as  the  darrping  is  small,  it  is  quite 
possible  (and  common  practice)  to  treat  the  peak-to-peak  stress  strain 
response  as  an  "equivalent-linear"  modulus ,  and  introduce  damping  as 
an  energy  dissipation  'over  time. 

Resonant  Column  Testing 

In  the  triaxial  resonant-column  testing  series,  the  damping  was 
measured  using  either  the  amplitude  decay  technique  or  the  multiplica¬ 
tion  factor  approach.  The  first  of  these  involves  the  excitation  of 
the  specimen  at  its  resonant  frequency  and  observing  the  acceleration 
response  as  a  function  of  time.  When  the  driving  force  is  abruptly 
stopped,  the  sample  will  continue  to  oscillate  freely,  but  the 
acceleration  amplitude  will  decay  with  time  because  of  damping.  The 
value  of  the  damping  may  be  calculated  as  follows: 


(4.22) 


where  A  =  amplitude  decay  danping 


n  =  number  of  cycles,  where  n  >  1 
aQ  =  acceleration  amplitude  of  the  O-th  cycle 
*  acceleration  amplitude  of  the  n  -  th  cycle 
This  is  illustrated  graphically  in  Figure  4-9(a). 

The  amplitude  decay  damping  is  related  to  the  ratio  of  critical 
damping,  D,  by  the  following  expression: 


A  = 


2  IT  D 


ST? 


(4.23) 


and,  for  small  values  of  damping,  D  may  be  expressed  in  percent  as 
follows : 


D  (%)  =  —  •  100% 
2tt 


(4.24) 


The  use  of  the  multiplication  factor  approach  requires  a  look  at 
the  response  spectrum  for  the  specimen  as  shown  in  Figure  4-9 (b) .  The 
resonant  frequency  is  obtained  as  described  in  Chapter  3;  then  frequencies 
and  amplitudes  may  be  determined  from  other  points  along  the  "bell 
shaped  curve,"  as  shown  in  Figure  4-9 (b). 

The  sharpness  of  the  "bell"  is  an  indication  of  the  energy 
dissipation  of  the  sample.  A  perfectly  linear,  elastic  material  would 
respond  at  one  precise  resonant  frequency,  and  would  have  a  response 
spectrum  which  is  just  a  vertical  line  at  the  resonant  frequency. 

The  values  of  frequency  and  response  shown  in  Figure  4-9 (b)  are 
convenient  because  they  lead  to  relatively  simple  expressions  for 
calculating  the  damping.  The  values  of  f^  and  f^  were  selected  such 
that  the  acceleration  response  at  those  frequencies  was  l//l  times  the 
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maximum  acceleration  response  occuring  at  or  near  resonance.  Using 
these  frequencies,  the  damping  may  be  expressed  as  follows: 


A  = 


it  (f  _  -  f  )  a 
2  1  max 


res 


(4.25) 


or 


(f. 


D  (%)  * 


V 


max 


2  f 


100% 


(4.26) 


The  value  a  was  selected  such  that  the  frequency  at  which  it 
occurred  is  equal  to  /2  times  the  resonant  frequency.  In  this  case 
the  damping  may  be  calculated  as  follows: 


\  *  P 

A  •  - 


”Aa 

a. 

amax 

1 

a 

Aa 

_  max  _ 

L  al  J 

(4.27) 


or 


D(%)  -  ^ 


"Aa 

a  1 

amax 

1 

a 

t> 

Q 

_  max  _ 

L  alj 

100% 


(4.28) 


where  Aa  =  driving  stress  at  a  response 
amax  max 

Aa^  =  driving  stress  at  a^  response 
If  the  driving  force  is  held  constant  throughout  the  response  spectrum, 
which  is  usually  the  case.  Equations  4.27  and  4.28  reduce  to  the 
following : 


7  a. 


A  = 


2a 


(4.29) 


max 


and 


D  = 


4a 


max 


100% 


(4.30) 


The  practical  advantage  to  using  Equations  4.27  and  4.28  in  the 


form  presented  is  that  the  factor 


rJi-1 

LAOaiJ 


does  not  change  appreciably 
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during  testing,  and  normally  needs  to  be  measured  only  once.  The 
damping  may  then  be  calculated  by  merely  knowing  the  acceleration 
response  and  the  driving  stress  at  resonance  throughout  testing.  The 
damping  equations  may  thus  be  rewritten  as  follows: 


A 


2tt6 

100% 


(4.31) 


and 


D (%)  =  6  • 


(4.32) 


where 


,  100%  r  ai  i 

6  =  I  4  LAcr  J  (4.33) 

I 

'  damping  calibration  factor 


Hollow  Cylinder  Testing 

In  the  thin-walled  hollow  cylinder  testing  series,  the  energy 
dissipation  is  seen  as  a  hysteresis  loop  stress-strain  response,  as 
shown  in  Figure  4-10.  The  "equivalent  linear"  modulus  is  shown  as  the 
ratio  of  peak-to-peak  stress  to  peak-to-peak  strain.  In  Figure  4-10, 
the  following  equations  may  be  written: 


A  0 

AB  = 

A '  B  1 

-  z 

=  ~T~ 

(4.34) 

Ae 

OB  = 

B'O 

z 

2 

(4.35) 

AB 

_  A'B' 

(4.36) 

‘PS  " 

OB 

B'O 

MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS  1%  ’  A 


The  damping  may  be  expressed  as  follows: 


a  +  a 
ABO  A'B’O 

where  A^  »  area  of  the  loop 

AABO  *  ****  the  triangle  ABO 

AA'B'0  *  area  the  triangle  A'B'O 

This  value  of  hysteretic  damping  is  related  to  the  ratio  of  critical 
damping  by  the  following  relationship: 


,  2  TTD 

or,  for  small  values  of  damping. 


_X_ 

27T 


100% 


(4 


D  (%) 


(4 


Low  Strain  Combined  Cyclic  Loading 

Introduction 

In  this  chapter  the  actual  laboratory  test  results  reduced  from  the 
triaxial  resonant  column  series  will  be  presented.  Shown  in  Figure  5-1 
is  a  photograph  of  the  resonant  column  testing  device  with  a  specimen  in 
place.  Listed  in  Table  5-1  is  a  summary  of  the  seventeen  tests  which 
have  been  analyzed  for  this  presentation. 

As  discussed  previously,  the  moduli  during  dynamic  combined  loading 
will  change  both  with  vertical  strain  amplitude  and  with  compression- 
shear  interaction.  The  data  will  be  presented  in  such  a  way  that  the  two 
effects  may  be  seen  both  separately  and  together.  The  results  will  also 
be  presented  as  a  function  of  the  octahedral  shearing  strain. 

Vertical  Loading  Alone 

The  dynamic  compression  modulus  is  shown  as  a  function  of  the  cyclic 

vertical  strain  under  conditions  of  vertical  loading  alone  in  Figure  5-2. 

It  should  be  noted  that  for  all  values  of  strain  referred  to  herafter, 

the  single-amplitude  values  are  used.  The  cyclic  stress  amplitudes 

shown  on  drawings,  however,  are  peak-to-peak  values.  This  Figure  includes 

data  from  a  variety  of  samples  at  various  densities  and  confining  pressures. 

Superimposed  upon  the  drawing  are  contour  lines  showing  values  of  equal 

vertical  normal  cyclic  stress ,  Ac  .  From  this  Figure  we  can  see  that 

z 

the  range  of  cyclic  vertical  normal  stress  used  in  this  testing  series 
was  from  approximately  0.02  psi  to  25  psi.  The  corresponding  vertical 
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FIGURE  5-2  DEGRADATION  OF  COMPRESSION  MODULUS  WITH  VERTICAL  STRAIN  FOR  SEVERAL 
TESTS  UNDER  CONDITION  OF  VERTICAL  LOADING  ALONE 


strain  amplitudes  were  in  the  range  of  from  approximately  10  5%  to  3X10  2%. 
It  is  interesting  to  note  the  steepness  of  the  stress  contours,  indicating 
that  relatively  small  variations  in  strain  amplitude  will  yield  relatively 
large  variations  in  moduli. 

To  eliminate  variations  in  the  modulus-strain  curves  due  to  sample 
peculiarities,  density  variations,  and  small  measurement  errors,  it  is 
desirable  to  "normalize"  these  curves  by  plotting  a  "relative  modulus"  in 
place  of  the  absolute  value  of  modulus.  This  technique  allows  a  direct 
evaluation  of  the  relative  variation  of  modulus  with  strain  amplitude. 

In  Figure  S-3  six  tests  are  normalized  to  the  value  of  the  modulus 
at  a  single  amplitude  vertical  strain  of  5X10-3%.  This  value  was  se¬ 
lected  to  allow  for  comparison  of  these  test  results  with  those  of  the 
higher  strain,  thin-walled  hollow  cylinder  testing  series  which  will  be 
presented  in  Chapter  6. 

From  Figure  5-3  it  can  be  seen  that  the  greatest  rate  of  change  of 
modulus  with  strain  amplitude  occurs  at  the  lowest  confining  pressures. 
Although  there  is  some  variation  in  this  rate  as  a  function  of  density, 
this  variation  is  small  when  compared  with  the  influence  of  the  con¬ 
fining  pressure  on  the  rate  of  variation. 

Torsional  Loading  Alone 

The  results  of  four  tests  in  which  pure  torsional  loading  was  ap¬ 
plied  to  samples  are  presented  in  Figure  5-4.  As  with  Figure  5-2,  the 
modulus  (in  this  case  the  Shear  Modulus)  is  seen  to  change  as  a  function 
of  the  strain  amplitude.  In  these  tests  the  range  of  the  cyclic  torsional 
shear  stress  was  from  approximately  .01  psi  to  15  psi,  very  similar  to 


SEVERAL  TESTS  UNDER  CONDITION  OF  VERTICAL  LOADING  ALONE 
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the  load  range  in  vertical  compression.  The  corresponding  torsional 
shear  strain  amplitudes  were  in  the  range  of  from  approximately  4X10  5% 
to  4X10"2%. 

The  shear  modulus-shear  strain  curves  were  normalized  to  the  modulus 
at  a  shear  strain  of  5X10  3%  and  presented  in  Figure  5-5.  In  this  Figure 
it  can  be  seen  that  the  relative  variation  of  shear  modulus  with  strain 
amplitude  is  somewhat  greater  than  the  variation  of  compression  modulus 
shown  in  Figure  5-3  over  a  similar  strain  ramge. 

Combined  Loading 

Under  conditions  of  simultameous  vertical  and  torsional  dynamic 
loading,  the  modulus  will  vary  both  as  a  function  of  the  vertical  strain 
amplitude  and  as  a  function  of  the  compression-shear  interaction  effects. 
As  indicated  in  Tadsle  5-1,  six  tests  performed  under  simultameous  com¬ 
bined  loading  conditions  have  been  presented  and  amalyzed.  These  tests 
represented  two  different  densities  and  three  different  confining  pres¬ 
sures. 


Interaction  Effects 

The  effects  of  this  interaction  on  two  of  these  tests  are  shown  in 
Figure  5-6.  In  this  Figure  modulus-strain  amplitude  curves  are  plotted 
for  the  two  test  specimens,  with  contours  of  equal  shear  strain  ampli¬ 
tude  superimposed  on  the  drawing.  The  curves  labeled  Yq2  <  10  4%  are 
those  corresponding  to  the  "virgin"  modulus-strain  curve  in  which  no 
measurable  compression-shear  interaction  effects  are  present.  These 
are  the  solid  lines  in  this  Figure. 


5-6  DEGRADATION  OF  COMPRESSION  MODULUS  WITH  VERTICAL  STRAIN  UNDER  SEVERAL 
LEVELS  OF  SIMULTANEOUSLY  APPLIED  SHEAR  STRAIN 


Strain  Ratio  Effects 

These  interaction  effects  can  be  isolated  from  the  other  factors 

influencing  the  modulus  value  by  plotting  a  normalized  modulus  vs.  the 

ratio  of  average  shear  strain,  yq  ,  (see  Chapter  2  for  definition)  to  the 

uz 

vertical  normal  strain,  t  .  The  modulus  should  be  normalized  to  the 

z 

modulus  at  that  vertical  strain  amplitude  on  the  "virgin"  modulus-strain 
curve,  E<<y^. 

These  normalized  plots  were  prepared  for  all  of  the  combined  loading 
tests  in  this  series,  and  are  presented  in  Figures  5-7,  5-8,  and  5-9,  for 
confining  pressures  of  0.5,  2.0,  and  3.5  KSC,  respectively.  These  normal¬ 
ized  curves  are  combined  and  summarized  for  those  tests  with  a  relative  den¬ 
sity  of  95%  in  Figure  5-10;  and  for  all  combined  loading  tests  in  Figure  5-11. 

Octahedral  Shearing  Strain  Effects 

As  discussed  in  Chapter  2,  the  octahedral  shearing  strain,  y  ,  is 

vJVv  X 

a  very  useful  tool  in  evaluating  the  degradation  of  modulus  with  strain 
amplitude.  For  these  testing  series  a  degradation  of  modulus  has  been 
observed  both  with  increasing  dynamic  axial  strain,  and  with 

increasing  dynamic  shear  strain,  Ay^.  The  octahedral  shearing  strain 
combines  both  of  these  types  of  straining,  and  should  be  very  useful  in 
studying  the  combined  degradation  effects. 

It  should  be  noted  that  if  the  material  is  homogeneous  and  isotropic, 
then  the  degradation  of  modulus  with  octahedral  shearing  strain,  y  , 

UL  X 

will  be  independent  of  how  the  strain  was  developed;  i.e.,  from  axial  or 
torsional  straining. 

In  Figure  5-12  are  shown  the  "virgin"  normalized  compression  moduli, 
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FIGURE  5-7  DEGRADATION  OF  NORMALIZED  COMPRESSION  MODULUS  WITH  RATIO  OF  SHEAR  STRAIN  TO 
VERTICAL  STRAIN  AT  =  0.5  KSC  FOR  SEVERAL  COMBINED  LOADING  TESTS 


DEGRADATION  Of  NORMALIZED  COMPRESSION  MOOtllllS  WITH  RATIO  OF  SHEAR  STRAIN  TO 
VERTICAL  STRAIN  AT  .»,  =2.0  KSC  FOR  SEVERAL  COMBINED  LOADING  TESTS 


GRADATION  OF  NORMALIZED  COMPRESSION  MODULUS  WITH  RATIO  OF  SHEAR  STRAIN  TO 

:rtical  STRAIN  AT  =  3.5  KSC  for  several  combined  loaoing  tests 


SUWARY  OF  DF GRADATION  OF  NORMALIZED  COMPRESSION  MODULUS  WITH  RATIO  OF  SHEAR 
STRAIN  TO  VERTICAL  STRAIN  FOR  A  RELATIVE  DENSITY  OF  951 


DEGRADATION  OF  NORMALIZED  COMPRESSION  MODULUS  WITH  RATIO  OF  SHEAR 
VERTICAL  STRAIN  FOR  ALL  COMBINED  LOADING  TESTS 
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E/Emax.  plotted  vs.  Y  .  Because  these  curves  were  developed  by  applying 
vertical  loading  alone  to  the  soil  specimens ,  the  octahedral  shearing 
strain  in  that  Figure  is  developed  purely  from  vertical,  radial,  and 
tangential  normal  straining.  These  normal  strains  correspond  to  the 
principal  strains  in  this  case. 

In  the  combined  loading  case,  shear  strains  introduced  in  the  9 

z 

"plane”  result  in  a  reorientation  of  the  principal  strain  directions; 
thus,  the  octahedral  shearing  strain  in  this  case  is  developed  from  the 
combination  of  vertical,  radial,  and  tangential  normal  strains,  plus  9^ 
"plane"  shear  strains . 

The  variation  of  compression  moduli  with  the  virgin  compression 

moduli,  as  a  function  of  octahedral  shearing  strain  is  shown  in  Figures 

5-13,  5-14,  and  5-15  for  confining  pressures  of  0.5,  2.0,  and  3.5  KSC, 

respectively.  These  curves  are  combined  with  the  virgin  curves  from 

Figure  5-12  to  determine  the  effects  of  shear-compression  interaction, 

and  are  presented  in  Figures  5-16,  5-17,  and  5-18  for  the  same  three 

confining  pressures.  For  the  dashed  curves  the  increased  in  Y-*—,  arise 

from  increases  in  ya  with  E  held  constant  as  shown.  For  the  solid 
OZ  Z 

curves,  the  increases  in  YqCT  arise  from  increases  in  z^  with  Ygz  held 
as  small  as  possible,  typically  less  than  10  4%  as  shown.  The  intersection 
of  the  dashed  curves  and  the  solid  curve  occurs  when  the  contribution  of 
Y0Z  to  yqct  Is  small  compared  to  the  contribution  of  sz  for  these  tests. 
Therefore  the  value  of  YqCT  at  which  the  intersection  occurs  can  be 
computed  directly  from  for  practical  purposes. 

In  these  Figures  it  is  apparent  that  the  two  types  of  curves  are 
essentially  coincident,  and  thus  any  shear-compression  interaction  effect 
which  is  not  accounted  for  by  the  use  of  the  octahedral  shearing  strain 
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is  negligible. 

The  resultant  best-fit  normalized  moduli-octahedral  shearing  strain 
curves  are  presented  in  Figure  5-19. 

Conclusions 

The  various  test  results  were  presented  in  different  ways  to  evalu¬ 
ate  the  relative  compression-shear  interaction  under  combined  loading 
conditions . 


Strain  Ratio  Effects 

From  the  results  of  the  presentation  of  the  relative  normalized  com¬ 
pression  moduli  plotted  against  the  ratio  of  average  shear  strain,  yQ  , 

uz 

to  vertical  normal  strain,  £  ,  in  Figure  5-11,  it  is  apparent  that  there 

z 

exists  a  "threshold"  strain  ratio  below  which  negligible  interaction  ef¬ 
fects  are  observed.  It  appears  that  the  minimum  threshold  strain  ratio 
in  this  testing  program  was  approximately  1,  where  the  average  shear 
strain  equaled  the  vertical  normal  strain.  This  means  that  when  the 
shear  strain,  Yqz»  was  numerically  smaller  than  the  vertical  normal 
strain  the  degradation  of  modulus  due  to  ya  was  negligible. 

HZ 

Also  apparent  is  the  fact  that  the  moduli-strain  ratio  curve  and  the 
threshold  strain  ratio  depend  primarily  upon  the  vertical  normal  strain 
amplitude.  The  effects  of  confining  pressure  and  relative  density  upon 
these  curves  appears  to  be  very  small. 

Octahedral  Shearing  Strain  Effects 


From  the  results  of  the  presentation  of  normalized  compression 
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moduli  plotted  against  the  octahedral  shearing  strain  in  Figure  5-19  it 
appears  that  the  compression-shear  interaction  is  directly  predictable  or 
"expressible"  by  use  of  the  octahedral  shearing  strain. 

The  maximum  octahedral  shearing  strain  developed  in  this  testing 
series  was  approximately  3X10  2%. 
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Chapter  6 

High  Strain  Combined  Cyclic  Loading 

Introduction 

In  this  chapter  the  results  of  the  thin-walled  hollow  cylinder  testing 
program  will  be  presented.  A  photograph  of  the  testing  apparatus  is  shown 
in  Figure  6-1.  A  summary  of  the  twenty-one  tests  which  make  up  this 
testing  series  is  presented  in  Table  6-1. 

As  with  the  low  strain  triaxial  resonant  column  testing  program,  the 
data  is  presented  in  such  a  manner  that  changes  in  the  dynamic  compression 
modulus  due  to  vertical  strain  amplitude  and  due  to  combined  loading  in¬ 
teraction  may  be  evaluated  separately  and  independently  of  one  another. 
Additionally,  the  vertical  compression  stress-strain  response  under  con¬ 
ditions  of  very  high  shear  loading  will  be  presented  and  discussed. 

Vertical  Loading  Alone 

In  Figure  6-2  the  dynamic  compression  modulus  is  shown  as  a  function 
of  the  vertical  normal  strain  for  a  value  of  the  Poisson's  ratio  of  0.34. 
The  curves  shown  are  "virgin"  curves,  with  no  simultaneously  applied  shear 
loading.  The  value  selected  for  the  Poisson's  ratio  is  near  the  center 
of  the  normal  range  of  from  approximately  0.2  to  0.5.  Variation  of  the 
Poisson's  ratio  within  that  range  will  have  little  effect  upon  these 
curves,  since  the  calculated  value  of  the  compression  modulus  is  not  very 
sensitive  to  the  Poisson's  ratio. 

Superimposed  upon  Figure  6-2  are  contours  of  equal  values  of  cyclic 
vertical  normal  stress,  Ao^.  The  range  of  Aaz  in  this  testing  series  was 
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TABLE  6-1 

Summary  of  Hollow  Cylinder  Tests 


1 

Yd  (g/cc) 

! 

V-VERTICAL 

TEST  NO. 

Dr  (%)  j 

a,  (ksc) 

3c 

MODE 

C-COMBINED 

B10-1 

1.68 

91.8 

0.5 

C 

Bll-1 

1.67 

88.8 

0.5 

C 

B12-2 

1.69 

94.9 

2.0 

c 

B13-3 

1.69 

94.9 

3.5 

c 

AB15-2 

1.71 

100.0 

2.0 

V  &  c 

B16-2 

1.69 

94.9 

2.0 

c 

B17-3 

1.69 

94.9 

3.5 

c 

B18-1 

1.70 

97.9 

0.5 

c 

B19-2 

1.68 

91.8 

2.0 

c 

B20-3 

1.67 

88.8 

3.5 

c 

B21-1 

1.67 

88.8 

0.5 

c 

B22L-1 

1.58 

59.3 

0.5 

c 

B23L-2 

1.58 

59.3 

2.0 

C 

B24L-3 

1.55 

48.7 

3.5 

c 

A25L-1 

1.59 

62.7 

0.5 

V 

A26L-2 

1.56 

52.3 

2.0 

V 

A27L-3 

1.55 

48.7 

3.5 

V 

A28-3 

1.69 

94.9 

3.5 

V 

B29-3 

1.68 

91.8 

3.5 

c 

B30L-3 

1.58 

59.3 

3.5 

c 

B31-1 

1.68 

91.8 

0.5 

c 

*  dry  density  of  soil 
D  *  relative  density 
a3c  *  lateral  confining  pressure 


from  approximately  1  psi  to  100  psi.  For  low  confining  pressure  tests 
particularly,  these  values  of  Aa^  led  to  significant  principal  stress 
ratios  during  cyclic  loading.  These  principal  stress  ratios  are  sum¬ 
marized  in  Table  6-2.  Not  all  of  the  values  in  Table  6-2  were  actually 
reached  in  the  test  specimens.  Figure  6-2  shows  the  combinations  of  o^c 
and  Aaz  actually  attained  during  testing. 

Vertical  strain  amplitudes  achieved  during  this  testing  series  were 
in  the  range  of  from  approximately  2X10""3%  to  2X10  J%. 

The  compression  moduli  shown  in  Figure  6-2  were  normalized  against 
the  moduli  at  a  vertical  normal  strain  of  5X10  3%,  as  shown  in  Figure  6-3. 
As  with  the  low  strain  resonant  column  testing  program,  the  variation  of 
relative  compression  modulus  with  strain  amplitude  appears  to  be  primarily 
a  function  of  the  confining  pressure,  with' a  lower  confining  pressure  cor¬ 
responding  to  a  greater  rate  of  degradation. 

Combined  Loading 

As  with  the  low  strain  resonant  column  testing  program,  cyclic  loading 
tests  under  conditions  of  simultaneous  vertical  and  torsional  loading  were 
performed  on  a  number  of  hollow  cylinder  samples.  The  specimens  tested 
represented  two  different  densities  and  three  different  confining  pres¬ 
sures. 


Strain  Measurements 

In  this  testing  series  the  value  of  the  relative  rotational  displace¬ 
ment  of  the  top  of  the  specimen  (and  thus  the  shear  strain)  was  measured 
directly  with  an  LVDT  mounted  on  the  side  of  the  specimen.  Although  the 


relatively  larger  measurements  were  consistent,  the  relatively  smaller 


shear  strain  measurements  were  somewhat  erratic.  This  was  found  to  be 
primarily  due  to  the  LVDT  core  friction,  which  was  significant  at  very 
small  displacements,  causing  the  LVDT  to  "hang-up"  and  under-measure  the 
true  shear  strain.  This  difficulty  was  not  overcome  by  great  care  in 
setting  up  the  LVDT  to  avoid  core  contact,  because  consolidation  and 
dynamic  straining  quickly  re-created  the  core  friction  problem. 

The  problem  was  further  complicated  by  the  fact  that  the  calculated 
shear  modulus  values  were  extremely  sensitive  to  the  very  small  shear 
strains.  A  relatively  minor  variation  in  the  measured  shear  strain  could 
result  in  a  large  error  in  the  calculated  modulus,  and  would  result  in 
unreasonable  predicted  values  for  Poisson's  ratio. 

It  was  concluded  that  a  more  accurate  estimate  of  the  actual  shear 
strain  could  be  obtained  using  the  measured  stresses  and  vertical  strain 
together  with  the  elastic  theory  formulations  developed  in  Chapters  2,  3, 
and  4. 

The  measured  values  of  strain  and  stress  provided  a  redundancy  of 
data  from  which  a  Poisson's  ratio  could  be  computed.  Because  the  shear 
strain  is  relatively  insensitive  to  the  value  of  the  Poisson's  ratio,  the 
computed  values  of  the  Poisson's  ratio  are  extremely  sensitive  to  small 
errors  in  shear  strain.  By  back-calculating  values  of  the  shear  strain 
for  various  values  of  Poisson's  ratio,  a  much  more  accurate  evaluation  of 
the  shear  strain  may  be  made. 

The  measured  values  of  the  vertical  and  shear  stress  were  obtained 
with  calibrated  differential  pressure  transducers  attached  to  the  loading 
cylinders,  and  are  relatively  accurate  and  consistent  in  determining 
those  stress  values.  The  vertical  strain  was  measured  with  a  proximeter 
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device,  which  accurately  measures  displacements  without  physical  contact, 
avoiding  the  core  friction  problem.  These  data  measurement  devices  are 
discussed  in  more  detail  in  Appendix  A-2. 

When  values  of  the  shear  strain  were  calculated  from  the  measured 
stresses  and  vertical  strain  values,  and  for  various  values  of  Poisson's 
ratio,  as  described  above,  they  were  found  to  be  in  the  same  relative 
range  as  the  measured  values,  but  were  generally  slightly  larger. 

In  the  example  of  test  results,  presented  in  Appendices  B-4  through 
B-6,  both  the  measured  and  the  calculated  values  of  the  shear  strain  are 
presented.  The  calculated  values  of  the  shear  strain  were  used  in 
developina  the  values  of  the  shear  modulus,  G,  and  the  strain  ratios 
and  octahedral  shearinq  strains  used  in  this  chapter.  The  complete  test 
results  are  reported  elsewhere  (Griffin,  1980). 

Strain  Ratio  Effects 

Normalized  plots  showing  the  degradation  from  the  "virgin"  modulus 
as  a  function  of  the  ratio  of  shear  strain  to  vertical  normal  strain  were 
prepared  for  three  different  confining  pressures.  The  resulting  plots 
are  shown  in  Figures  6-4,  6-5,  and  6-6  for  lateral  confining  pressures  of 
0.5,  2.0,  and  3.5  KSC,  respectively. 

The  plotted  modulus  reduction  curves  shown  in  Figures  6-4  through 
6-6  are  combined  and  summarized  in  Figure  6-7  for  two  values  of  vertical 
normal  strain:  10  ’%  and  10  2%. 

The  effect  of  varying  the  Poisson's  ratio  is  to  vary  the  strain  ratio 
slightly.  This  effect  is  shown  in  Figure  6-8. 


DEGRADATION  OF  NORMAL IZEO  COMPRESSION  MODULUS  WITH  RATIO  OF  SHEAR  STRAIN  TO  VERTICAL 
STRAIN  AT  o3c  =  0.5  KSC  FOR  SEVERAL  COMBINED  LOADING  TESTS 


FIGURE  6-5  DEGRADATION  OF  NORMALISED  COMPRESSION  MODULUS  WITH  RATIO  OF  SHEAR  STRAIN  TO 
VERTICAL  STRAIN  AT  n-  =  2.0  KSC  FOR  SEVERAL  COMBINED  LOADING  TESTS 


10’ 2  4 


FIGURE  6-6  DEGRADATION  OF  NORMALIZED  COMPRESSION  MODULUS  WITH  RATIO  OF  SHEAR  STRAIN  TO  VERTICAL 
STRAIN  AT  it-  =3.5  KSC  FOR  SEVERAL  COMBINED  LOADING  TESTS 
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FIGURE  6-8  EFFECT  OF  POISSON’S  RATIO,  u,  ON  THE  DEGRADATION  OF  THE  NORMALIZED  COMPRESSION 
MODULUS  WITH  THE  RATIO  OF  SHEAR  STRAIN  TO  VERTICAL  STRAIN  FOR  SEVERAL  TESTS 
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Octahedral  Shearing  Strain  Effects 

As  with  the  presentation  in  Chapter  5,  the  usefulness  of  the  octahe¬ 
dral  shearing  strain  in  predicting  compression-shear  interaction  was 
evaluated. 

In  Figure  6-9,  the  "virgin"  normalized  compression  moduli  curves  are 
shown  plotted  against  the  octahedral  shearing  strain.  For  the  data  in 
these  curves  the  principal  strain  directions  are  coincident  with  the 
vertical,  radial,  and  tangential  directions,  and  do  not  rotate  during 
loading. 

In  the  combined  loading  tests,  the  principal  strain  axes  in  the  6z 
"plane"  rotate  continuously  during  loading,  and  the  octahedral  shearing 
strain  arises  from  the  combination  of  normal  and  shear  strains  in  the 
vertical,  radial,  and  tangential  directions.  The  variation  of  the  nor¬ 
malized  compression  modulus  with  octahedral  shearing  strain  under  com¬ 
bined  loading  conditions  for  a  vertical  normal  strain  amplitude,  ,  of 
0.1%,  are  presented  in  Figures  6-10  through  6-12  for  lateral  confining 
pressures  of  0.5,  2.0,  and  3.5  KSC,  respectively.  These  data  are  com¬ 
bined  for  comparison  purposes  in  Figure  6-13. 

It  is  apparent  from  the  data  presented  in  Figure  6-12  that  the  mod¬ 
ulus  degradation  curve  for  a  vertical  normal  strain  of  0.1%  is  influ¬ 
enced  primarily  by  the  octahedral  shearing  strain  amplitude,  and  is  es¬ 
sentially  independent  of  the  confining  pressure.  It  is  therefore  pos¬ 
sible  to  draw  one  best-fit  curve  through  the  data  as  shown  in  Figure  6-13. 

Similarly,  a  best-fit  curve  was  drawn  through  the  data  for  a  vertical 
normal  strain  of  approximately  2X10  2%  as  shown  in  Figure  6-14. 


'  A 


OCTAHEDRAL  SHEARING  STRAIN  -  y0CT  (*') 

FIGURE  6-9  DEGRADATION  OF  NORMALIZED  COMPRESSION  MODULUS  WITH 
OCTAHEDRAL  SHEARING  STRAIN  FOR  CONDITION  OF  VERTICAL 
LOADING  ALONE 
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OCTAHEDRAL  SHEARING  STRAIN  -  Y()CT  («) 

FIGURE  6-11  DEGRADATION  OF  NORMALIZED  COMPRESSION  MODULUS  WITH 


OCTAHEDRAL  SHEARING  STRAIN  AT  03c  =  2.0  KSC  FOR 
SEVERAL  COMBINED  LOADING  TESTS  WITH  A  CONSTANT 
VERTICAL  STRAIN  OF  10'Io< 
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OCTAHEDRAL  SHEARING  STRAIN  -  yoct  (") 

FIGURE  6-14  SUMMARY  OF  DEGRADATION  OF  NORMALIZED  COMPRESSION 
MODULUS  WITH  OCTAHEDRAL  SHEARING  STRAIN  FOR  ALL 
COMBINED  LOADING  TESTS  AT  A  CONSTANT  VERTICAL 
STRAIN  OF  2 XI O- 2 % 
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These  relative  moduli-octahedral  shearing  strain  curves  were  then 
combined  with  the  virgin  curves  presented  in  Figure  6-9  to  obtain  the 
curves  presented  in  Figures  6-15  through  6-17  using  the  procedure  de¬ 
scribed  for  Figures  5-16  through  5-18.  These  plots  were  then  combined 
for  comparison  purposes  in  Figure  6-18. 

Unlike  the  results  of  the  low  strain  resonant  column  presentation  in 
Figures  5-16,  5-17,  and  5-18,  the  curves  summarized  in  Figure  6-18  show  a 
distinct  variation  depending  on  how  y^™  was  developed;  e.i.,  from  verti- 

Ut,  1 

cal  compression  or  from  shear  straining.  Of  course  part  of  the  variation 
shown  in  Figure  6-18  is  due  to  variations  in  confining  pressure,  which 
was  illustrated  in  Figures  6-15  through  6-17.  For  low  confining  pres¬ 
sures,  the  modulus  appears  to  degrade  more  rapidly  as  a  result  of  dynamic 
normal  vertical  straining  than  shear  straining.  For  higher  confining 
pressures,  the  reverse  appears  true.  Superimposed  upon  Figure  6-18  is  an 
interpolated  curve  where  it  appears  that  the  degradation  of  modulus  is 

independent  of  the  source  of  Y,^—.  This  curve  has  been  labeled  the  "no 

OCT 

interaction  line,"  and  appears  at  a  lateral  confining  pressure  of  appro¬ 
ximately  1.0  KSC. 

To  further  summarize  the  results  of  this  testing  series,  a  set  of 
best-fit  curves  have  been  presented  in  Figure  6-19  to  show  the  degrada¬ 
tion  of  normalized  compression  modulus  with  YQCT  irrespective  of  the 
strain  path.  It  is  noted  that  the  data  varies  from  these  curves  to  a 
greater  degree  than  the  comoarable  curves  in  Figure  5-19. 

High  Shear  Strain  Effects 

Several  combined  loading  tests  were  performed  under  conditions  of 


OCTAHEDRAL  SHEARING  STRAIN  -  Y0CT(") 

FIGURE  6-15  DEGRADATION  OF  NORMALIZED  COMPRESSION  MODULUS  WITH 
OCTAHEDRAL  SHEARING  STRAIN  AT  03c  =0.5  KSC  FOR 
BOTH  VERTICAL  AND  COMBINED  LOADING  TESTS 


OCTAHEDRAL  SHEARING  STRAIN  -  >0CT  (•',) 

FIGURE  6-17  DEGRADATION  OF  NORMALIZED  COMPRESSION  MODULUS  WITH 
OCTAHEDRAL  SHEARING  STRAIN  AT  a3c  =3.5  KSC  FOR 
BOTH  VERTICAL  AND  COMBINED  LOADING  TESTS 
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OCTAHEDRAL  SHEARING  STRAIN  -  y  {'■) 

FIGURE  6-18  DEGRADATION  OF  NORMALIZED  COMPRESSION  MODULUS  WITH 
OCTAHEDRAL  SHEARING  STRAIN  FOR  ALL  VERTICAL  AND 
COMBINED  LOADING  TESTS 


OCTAHEDRAL  SHEARING  STRAIN  -  vQCT  (c  ) 

FIGURE  6-19  SUMMARY  OF  DEGRADATION  OF  NORMALIZED  COMPRESSION 
MODULUS  WITH  OCTAHEDRAL  SHEARING  STRAIN 


relatively  very  high  shear  strains,  but  it  was  very  difficult,  if  not 
impossible,  to  determine  accurate  modulus  values  from  these  test  results 
because  of  the  shape  of  the  hysteretical  stress-strain  plots.  Several 
typical  plots  are  shown  in  Figure  6-20. 

The  odd  shapes  of  these  hysteretic  stress-strain  plots  provide  a  clue 
to  the  behavior  of  soils  under  tnese  loading  conditions.  It  is  noted,  for 
example,  that  several  plots  closely  approximate  a  "Figure-8"  shape,  having 
two  distinct  loops.  A  Lissajous-figure  analysis  of  a  standing  "Figure-8" 
would  indicate  an  out-of-phase  relationship,  with  the  vertical  axis  at 
half  the  frequency  of  the  horizontal  axis.  In  the  case  of  the  hollow 
cylinder  specimens  it  indicates  that  the  soil  sample  is  displaying  two 
cycles  of  vertical  strain  for  each  cycle  of  vertical  stress. 

If  we  consider  the  case  of  a  uniformly  graded  granular  material  with 
a  relatively  dense  packing,  the  matrix  of  grains  might  look  something  like 
that  shown  in  Figure  6-21 (a) .  As  large  shear  stresses  are  applied  to  the 
material  in  one  direction,  one  grain  would  tend  to  "ride-over"  another, 
resulting  in  a  vertical  displacement.  As  the  shear  stress  reversed  its 
direction  in  time,  the  grain  would  return  essentially  to  its  original 
packing  position,  then  "ride-up"  again  upon  another  grain  as  shown  in 
Figure  6-21 (b).  This  effect  would  evidence  as  two  cycles  of  vertical 
straining  for  each  shear  stress  cycle,  as  shown  in  Figure  6-21 (c). 

In  the  case  of  relatively  very  high  shear  stress  loading,  the  verti¬ 
cal  straining  resulting  from  this  process  might  be  large  in  comparison 
with  the  vertical  strain  induced  by  vertical  loading.  If  this  is  the 
case,  then  a  "Figure-8"  stress-strain  plot  will  result,  as  was  observed 


in  several  tests.  The  other  observed  odd-shaped  plots  may  be  explained 


6-20  TYPICAL  STRESS  VS.  STRAIN  CURVES  IN  VERTICAL  COMPRESSION 

UNDER  CONDITION  OP  SIMULTANEOUS  VERY  HIGH  RELATIVE  SHEAR  STRAIN 
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by  similar  reasoning,  considering  phase  effects  and  various  combinations 
of  vertical  strain  induced  by  vertical  loading  and  shear  loading. 

It  should  be  noted  at  this  point  that  the  consideration  of  a  shear 
stress  producing  a  normal  strain  is  a  coupling  effect  and  is  thus  contrary 
to  the  assumptions  of  the  applicability  of  elastic  theory  and  Hooke’s  Law 
to  this  study. 

This  postulated  theory  of  behavior  is  introduced  only  to  aid  in  the 
understanding  of  soil  behavior  under  conditions  of  very  high  shear  straining, 
where  inelastic,  non-linear  effects  are  apparent. 

Conclusions 

As  in  Chapter  5,  the  test  results  from  this  testing  series  were  pre¬ 
sented  in  different  ways.  The  degradation  of  modulus  was  presented  both 
as  a  function  of  the  strain  ratio,  yQ  /£  ,  and  as  a  function  of  the  octa- 

DZ  Z 

hedral  shearing  strain,  y  .  High  shear  strain  effects  were  also  ex- 

OC  a 

plored. 

Strain  Ratio  Effects 

The  result  of  the  first  presentation,  summarized  in  Figure  6-6,  in¬ 
dicates  that  there  exists  a  relatively  consistent  threshold  strain  ratio 
below  which  interaction  effects  are  negligible.  In  the  case  of  the  high 

vertical  normal  strain  tests,  where  e  =0.1%,  the  threshold  strain  ratio 

z 

fell  to  as  low  as  0.5,  indicating  that  shear  strains  below  0.05%  would 
have  negligible  effect  in  further  degradation  of  modulus. 

As  with  the  low  strain  resonant  column  testing  series,  the  normalized 
modulus  curves  were  found  to  be  relatively  independent  of  relative  density 
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within  the  range  tested  but  were  moderately  dependent  upon  confining  pres¬ 
sure.  The  magnitude  of  the  interaction  effects  appears  to  be  greater  at 
lower  confining  pressures.  This  effect  may  be  explained  by  the  fact  that 

the  principal  stress  ratios,  (0,/cr)  during  cyclic  loading  increase 

1  J  MAX 

very  rapidly  at  lower  confining  pressures  as  shown  in  Table  6-2;  so  that 
the  soil  material  more  closely  approaches  a  failure  condition  in  these 
tests.  If  the  incremental  modulus  degrades  more  rapidly  with  higher 
stress  increments,  which  was  observed  in  the  hysteretic  stress-strain 
plots  for  these  tests,  then  greater  straining  (and  greater  interaction) 
is  to  be  expected. 

Octahedral  Shearing  Strain  Effects 

The  second  presentation  was  of  the  degradation  of  modulus  with  the 

octahedral  shearing  strain,  y  ,  summarized  in  Figures  6-18  and  6-19. 

OCT 

Unlike  the  low  strain  resonant  column  testing  series  results,  there  was 
some  variation  in  the  degradation  curves  depending  upon  how  the  y  was 
developed. 

This  variation  is  proportional  to  the  interaction  effect  not  ac¬ 
counted  for  by  using  yocT*  However,  it  was  possible  to  select  best-fit 
curves  of  normalized  modulus  vs.  y„__  (Figure  6-19)  for  which  the  varia- 
tions  from  the  best-fit  curves  is  quite  acceptably  small,  except  where 
Yqct  becomes  very  large.  When  Y0CT  becomes  very  large  and  the  confining 
pressure  is  relatively  low,  the  assumptions  on  which  the  various  deriva¬ 
tions  were  made  become  invalid  as  discussed  earlier. 

Thus  the  accuracy  of  the  normalized  curves  presented  in  Figure  6-19 
decreases  progressively  with  increasing  y 0CT»  as  can  be  judged  by  com- 


paring  them  with  the  curves  of  Figure  6-18  from  which  they  were  derived. 

The  highest  octahedral  shearing  strain  used  in  this  testing  series 
was  approximately  0.5%. 

High  Shear  strain  Effects 

An  analysis  of  soil  behavior  involving  the  modeling  of  uniformly 
graded  granular  soil  as  a  collection  of  discrete  particles  was  presented 
to  help  explain  the  observed  soil  response  to  high  cyclic  shear  strains 
imposed  simultaneously  with  relatively  low  amplitude  vertical  loading. 
The  test  data  presented  in  Figure  6-20  indicated  inelastic,  non-linear 
response  of  the  soil  material;  specifically,  significant  vertical  normal 
strains  were  observed  to  result  from  the  application  of  high  horizontal 
shear  stresses.  When  such  coupling  response  is  significant,  elastic 
theory  and  Hooke's  Law  would  not  apply. 
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Chapter  7 

Comparison  of  Laboratory  Test  Results 

Introduction 

In  both  the  triaxial  resonant  column  testing  series  and  the  thin- 
walled  hollow  cylinder  testing  series,  the  effects  of  combined  compression 
and  shear  loading  were  explored.  The  two  testing  series  were  conducted 
differently,  producing  different  stress  histories  and  different  ranges  of 
strain  amplitudes.  It  is  the  purpose  of  this  chapter  to  further  develop 
the  data  presented  from  these  two  testing  series  in  Chapters  5  and  6,  and 
to  develop  a  set  of  curves  showing  the  loading  effects  over  the  full  range 
of  strains  observed  during  testing. 

Vertical  Loading  Alone 

The  compression  modulus  was  plotted  against  the  vertical  strain  under 

conditions  of  vertical  compression  loading  alone  in  Figures  5-2  and  6-2. 

These  figures  were  then  normalized  to  the  value  of  compression  modulus  at 

a  vertical  strain  of  £  =  5  X  10  3%,  with  the  resulting  curves  shown  in 

z 

Figures  5-3  and  6-3.  The  value  of  selected  for  normalization  was 
chosen  to  maximumize  the  amount  of  overlap  in  the  data  and  to  provide  a 
simple  means  of  combining  the  two  sets  of  results. 

The  data  from  Figures  5-3  and  6-3  were  combined  and  presented  in 
Figure  7-1.  These  curves  were  then  replotted  in  Figure  7-2  normalized  to 
the  maximum  compression  modulus  at  a  vertical  strain  of  =  10  5%. 
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Mean  Confining  Stress 


As  indicated  in  these  figures,  the  combined  results  were  plotted  for 
different  values  of  lateral  confining  pressure,  0  .  During  the  triaxial 
resonant  column  testing  series,  samples  were  consolidated  in  such  a  way 


that  normal  stresses  in  every  direction  were  equal  to  a„  ;  so  that  the 

3c 

mean  confining  stress  during  consolidation ,  a  ,  may  be  written: 


a 

m 


(7.1) 


With  the  thin-walled  hollow  cylinder  test  specimens,  however,  the 
principal  stresses  are  not  all  equal  during  consolidation,  but  are  as 
follows: 


o  =  1.85  •  0 

lc  3c 


=  U  • 


,0le  *  °3c>' 


And  where  u  =  0.34,  cr  may  be  written: 

tn 

O  =  1.273  •  0,  (7.6) 

m  3c 

Because  the  degradation  of  modulus  with  vertical  strain  for  the  two 
testing  series  may  be  better  represented  as  a  function  of  the  mean  con¬ 
fining  stress  than  the  lateral  confining  stress  alone,  it  is  desirable  to 
replot  the  data  from  Figure  6-3  for  values  of  0^  of  0.5,  2.0,  and  3.5  KSC. 
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In  order  to  accomplish  this  replot,  the  variation  of  the  normalized  com¬ 
pression  modulus  was  first  plotted  auainst  the  mean  confining  pressure, 

3  ,  for  a  constant  range  of  strain.  This  plot  is  shown  in  Figure  7-3,  with 

the  vertical  strain  range  of  from  e  =  5X10  H  to  £  =  10  '%.  The  circled 

z  z 

data  points  in  this  figure  come  directly  from  the  data  in  Figure  6-3, 
where  is  related  to  a^c  by  Equation  7.6.  The  triangles  represent  in¬ 
terpolated  data  points  for  0  values  of  0.5,  2.0,  and  3.5  KSC. 

m 

With  the  data  represented  by  the  triangles  in  Figure  7-3,  it  was  pos¬ 
sible  to  replot  Figure  6-3  to  reflect  the  desired  values  of  0  .  This  plot 

m 

is  shown  in  Figure  7-4. 

From  Figures  7-4  and  5-2  it  was  possible  to  replot  the  normalized 

compression  modulus-vertical  strain  curves  in  Figure  7-2  in  terms  of  the 

mean  confining  stress,  O  .  This  new  plot  is  shown  in  Figure  7-5.  It  is 

m 

noteworthy  that  the  difference  between  this  figure  and  Figure  7-2  is  small. 

There  is  a  somewhat  more  rapid  degradation  of  modulus  with  increasing 

strain  for  the  a  curves  than  for  the  comparable  o„  curves. 

m  3c 

Octahedral  Shearing  Strain 

In  Chapter  5  it  was  noted  that  the  normalized  compression  modulus 
degraded  with  the  octohedral  shearing  strain  relatively  independently  of 
how  that  strain  was  developed:  e.i.,  whether  from  vertical  loading  alone 
or  from  combined  vertical  and  torsional  loading.  It  was  concluded  there¬ 
fore,  that  for  the  resonant  column  data,  compression-shear  interaction 
could  be  completely  accounted  for  by  using  the  octahedral  shearing  strains. 

In  Chapter  6  some  compression-shear  interaction  effects  (for  hollow 
cylinder  specimens)  were  noted  which  were  not  completely  accounted  for  by 
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MEAN  CONFINING  PRESSURES 


NORMALIZED  COMPRESSION  MODULUS 


VERTICAL  STRAIN  cz  (?) 

FIGURE  7-4  DEGRADATION  OF  NORMALIZED  COMPRESSION  MODULUS  WITH  VERTICAL 
STRAIN  FOR  HOLLOW  CYLINDER  TESTS  UNDER  CONDITION  OF  VERTICAL 
LOADING  ALONE  FOR  MEAN  CONFINING  PRESSURE,  c  ,  OF  0.5,  2.0, 
AND  3.5  KSC  m 
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FIGURE  7-5  DEGRADATION  OF  NORMALIZED  COMPRESSION  MODULUS  WITH  VERTICAL  STRAIN  UNDER  CONDITION  OF 
VERTICAL  LOADING  ALONE  AT  MEAN  CONFINING  PRESSURES,  a  ,  OF  0.5,  2.0,  AND  3.5  KSC 
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the  use  of  the  octahedral  strain.  Nevertheless,  it  was  possible  to 
develop  a  family  of  "best  fit"  octahedral  strain  curves  which  represent 
an  excellent  approximation  of  the  variation  of  normalized  compression 
modulus  with  octahedral  shearing  strain  under  a  wide  range  of  combined 
loading  conditions. 

The  octahedral  shearing  strain  curves  were  presented  in  Figures  5-19 

and  6-19  for  the  two  testing  series.  Using  the  technique  described  above 

to  develop  a  series  of  curves  for  values  of  a  of  0.5,  2.0,  and  3.5  KSC, 

m 

these  two  figures  were  combined  as  shown  in  Figure  7-6.  This  figure  may 
be  used  with  reasonable  accuracy  to  predict  the  degradation  of  modulus 
with  octahedral  shearing  strain  for  Monterey  No.  0  sand  or  a  similar  sand. 

Strain  Ratio  Effects 

Another  way  of  evaluating  the  compression-shear  interaction  effects 
upon  the  degradation  of  modulus  with  strain  was  discussed  in  Chapters  5 
and  6,  and  was  referred  to  as  the  strain  ratio  effect.  Simply  stated, 
the  method  involved  the  use  of  the  strain  ratio  degradation  curves  pre¬ 
sented  in  Figures  5-10  and  6-7,  where  the  degree  of  degradation  was  shown 
as  a  function  of  the  ratio  of  shear  strain  to  vertical  normal  strain  for 
specific  values  of  vertical  normal  strain.  To  predict  the  degradation 
of  compression  modulus  under  any  straining  condition,  one  would  first 
determine  the  degradation  resulting  from  the  vertical  straining  alone; 
then  that  degradation  would  be  further  reduced  by  the  reduction  factor 
determined  from  the  strain  ratio  figures. 

The  initial  degradation  of  modulus  resulting  from  vertical  straining 
alone  may  now  be  obtained  directly  from  Figure  7-5.  To  provide  a  summary 


1 


0.34 


figure  for  obtaining  the  additional  degradation  resulting  from  the  exis¬ 
tence  of  simultaneous  shear  straining,  a  set  of  strain  ratio  curves  were 
plotted  for  values  of  vertical  strain  of  10  1 ,  10  2 ,  10  3,  and  10  4%. 

These  curves  appear  in  Figure  7-7. 

Conclusions 

The  presentations  made  in  Chapters  5  and  6  regarding  the  degradation 
of  modulus  with  strain  under  combined  compression-shear  loading  conditions 
were  further  developed  to  provide  data  curves  over  the  entire  strain  range 
observed  during  the  two  testing  series.  These  curves  may  be  used  to  pre¬ 
dict  degradation  of  modulus  under  conditions  of  simultaneous  compression- 
shear  loading  with  reasonable  accuracy. 

Octahedral  Shearing  Strain  Effects 

The  degree  of  degradation  of  modulus  may  be  predicted  by  calculating 
the  octahedral  shear  strain  from  the  complete  strain  state,  and  interpo¬ 
lating  from  the  curves  in  Figure  7-6.  This  degradation  may  then  be  ap¬ 
plied  to  an  estimate  of  the  maximum,  low-strain  modulus  to  determine  the 
modulus  corresponding  to  a  given  state  of  strain.  The  maximum  modulus 
value  may  be  obtained  either  from  a  laboratory  or  field  test,  such  as  a 
seismic  survey  to  obtain  P-wave  velocity,  or  from  published  data  showing 
typical  values  of  low-strain  compression  moduli  at  various  mean  confining 
stresses.  The  values  of  low-strain  moduli  observed  in  these  testing  series 
are  shown  in  Figures  5-2  and  6-2,  and  are  available  in  more  detail  else¬ 


where  (Griffin,  1980) . 
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Strain  Ratio  Effects 

The  degree  of  degradation  of  modulus  may  also  be  estimated  if  the 
vertical  normal  strain  and  the  horizontal-vertical  shear  strain  are  known. 
The  degradation  resulting  from  vertical  straining  alone  is  first  deter¬ 
mined  from  Figure  7-5,  then  the  additional  degradation  resulting  from 
shear  straining  is  obtained  by  interpolating  from  the  curves  in  Figure 
7-7.  Because  of  the  difficulty  in  making  the  logarithmic  interpolations 
in  Figure  7-7  the  octahedral  shearing  strain  technique  described  earlier 
appears  more  promising  for  estimating  compression  modulus  values  for 
analysis  purposes. 
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Chapter  8 

Large  Scale  Model  Studies 

Introduction 

It  is  the  purpose  of  this  chapter  to  present  the  observed  soil  re¬ 
sponse  for  a  series  of  large  scale  dynamic  model  studies  on  slopes  of 
sand.  Because  the  significant  response  is  primarily  due  to  yielding  within 
the  soil  slope,  which  is  beyond  the  elastic  range  of  loading  for  the  soil 
material,  conventional  elastic  response  analysis  techniques  are  of  little 
value.  An  attempt  will  be  made  to  predict  the  point  at  which  yielding 
begins  using  an  analysis  technique  proposed  by  Seed  and  Goodman  (1964) . 

In  earlier  chapters,  compression-shear  interaction  effects  were 
studied  under  a  variety  of  combined  loading  conditions  on  laboratory  sam¬ 
ples  of  sand.  Interaction  effects  were  presented  in  such  a  manner  as  to 
show  their  influence  upon  the  dynamic  moduli  of  the  sands  under  various 
conditions  of  density,  confining  pressure,  and  strain  amplitude.  One 
major  advantage  of  this  form  of  presentation  was  that  there  exist  a  large 
number  of  elastic  analysis  techniques  which  employ  the  dynamic  moduli  to 
predict  soil  response,  thus  allowing  for  relatively  simple  application  of 
these  research  findings  in  conventional  elastic  analysis  methods. 

In  order  to  evaluate  any  interaction  effects  occurring  beyond  the 
elastic  loading  range  of  the  soils,  it  was  necessary  to  perform  a  soil 
loading  test  in  which  yielding  occurs  with  simultaneous  shear-compression 
dynamic  excitation,  and  for  which  an  accurate  analysis  technique  has  been 
developed.  These  slope  model  studies  appeared  to  meet  those  criteria. 
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Test  Set-up 

The  model  studies  were  conducted  on  the  20  ft  x  20  ft  shaking  table 
located  at  the  Earthquake  Simulator  Laboratory  at  the  Richmond  Field  Sta¬ 
tion,  Earthquake  Engineering  Research  Center,  University  of  California, 
Berkeley,  California.  This  table  can  move  in  one  horizontal  direction 
and  the  vertical  direction.  It  was  designed  to  exactly  reproduce  dynamic 
accelerations,  velocities,  and  displacements,  within  certain  limits,  from 
analog  records  stored  on  either  magnetic  tape  or  disk.  A  detailed  descrip¬ 
tion  of  this  shaking  table  has  been  provided  by  Rea  (1972)  and  by  Rea  and 
Penzien  (1972) . 

The  cross-section  of  a  typical  soil  slope  specimen  is  shown  in 
Figure  8-1.  In  this  figure,  four  accelerometers  and  four  DC  linear  vari¬ 
able  displacement  transformers  (labeled  DCDT)  are  shown.  In  addition  to 
these,  the  table  itself  is  heavily  instrumented  with  accelerometers  and 
displacement  transducers . 


Test  Box 

The  inside  dimensions  of  the  test  box  in  which  the  soil  slope  speci¬ 
mens  are  formed  are  approximately  84.2  in  x  42.5  in  x  21.7  in,  with  a  tri¬ 
angular  spacer  in  one  corner  approximately  15.8  in  on  each  side,  placed 
at  an  angle  of  approximately  45°  with  the  horizontal.  The  test  box  was 
bolted  securely  to  the  shaking  table,  and  a  layer  of  high-strength 
"hydrostone"  cement  mortar  was  placed  between  the  box  and  the  table,  and 
in  the  base  of  the  box  to  assure  rigid,  solid  contact  between  the  two. 

The  test  box  is  shown  in  Figure  8-2 (a) ,  and  the  bracing  behind  the 


tIL  OF  TEST  B' 
NO  SPACER  BO, 
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spacer  board  is  shown  in  Figure  8-2  (b) .  The  bracing  was  installed  to  in¬ 
sure  a  stiff  response  of  the  spacer  board  during  loading,  so  that  its  de¬ 
flection  during  loading  would  be  very  small  when  compared  with  the  dis¬ 
placements  of  the  soil  slope. 

Before  the  test  box  was  used  for  this  testing  series,  a  thin  film  of 
high-strength  epoxy  was  sprayed  upon  the  back,  spacer  board,  and  bottom, 
and  a  thin  layer  of  the  test  sand  placed  upon  those  surfaces. 

The  sand  used  in  this  testing  series  was  the  same  Monterey  No.  0  sand 
used  in  the  hollow  cylinder  and  resonant  column  testing  series. 

Forming  Test  Specimens 

Slope  specimens  were  formed  by  placing  sand  into  a  "mold"  formed 
within  the  test  box,  and  vibrating  the  sand,  one  lift  at  a  time,  until  the 
desired  slope  height  was  obtained.  Vibration  was  done  by  the  shaking 
table.  The  forming  "mold"  assembly  is  shown  in  Figures  8-3(a)  and  (b) . 

A  typical  completed  soil  specimen  is  shown  in  Figures  8-4  (a)  and  (b) , 
with  instrumentation  in  place.  A  closeup  of  the  placement  of  the  DCDT's 
on  the  soil  slope  surface  is  shown  in  Figure  8-5  (a).  As  placed,  the  core 
of  the  DCDT's  have  a  very  thin  (.01  in  diameter),  stiff  2  in  length  of 
piano  wire  protruding  from  them.  This  piano  wire  is  coated  with  a  very 
thin  film  of  epoxy,  and  carefully  inserted  into  the  soil  slope  surface. 

The  DCDT's  are  then  zeroed  by  adjusting  their  position  relative  to  their 
cores  in  the  mounting  bracket.  Once  the  epoxy  film  hardens,  the  sand  im¬ 
mediately  surrounding  the  piano  wire  protrusion  is  cemented  to  the  wire, 
but  its  fabric  remains  relatively  undisturbed.  A  closeup  of  this  wire, 
epoxy,  sand  cementation  is  shown  in  Figure  8-5 (b). 


■  V4  -'-V1'' 


FIGURE  8-5  DETAIL  OF  PLACEMENT  OF  DCDT  CORE  PROBES 
ON  THE  SURFACE  OF  A  SLOPE  SPECIMEN  (a), 
AND  OF  SAND-EPOXY  BOND  ON  WIRE  PERTUSION 
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Summary  of  Tests 

The  tests  performed  in  this  testing  series  are  summarized  in  Table 

8-1. 

Soil  Strength  Characteristics 

Because  the  yielding  which  occurs  along  the  outer  surface  of  the  soil 
slopes  is  of  primary  importance  to  this  study,  it  is  important  to  know  ac¬ 
curately  the  shear  strength  characteristics  of  the  outermost  soil  slope 
materials . 

The  soil  slopes  were  formed  by  vibrating  soils  within  a  forming  mold, 
which  was  later  removed.  It  is  highly  probable  that  passive  pressures  and 
other  stress  concentrations  developed  during  sample  formation  resulted  in 
non-uniform  distribution  of  densities  within  the  slopes.  It  is  therefore 
necessary  to  determine  the  shear  strength  characteristics  of  the  outer¬ 
most  soil  materials  by  methods  other  than  conventional  shear-strength  vs 
density  correlations. 

Static  Slope  Failure 

In  Test  No.  9,  shown  in  Table  8-1,  a  marginal  static  slope  failure 
occurred  soon  after  the  slope  was  formed.  This  observed  failure,  which 
was  a  shallow  sloughing  failure,  provides  a  probable  upper  bound  for  the 
angle  of  internal  friction,  <j>,  of: 

afailure 


MAX 


37.6 


(8-1) 


TABLE  8-1 


Summary  of  Soil  Slope  Model  Tests 
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where  a  is  the  angle  of  the  slope, 

D  is  the  relative  density  of  the  soil, 

R 

av  is  the  peak  vertical  acceleration, 

a,  is  the  peak  horizontal  acceleration,  and 
h 

oo  is  the  angle  of  application  of  the  resultant  acceleration 


Direct  Shear  Tests 


During  the  testing  series,  two  groups  of  special  direct-shear  tests 
were  performed  on  the  outermost  soils  of  the  slopes  after  they  had  been 
formed  and  were  ready  for  testing.  These  direct-shear  tests  were  per¬ 
formed  by  placing  a  lightweight  block ,  specially  prepared  with  a  thin 
layer  of  sand  epoxied  to  it,  upon  the  surface  of  the  slope.  The  block 
was  then  pulled  upslope  with  a  string  until  yielding  occurred.  As  shown 
in  Figure  8-6,  if  the  weight  of  the  block  and  the  string  tension  are 
known,  the  normal  and  shear  stresses  on  the  failure  surface  at  yielding 
are  easily  calculated. 

In  Tables  8-2  and  8-3  the  calculated  values  of  N  and  S  and  the  normal 
and  shear  stresses  on  the  failure  surface  at  yielding  are  shown  for  aver¬ 
age  relative  densities  of  95%  and  80%  respectively.  From  these  tables, 
the  shear  strength  plots  shown  in  Figures  8-7  and  8-8  have  been  developed. 
From  these  figures  it  appears  that  the  shear  strength  characteristics  of 
the  outermost  slope  materials  are  approximately  the  same  even  when  the 
average  soil  densities  are  not,  and  are  defined  by  the  following  equation: 

S  =  s.  •  £  +  N  •  tan(}>  (8.2) 

l 

where  s^  is  the  shear  strength  intercept,  1  is  the  length,  and  <J>  is  the 
angle  of  internal  friction.  Equation  8.2  may  be  rewritten  as  follows: 

S  =  0.029  •  i  +  N  •  tan  (35.4°)  (8.3) 

It  is  interesting  to  note  at  this  point  that  the  value  of  s^  is  very 
small,  and  that  the  values  of  s^  and  <J>  are  much  smaller  than  those  used 
by  Seed  and  Goodman  in  their  studies  involvinq  Monterey  sand.  There  are 
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(a) 


FIGURE  8-6  DIAGRAM  OF  FORCES  ACTING  ON  THE  SLIDING 
BLOCK  USED  IN  SURFACE  DIRECT  SHEAR  TESTS 
PERFORMED  ON  SLOPE  SPECIMENS 


152 


TABLE  8-2 


N  and  S  for  DR  -  95% 


Test 

ot(°) 

W  ( lb . ) 

W  cosa 

A 

N  (psi) 

W  sina 

A  _ 

-Sw(psi)  +  T(psi) 

=  S (psi) 

060576.0 

33.5 

0.1552 

.0032 

.0021 

.0040 

.0019 

0.5847 

.012 

.008 

.0164 

.0084 

1.2339 

.0254 

.0168 

.0333 

.0165 

0.3849 

.0079 

.0052 

.0098 

.0046 

0.1552 

.0032 

.0021 

.0047 

.0026 

0.5847 

.012 

.008 

.0166 

.0086 

0.1552 

.0032 

.0021 

.0041 

.002 

0.5847 

.012 

.008 

.0151 

.0071 

100576.1 

35.3 

0.1552 

.0031 

.0022 

.0047 

.0025 

0 . 3849 

.0078 

.0055 

.0121 

.0066 

0.5847 

.0118 

.0083 

.018 

.0097 

0.8144 

.0164 

.0116 

.0249 

.0133 

0.1552 

.0031 

.0022 

.0046 

.0024 

0.3849 

.0078 

.0055 

.0112 

.0057 

0.5847 

.0118 

.0083 

.0179 

.0096 

0.8144 

.0164 

.0116 

.0257 

.0141 

0.1552 

.0031 

.0022 

.0044 

.0022 

0.3849 

.0078 

.0055 

.0105 

.005 

0.5847 

.0118 

.0083 

.0163 

.008 

0.8144 

.0164 

.0116 

.0238 

.0122 

0.1552 

.0031 

.0022 

.0043 

.0021 

0.3849 

.0078 

.0055 

.0106 

.0051 

0.5847 

.0118 

.0083 

.0156 

.0073 

0.8144 

.0164 

.0116 

.0216 

.01 

0.3849 

.0078 

.0055 

.011 

.0055 

0.8144 

.0164 

.0116 

.0238 

.0122 

110576.1 

35.7 

0.1552 

.0031 

.0022 

.0047 

.0025 

1.2339 

.0247 

.0178 

.0339 

.0161 

0.1552 

.0031 

.0022 

.0043 

.0021 

0.5847 

.0117 

.0084 

.017 

.0086 

1.2339 

.0247 

.0178 

.0339 

.0161 

2.7893 

.0559 

.0402 

.0741 

.0339 

l  -  -  ■  .  ,  -  - _ 

1.2339 

.0247 

.0178 

.0336 

.0158 

a  =  angle  of  slope 


A  =  cross-sectional  area  of  block 
T  =  tension  in  string 
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TABLE  8-3 

N  and  S  for  DR  -  80% 


Test 

a(°) 

W (lb . ) 

W  cosa 

A 

w  sina 

A 

+  Tvpsi) 

=  S (psi ) 

N  (psi) 

-S  (psi) 
w 

140576.1 

35.2 

.1552 

.0031 

.0022 

.0044 

.0022 

.5847 

.0118 

.0083 

.0164 

.0081 

.8045 

.0162 

.0115 

.0226 

.0111 

1.234 

.0249 

.0176 

.0345 

.0169 

2.789 

.0563 

.0397 

.0774 

.0377 

.5847 

.0118 

.0083 

.0162 

.0079 

1.234 

.0249 

.0176 

.0351 

.0175 

2.789 

.0563 

.0397 

.0765 

.0368 

.1522 

.0031 

.0022 

.0044 

.0022 

.5847 

.0118 

.0083 

.0173 

.0090 

.8045 

.0162 

.0115 

.0234 

.0119 

1.234 

.0249 

.0176 

.0360 

.0184 

2.789 

.0563 

.0397 

.0785 

.0388 

170576.1 

35.4 

.1552 

.0031 

.0022 

.0045 

.0023 

.5847 

.0118 

.0083 

.0163 

.0080 

1.234 

.0248 

.0177 

.0347 

.0170 

2.789 

.0561 

.0399 

.0761 

.0362 

.8045 

.0162 

.0115 

.0226 

.0111 

.5847 

.0118 

.0083 

.0166 

.0083 

1.234 

.0248 

.0177 

.0348 

.0171 

2.789 

.0561 

.0399 

.0812 

.0413 

.1552 

.0031 

.0022 

.0044 

.0022 

.5847 

.0118 

.0083 

.0170 

.0087 

.8045 

.0162 

.0115 

.0241 

.0126 

1.234 

.0248 

.0177 

.0351 

.0174 

2.789 

.0561 

.0399 

.0797 

.0398 

.5847 

.0118 

.0083 

.0170 

.0087 

1.234 

.0248 

.0177 

.0358 

.0181 

190576.1 

35.3 

.1552 

.0031 

.0022 

.0045 

.0023 

.5847 

.0118 

.0083 

.0163 

.0080 

.8045 

.0162 

.0115 

.0239 

.0124 

1.2  34 

.0249 

.0176 

.0352 

.0176 

2.789 

.0562 

.0398 

.0812 

.0414 

1.2  34 

.0249 

.0176 

.0345 

.0169 

2.789 

.0562 

.0398 

.0788 

.0390 

.1552 

.0031 

.0022 

.0045 

.0023 

.5847 

.0118 

.0083 

.0164 

.0081 

.8045 

.0162 

.0115 

.0222 

.0107 

1.234 

.0249 

.0176 

.0344 

.0168 

2.789 

.0562 

.0398 

.0752 

.0354 

3.439 

.0693 

.0491 

.0910 

.0419 
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TABLE  8-3  continued 


N  and  S  for  DR  =  80% 


Test 

ot(°) 

W  ( lb . ) 

W  cosa 

W  sina 

+  T(psi) 

=  S (psi ) 

A 

A 

N (psi ) 

-S  (psi) 
w 

250576.1 

35.4 

.1552 

.0031 

.0022 

.0045 

.0023 

.5847 

.0118 

.0083 

.0170 

.0087 

1.234 

.0248 

.0177 

.0370 

.0193 

2.789 

.0561 

.0399 

.0811 

.0412 

.1552 

.0031 

.0022 

.0049 

.0027 

.5847 

.0118 

.0083 

.0170 

.0087 

1.234 

.0248 

.0177 

.0183 

.1552 

.0031 

.0022 

.0045 

.0023 

.5847 

.0118 

.0083 

.0173 

.0090 

1.234 

.0248 

.0177 

.0168 

2.789 

.0561 

.0399 

.0810 

.0411 

155 


OF  N  AND  S  FOR  SLOPE  SPECIMEN  SANDS  HIT H  AN  AVERAGE  RELATIVE  DENSITY, 


ART AT  I ON  OF  N  ANO  S  TOR  SLOPE  SPECIMEN  SANDS  WITH  AN  AVERAGE  RCLATIVE  DENSITY 


several  reasons  for  this  difference,  one  of  which  is  the  fact  that  dif¬ 
ferent  gradations  of  Monterey  sand  were  used  in  the  Seed  and  Goodman 
study  than  were  used  in  this  testing  program.  A  gradation  analysis  of  the 
Monterey  No.  20  sand  used  by  Seed  and  Goodman  is  shown  alongside  that  for 
the  Monterey  No.  0  sand  used  in  this  study  in  Figure  8-9. 

The  low  value  of  the  shear  strength  intercept,  s^,  is  expected  both 
because  of  the  low  mean  grain  size  of  the  sand,  and  because  the  outer 
slope  material  is  probably  at  lower  density  than  the  average  slope  material 
density  because  of  arching  and  stress  concentration  effects  during  the  sam¬ 
ple  formation  process.  The  effect  of  mean  grain  size  on  the  value  of  s. 
is  shown  in  Figure  8-10,  as  developed  by  Seed  and  Goodman.  The  mean  grain 
size  of  Monterey  No.  0  sand  (.015  millimeters)  would  correspond  to  a  very 
low  value  for  . 

The  larger  degree  of  uniformity  of  the  gradation  analysis  curve  for 
Monterey  No.  0  sand  would  support  a  lower  friction  angle,  <p,  for  that  sand. 

Shaking  Table  Tests 

As  described  earlier,  slope  specimens  were  subjected  to  sinusoidal 
dynamic  loading  in  either  the  horizontal  direction  only,  or  in  combined 
horizontal  and  vertical  excitation.  Combined  horizontal  and  vertical 
loading  was  accomplished  in  phase,  as  illustrated  in  Figure  8-11.  The  re¬ 
sultant  excitation  which  the  samples  experience  is  as  shown  in  Figure 
8-11 (c),  inclined  at  an  angle  to  with  the  horizontal  plane.  The  two 
loading  conditions  are  further  illustrated  in  Figure  8-12.  The  horizon¬ 
tal  component  of  loading  was  held  constant  for  the  different  testing  se¬ 
ries,  so  that  the  resultant  acceleration  amplitude  in  the  combined  loading 


SIEVE  SIZES 


FIGURE  8-9  COMPARISON  OF  GRADATION  ANALYSES  FOR 

MONTEREY  NO.  20  AND  MONTEREY  NO.  0  SANDS 
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(a) 


(b) 


FIGURE  8-12  GRAPHICAL  ILLUSTRATION  OF  APPLICATION  OF 
HORIZONTAL  ACCELERATION  ALONE  (a)  AND 
COMBINED  HORIZONTAL  AND  VERTICAL  ACCELERATION 
(b)  ON  SLOPE  SPECIMENS 
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tests  was  necessarily  higher  than  in  the  horizontal  loading  tests. 

Test  Records 

The  two  basic  testing  records  utilized  in  this  testing  series  were 
records  A  and  B  as  shown  in  Figure  8-13.  Record  A  consists  of  four  ac¬ 
celeration  amplitude  steps  of  ten  cycles  each  increasing  in  equal  incre¬ 
ments  to  a  maximum  acceleration  amplitude  of  0.35  g.  Record  B  contains 
fifteen  acceleration  amplitude  steps  of  five  cycles  each  increasing  in 
equal  increments  to  a  maximum  acceleration  amplitude  of  0.175  g. 

Analysis 

In  accordance  with  the  analysis  technique  proposed  by  Seed  and  Good¬ 
man  (1964) ,  the  sliding  surface  at  which  limiting  equilibrium  exists  when 
the  critical  yield  acceleration  is  applied  is  approximated  as  shown  with 
the  dashed  line  in  Figure  8-12,  at  a  depth  d  below  the  surface  of  the 
slope.  The  forces  acting  upon  an  element  of  soil  along  this  sliding  sur¬ 
face  are  illustrated  in  Figure  8-14 (a).  To  satisfy  equilibrium,  Fu  must 
equal  Fd  in  this  figure,  and  the  remaining  forces  may  be  evaluated  as 
follows : 

N  =  W*cosa  -  ky>W'sin(a  +  to)  ,  and  (8.4) 

S  =  W-sino  +  ky.W.cos  (a  +  to)  (8.5) 

But,  from  Equation  8.2  we  know  that: 

S  =  sA  .  I  +  N-tan4>  (8.2) 

Equations  8.4  and  8.5  may  now  be  written: 
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S  =  s  •  b  +  Iw.  cosa  -  ky-W-sin(a  +  w)]-tan4> 
1 


=  w.sina  +  [ky-W*cos(a  +  w) ] 
Collecting  terms  and  writing: 

w  *  b  •  d  *  Y<j 

it  is  found  that: 


s . 

i 


d-Y 


+  cosa*  tarup  -  sina 


ky  = 


d 


costa  +  u>)  +  sin  (a  +  a>)‘tan4> 


(8.6) 


(8.7) 


(8.8) 


or: 


ky  =  + 


s. 

1 


cos(4>  -  a  -  oo)  d*Yj  •  [cos(a+oo)  +  sin(a+a))  •  tan4>] 
a 


(8.9) 


For  the  slope  of  finite  length  used  in  this  testing  series,  a  passive 
wedge  will  form  at  the  toe  of  the  slope  as  shown  in  Figure  8-14 (b).  The 
primary  effect  of  this  passive  wedge  is  to  resist  sliding  not  only  by  the 
shear  resistance  along  the  base  of  the  sliding  mass,  S,  but  also  by  the 
force,  P,  at  the  toe  of  the  slope. 

Summing  the  forces  perpendicular  to  the  resultant  force,  R,  gives 
the  following  expression: 


[P  •  cos(2$  -  9)  1  +[ky-W-cos  ($  -  to) )  = 

W-sin$  +  [s.  •  )L_  •  cosi)>]  +  [s.  •  l  •  sin (<j>  -  0)] 
1  Ao  X  BC 


Noting  that. 

z 

d  . 

*  *, ■■  .  ana 

AS 

sina  ' 

*BC 

d 

cos(a  +  0)  ' 

and  solving  for  P  leads 

to  the  following 

(8.10) 


(8.11) 

(8.12) 
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P  =  W-  tsin0  -  ky*cos(0  -  co)  1  +  +  sin(({)  9)1 

_ i  Lsing  cos (g  +  Q)J 


Writing 


W  = 


cos (2$  -  6) 

Yd  •  d  •  v 


(8.13) 


(8.14) 


where 


2-ac  =  d  *  [tana  +  tan(a  +  6)] 


(8.15) 


and  combining  and  rewriting  Equation  8.13,  the  following  is  obtained: 


d2  *Y 


P  = 


2 cos  (20 
2s 


-©>  [{* 


tana  +  tan(a+0)|  *|sin0  -  kycos (4>— to) 


ityT ' cos<t>  ‘ 


sin (0-9) 
sina  ‘  cos(a+9)  *cos0 


(8.16) 


As  pointed  out  by  Seed  and  Goodman,  it  is  simply  necessary  to  substi¬ 
tute  the  appropriate  values  of  a,  s  0  and  y  ,  and  a  close  approximation 

d 

for  ky,  to  obtain  a  relationship  for  P  in  terms  of  d  and  0.  By  means  of 
several  trials,  the  values  of  9  giving  the  minimum  value  of  P  for  differ¬ 
ent  values  of  d  can  then  be  determined. 

The  influence  of  the  force,  P,  can  then  be  represented  approximately 
by  an  equivalent  shear  resistance,  s  ,  as  follows: 


se  =  P/L 


(8.17) 


where  L  is  the  length  of  the  slope.  This  is  analogous  to  increasing  the 
shear  strength  intercept,  s^,  to  a  value  of  s^  +  s^,  so  that  the  cor¬ 
responding  value  of  ky  is  given  by: 


s  s 

sin(0  -  g) _ i  +  e _ 

cos (0  —  a  -  oj)  d*Yd*[cos (a+w)  +  sin(a-Ko)  'tan0] 


(8.18) 
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By  means  of  several  trials,  the  value  of  d  corresponding  to  the  mini¬ 
mum  ky  may  be  determined. 

The  equations  presented  by  Seed  and  Goodman  differ  from  those  above 
only  in  that  their  presentation  dealt  with  the  condition  of  horizontal  ac¬ 
celerations  only,  whereas  these  equations  were  developed  for  combined 
vertical  and  horizontal  accelerations.  The  vertical  component  of  accelera¬ 
tion  appears  in  Equations  8.16  and  8.18  in  the  form  of  the  angle  fa). 

The  influence  of  the  angle,  co,  may  be  seen  in  its  effect  upon  the 
value  of  ky  in  Equation  8.18.  The  influence  of  the  angle,  a),  upon  the 
first  term  of  that  equation,  and  upon  the  denominator  of  the  second  term 
is  illustrated  graphically  in  Figure  8-15.  It  is  apparent  from  this 
figure  that  ky  does  not  vary  greatly  as  a  function  of  the  angle,  oj ,  for 
the  soil  and  slope  characteristics  of  the  current  testing  program. 

For  these  slope  specimens,  ky  at  a  value  of  fa)  *  33.7°  will  be  ap¬ 
proximately  1.2  times  the  value  at  fa)  *  0.  This  fact  is  of  value  to  this 
analysis,  as  it  allows  the  use  of  the  simplified  analysis  technique  for¬ 
mulated  by  Seed  and  Goodman  to  predict  yield  accelerations  for  the  current 
testing  program.  This  eliminates  the  need  to  perform  the  trials  and 
iterations  as  described  earlier. 

Simplified  Analysis 

In  the  simplified  analysis  proposed  by  Seed  and  Goodman,  the  yield 
acceleration  is  calculated  as: 

ky  *  tan($  -  a  +  $„.)  (8.19) 

SL 


dHK »«•«**** 
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where  <j>  is  a  function  only  of  the  shear  strength  intercept,  s.,  and  the 

o  Xj  X 

slope  length,  l.  The  values  suggested  for  <<>  are  presented  in  Figure 

SL 

8-16.  For  the  current  testing  program,  is  approximately  equal  to  0.029 

psf  (from  Figures  8-7  and  8-8) ,  and  the  slopes  are  approximately  3  ft  long. 

As  an  example  of  how  this  simplified  procedure  may  be  used,  consider 

the  case  of  a  slope  with  a  =  33.6°  and  u  =  36.9°.  From  Figure  8-16,  with 

s.  -  0.03  and  t  *  3  ft,  4>  -  0.8®.  Now,  from  Equation  8.19, 

X  SL 

ky  =  tan(35.4°  -  33.6®  +  0.8°)  (8.20) 

=  tan(2.6°)  =  0.045 

This  value  of  ky  is  that  for  0)  =  0.  To  obtain  the  value  of  ky  when  u)  = 
36.9°,  reference  is  made  to  Figure  8-15.  The  value  of  the  first  term  in 
Equation  8.18,  illustrated  in  Figure  8-15 (b),  is  approximately  1.22  times 
the  value  when  u)  =  0.  Similarly,  the  value  of  the  second  term  in  Equation 
8.18,  whose  denominator  is  illustrated  in  Figure  8-15 (a),  is  also  approxi¬ 
mately  1.22  times  the  value  when  u>  =  0.  The  value  of  ky  from  Equation 
8-30  must  therefore  be  multiplied  by  1.22  to  obtain  the  value  when  oo  = 
36.9®.  This  new  value  is  as  follows: 

ky(w=36.9®)  =  (,045)  ’  (1,22)  =  -055  (8.21) 

Predicted  Yield  Accelerations 

This  computation  was  performed  for  the  various  test  conditions  en¬ 
countered  during  this  testing  series,  and  a  summary  of  the  predicted  yield 
acceleration  values  are  presented  in  Table  8-4. 


SHEAR  STRENGTH  INTERCEPT,  S.  -  LB  PER  SQ  FT 


FIGURE  8-16  RELATIONSHIP  BETWEEN  SLOPE  LENGTH, 

SHEAR  STRENGTH  INTERCEPT,  AND  ANGLE  <$SI 

(AFTER  SEED  AND  GOODMAN) 


TABLE  8-4 


Predicted  Values  of  Yield  Acceleration, 
ky,  for  Various  Test  Conditions 


Slope 

Angle,  a 

Angle  (d 

33.6° 

35.4° 

0° 

.045 

.016 

33.7° 

.054 

.019 

36.9° 

.055 

.020 

(TEST  8Y) 


31 


052 


018 


(TEST  8Z) 
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Test  Results 

As  discussed  earlier,  test  results  consisted  of  digitized  records  of 
accelerations  and  displacements  measured  during  the  course  of  dynamic 
loading.  For  each  test,  four  accelerometers  and  four  DCDT's  measured  the 
response  of  the  slope  specimen  and  test  box  to  the  base  excitation,  and 
numerous  accelerometers  mounted  upon  the  table  measured  the  average  hori¬ 
zontal  and  vertical  table  acceleration  during  loading.  These  digitized 
records  are  stored  on  magnetic  tape,  and  are  easily  processed  with  the  aid 
of  a  computer. 

Presented  in  Figures  8-17  and  8-18  are  several  typical  acceleration 
time  histories  which  have  been  plotted  from  these  dioitized  test  results. 
In  Figure  8-17  are  several  A  records  showing  both  table  and  top  of  slope 
accelerations  in  both  the  horizontal  and  vertical  directions.  Note  how 
closely  the  top  of  slope  acceleration  records  mirror  the  average  table 
acceleration  records.  Similarly,  Figure  8-18  includes  several  B  records 
showing  table  and  top  of  slope  accelerograms  for  horizontal  and  vertical 
loading.  Some  high-frequency  noise  appears  to  have  been  "picked-up"  in 
the  t  ,  of  slope  response  shown  in  these  results. 

The  plotted  time  histories  shown  in  these  and  the  following  figures 
were  drawn  by  the  CALCOMP  plotter  at  the  University  of  California  Com¬ 
puting  Center,  Berkeley.  A  complete  set  of  CALCOMP  plots  of  results  of 
this  testing  program  are  on  file  along  with  the  digitized  records  at  the 


Geotechnical  Engineering  Office  at  Berkeley. 
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IGURE  8-18  ACCELERATION  TEST  RECORDS  FOR  TESTS  7X  AND  8Y 
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Data  Reduction 

To  obtain  estimated  values  of  the  yield  acceleration  from  the  accel¬ 
eration  and  displacement  records,  it  is  necessary  to  carefully  study  the 
displacement  records  to  determine  the  time  in  the  record  when  yielding 
first  began. 

Figures  8-19  and  8-20  include  several  typical  acceleration  and  dis¬ 
placement  records  for  tests  1  and  4  respectively.  In  Figure  8-19  yielding 
appears  to  begin  at  approximately  2  seconds  into  DCDT  record  No.  2,  near 
the  top  of  the  slope,  and  some  time  later  in  DCDT  record  No.  3,  near  the 
bottom  of  the  slope.  At  2  seconds,  the  average  horizontal  table  accel¬ 
eration  is  changing  from  0  g  to  approximately  .085  g,  and  accelerometer 
No.  4  at  the  top  of  the  slope  is  changing  from  .015  g  to  approximately 
0.1  g.  Although  the  .015  g  acceleration  before  2  seconds  is  essentially 
all  high  frequency  noise,  it  is  being  applied  upon  the  sample  and  is 
therefore  a  lower  bound  on  the  yield  acceleration.  Because  yielding  was 
first  observed  at  2  seconds,  the  yield  acceleration  for  this  slope  must 
fall  somewhere  between  these  acceleration  values. 

In  Figure  8-20,  yielding  may  be  seen  at  2  seconds  into  DCDT  records 
1  and  4,  at  both  the  top  and  bottom  of  the  slope.  At  this  time  the  lower 
bound  values  of  the  yield  acceleration  are  0  g  and  .018  g,  and  the  upper 
bound  values  are  .106  g  and  .118  g  for  the  table  acceleration  and  top  of 
slope  acceleration  respectively.  These  accelerations  are  obtained  by 
noting  the  vertical  (av>  and  horizontal  (a^)  peak  acceleration  values  and, 
because  they  are  in  phase,  calculating  the  resultant  acceleration,  a',  as 
follows: 
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FIGURE  H-20  ACCELERATION  AND  DCDT  TEST  RECORDS  FOR  TEST  NO. 
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a'  =  /a£  +  a*  (8.22) 

This  value  of  a'  is  then  the  lower  or  upper  bound  on  ky  for  the  slope. 

In  the  example  shown  in  Figure  8-20,  the  peak  vertical  table  accel¬ 
eration  shown  after  2  seconds  is  approximately  0.056  g.  The  peak  horizon¬ 
tal  table  acceleration  occurring  at  the  same  time,  which  is  not  shown  in 
the  figure,  is  approximately  0.09  g.  The  resultant  acceleration  is: 

a'  =  / ( .056) 2  +  (.09) 2  =  0.106  (8.23) 

This  acceleration  is  an  upper  bound  on  the  yield  acceleration  measured 
from  the  bottom  of  the  slope. 

It  is  interesting  to  note  in  Figure  8-20  that  acceleration  record  No. 

3  shows  a  marked  shifting  of  center  acceleration  beginning  at  approximately 
seven  seconds  into  the  record.  This  occurred  because  the  soil  around  this 
accelerometer  began  to  move  significantly  at  this  time  and  the  "plane”  of 
the  accelerometer  rotated.  The  accelerometer  was  originally  oriented 
vertically,  so  that  the  center  axis  was  at  a  value  of  1  g.  If  the  center 
axis  moves  by  .2  g,  as  it  appears  to  at  eight  seconds,  then  the  new  center 
axis  becomes  0.8  g.  The  accelerometer  has  been  reoriented  by  the  angle  0 
from  the  vertical,  where  0  is  defined  by  the  following  expression 

0  =  arc  cos (0. 8)  =  36.9°  (8.24) 

Yield  Accelerations 

The  computations  described  above  for  determining  the  upper  and  lower 
bounds  on  the  yield  acceleration,  ky,  were  performed  for  all  tests  in 
this  testing  series.  The  results  are  tabulated  and  compared  with  the 
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predicted  values  of  ky  in  Table  8-5.  It  is  noteworthy  that  the  predicted 
values  of  ky  generally  fall  within  the  bounds  presented  in  that  table.  It 
should  also  be  noted  that  for  all  but  two  of  the  tests,  yielding  was  ob¬ 
served  within  2  to  3  seconds  into  the  record,  as  the  first  increment  of 
acceleration  was  applied. 

The  observed  values  of  ky  show  two  things: 

1.  There  is  often  a  significant  difference  between  the  amplitude  of 
the  accelerations  at  the  bottom  of  the  slope  (the  table)  and  the 
top  of  the  slope,  due  primarily  to  amplification  of  the  accelera¬ 
tions  ;  and 

2.  There  is  usually  a  significant  difference  between  the  observed 
upper  bound  and  observed  lower  bound,  due  primarily  to  the  fact 
that  the  increments  by  which  the  acceleration  amplitudes  were  in¬ 
creased  were  relatively  large:  typically  larger  than  ky  itself. 

This  latter  effect  wa3  minimized  somewhat  in  the  B  record  tests  because 
a  smaller  increment  of  increased  acceleration  amplitude  was  used. 

Due  to  these  facts  and  the  general  technical  difficulties  of  making 
accurate  measurements  of  accelerations,  it  is  difficult  to  obtain  highly 
precise  values  of  the  observed  ky  values.  However,  the  comparisons  shown 
in  Table  8-5  serve  at  least  to  support  the  conclusion  that  no  gross  er¬ 
rors  were  made  in  computing  the  predicted  ky  values.  Furthermore,  as  a 
corollary,  it  may  be  reasonable  to  conclude  that  ky  can  be  computed  as 
accurately  as  it  can  be  directly  measured. 

One  of  the  objectives  of  this  research  was  to  investigate  the 
validity  of  direct  vectorial  superposition  of  the  horizontal  and  vertical 
acceleration  in  arriving  at  the  yield  acceleration.  Because  of  the 


Time"  refers  to  the  time  in  a  test  record  where  yielding  was  first  observed. 
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difficulty  in  pinpointing  precise  values  of  the  observed  Icy  values,  the 
validity  of  this  superposition  cannot  be  verified  conclusively.  However, 
the  experimental  evidence  is  at  least  consistent  with  the  validity  of  this 
superposition. 

It  should  also  be  noted  that  the  calculated  (predicted)  value  of  ky 
is  quite  sensitive  to  the  near-surface  <p  value,  which  is  not  easily  mea¬ 
sured  with  great  accuracy.  In  view  of  this  fact,  it  appears  in  retrospect 
that  any  error  introduced  by  assuming  validity  of  vectorial  addition  of 
the  accelerations  is  much  smaller  than  the  error  introduced  by  using  im¬ 
precise  <j>  values  in  the  computations.  Furthermore,  it  appears  that  any 
error  incurred  by  direct  vectorial  superposition  is  smaller  than  that 
which  can  be  detected  using  an  experimental  program  similar  to  that  de¬ 
scribed  herein. 

In  view  of  these  conclusions  it  is  recommended  that  direct  super¬ 
position  of  simultaneous  horizontal  and  vertical  accelerations  be  made 
as  outlined  in  this  chapter  if  practical  cases  should  arise  where  yield 
accelerations  must  be  predicted. 

Qualitative  Observations 

In  addition  to  the  quantitative  tabulations  presented  in  Table  8-5, 
it  is  of  value  to  study  the  test  records  and  qualitatively  evaluate  the 
performance  of  the  slopes  during  dynamic  loading. 

The  four  DCDT  test  records  for  Test  No.  3  are  presented  in  Figure 
8-21.  Referring  to  the  positioning  of  these  instruments  upon  the  slope 
which  was  shown  in  Figure  8-1;  it  is  seen  that  DCDT's  1  and  2  are  lo¬ 
cated  in  the  upper  portion  of  the  slope,  and  DCDT's  3  and  4  are  in  the 


FIGURE  ft- 21  OCDT  TEST  RECORDS  FOR  TEST  NO. 
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lower  half.  DCDT's  1  and  3  are  vertically  oriented,  and  DCDT's  2  and  4 
are  horizontal. 

The  two  horizontal  DCDT's  reflect  the  acceleration  record  for  Test 
No.  3  which  was  shown  in  Figure  8-17.  Although  yielding  may  be  seen  from 
the  very  beginning  of  these  records ,  large  permanent  displacements  are  not 
seen  until  late  in  the  record.  This  is  because  yielding  may  occur  during 
only  a  small  portion  of  each  cycle  of  loading,  and  displacements  do  not 
accumulate  to  significant  values  until  the  larger  accelerations  occur. 

Shown  in  Figures  8-22  through  8-25  are  the  DCDT  test  records  for  the 
four  B  record  tests:  7X,  8X,  8Y,  and  8Z.  As  was  shown  in  Table  8-1,  the 
only  differences  between  these  tests  was  in  the  value  of  the  vertical  com¬ 
ponent  of  acceleration.  The  vertical  acceleration  was  zero  for  Test  7X, 
was  0.67  of  the  horizontal  for  Test  8X,  etc.  Because  the  horizontal  com¬ 
ponent  remained  unchanged,  the  resultant  acceleration  amplitude  steps 
were  different  for  the  different  tests. 

The  DCDT  test  records  for  Test  7X  shown  in  Figure  8-22  show  con¬ 
siderably  less  permanent  displacement  than  the  other  three  tests,  as  would 
be  expected.  The  DCDT  test  records  for  Tests  8X,  8Y,  and  8Z  show  remark¬ 
able  similarity,  both  in  their  shapes  and  in  the  amplitudes  of  the  dis¬ 
placements.  Theoretically,  the  displacements  should  be  slightly  greater 
for  Test  8Y,  and  less  for  Test  8Z,  since  the  vertical  accelerations  are 
slightly  higher  and  lower  for  those  tests.  It  is  difficult  to  distinguish 
significant  variation  for  these  records,  however. 

Additional  Observations 

It  should  be  emphasized  that  the  placement  of  DCDT's  on  the  face  of 
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FIGURE  8-23  DCDT  TEST  RECORDS  FOR  TEST  NO.  8X 
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the  test  slopes  was  intended  to  assist  in  determining  the  point  at  which 
yielding  began,  and  not  to  accurately  measure  displacements  or  strains. 
Because  of  the  manner  in  which  the  DCDT  core  probes  were  carefully  in¬ 
serted  into  the  soil,  minimizing  the  disturbance  to  the  soil  structure, 
they  are  only  capable  of  accurately  measuring  the  very  smallest  displace¬ 
ments.  As  the  sand  grains  become  mobilized,  they  quickly  flow  around  the 
core  probes,  thus  preventing  accurate  measurements. 

As  an  additional  tool  in  this  testing  series,  several  tests  were 
filmed  with  a  high-speed  data  camera  on  16  millimeter  film.  These  films 
were  synchronized  to  the  base  acceleration  record  and  contained  a  minimum 
of  40  frames  per  cycle  of  loading.  Thin  horizontal  lines  were  highlighted 
along  the  sides  of  those  slopes  which  were  photographed  to  provide  a  clear 
view  of  any  large  soil  displacements  and  sand  flows  during  dynamic  loading. 
Additionally,  a  triangular  "window"  was  cut  out  of  one  side  of  the  test 
box  to  provide  a  better  wide-angle  view  of  the  slope  samples.  These 
films  are  on  file  along  with  the  other  test  results  at  the  Geotechnical 
Engineering  Office  at  the  University  of  California,  Berkeley. 

Analysis  of  Real  Slopes 

The  analyses  of  real  slopes  of  cohesionless  materials  in  the  field, 
both  natural  and  man-made,  are  somewhat  more  complex  than  for  the  care¬ 
fully  controlled  laboratory  test  slopes  of  standard  sand  analyzed  in  this 
study.  As  discussed  earlier,  the  predicted  yield  accelerations  of  slopes 
are  extremely  sensitive  to  the  shear  strength  characteristics  of  the 
outermost  slope  material.  The  in-situ  outer  surface  shear  strengths  may 
vary  significantly  in  real  slopes,  and  will  be  very  difficult  to  accurately 
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measure.  Nevertheless,  it  is  believed  that  the  analysis  techniques  de¬ 
scribed  herein  may  be  used  with  appropriate  precautions  and  careful  judg¬ 
ment  to  predict  the  yielding  response  of  real  slopes  to  real  earthquakes. 

When  selecting  the  shear  strength  characteristics  of  the  outer  slope 
materials  for  analysis  of  a  real  slope,  care  should  be  taken  to  use  the 
in-situ  strengths.  The  shear  strength  intercept,  s^,  will  tend  to  be 
higher  in-situ,  particularly  for  older  slopes  or  those  with  surface  vege¬ 
tation.  The  primary  effect  of  a  higher  value  of  s^  is  to  move  the  criti¬ 
cal  sliding  surface  down  into  the  slope,  thus  increasing  the  yield  accel¬ 
eration. 

For  some  slopes,  special  direct  shear  tests,  similar  to  those  per¬ 
formed  in  this  study,  may  be  performed.  In  other  slopes  some  other  special 
field  test  may  be  devised  to  ascertain  the  in-situ  shear  strength  charac¬ 
teristics.  In  either  case,  careful  judgment  should  be  exercised  in 
selecting  the  values  of  s^  and  <f>  from  the  field  tests  to  insure  that  the 
lowest  measured  values  (and  thus  the  most  critical  zones)  are  given  maxi¬ 
mum  consideration  in  selecting  weighted  average  values  for  analysis  pur¬ 
poses. 


Displacements 


Seed  and  Goodman  (1964)  have  also  proposed  a  procedure  for  estimating 
total  permanent  displacements  in  slopes  from  yield  acceleration  values  and 
dynamic  acceleration  records.  This  procedure  is  based  on  suggestions  by 
Newmark  (1963) .  It  has  been  utilized  and  further  refined  by  Makdisi  and 
Seed  (1977)  and  others.  The  procedure,  simply  stated,  provides  for  the 


double  integration  of  dynamic  accelerations  in  excess  of  the  yield  accel- 
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eration  to  obtain  values  of  dynamic  displacement.  Because  displacement 
usually  occurs  during  only  one  half  of  each  dynamic  cycle  of  loading,  total 
residual  displacement  may  be  estimated  by  summing  the  incremental  dynamic 
displacements  accumulated  for  a  particular  earthquake  record. 

When  this  method  is  applied  to  real  slopes,  special  precautions  should 
be  taken  to  consider  the  possibility  that  the  predicted  yield  acceleration 
may  decrease  with  progressively  increasing  displacements.  Because  the  cal¬ 
culated  yield  acceleration  is  highly  sensitive  to  the  value  used  for  the 
shear  strength  intercept,  s^;  the  fact  that  the  true  value  of  s^  may  change 
with  straining  should  be  taken  into  account.  The  variation  of  s^  with 
strain  may  also  be  evaluated  from  a  field  testing  program. 

It  should  also  be  noted  that  for  the  tests  herein  the  horizontal  and 
vertical  accelerations  were  in  phase  and  at  the  same  frequency.  Thus  '*cd" 
was  constant  for  a  given  test.  By  contrast,  in  real  slopes  during  an 
actual  earthquake,  the  horizontal  and  vertical  accelerations  would  be  ex¬ 
pected  to  have  different  and  variable  frequencies.  Thus  a  numerical  pro¬ 
cedure  using  a  computer  program  must  be  used  to  perform  the  double  inte¬ 
gration  of  the  accelerations  in  excess  of  the  yield  acceleration.  This 
integration  must  be  performed  over  very  small  time  steps,  and  the  yield 
acceleration  must  be  recomputed  for  each  time  step,  as  a  function  of  oj. 

The  double  integration  process  should  be  repeated  for  various  rela¬ 
tive  "starting  points"  for  the  horizontal  and  vertical  acceleration  re¬ 
cords  in  order  to  obtain  the  range  of  residual  displacement  values  which 
corresponds  to  these  records. 
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Conclusions 

A  large  scale  slope  model  testing  program  was  developed  to  explore 
the  effects  of  combined  vertical  and  horizontal  dynamic  loading  on  the 
yield  acceleration  of  sand  slopes.  The  equipment,  sample  formation  and 
instrumentation  used  was  discussed  in  the  first  part  of  this  chapter  and 
in  Appendix  A- 3. 

The  strength  characteristics  of  the  Monterey  No.  0  sand  were  evalu¬ 
ated  with  special  direct  shear  tests  performed  on  the  test  slopes.  The 
results  of  this  evaluation  are  summarized  in  Figures  8-7  and  8-8. 

The  yield  accelerations  for  the  various  dynamic  test  conditions  were 
predicted  by  using  an  extension  of  an  analysis  method  originally  proposed 
by  Seed  and  Goodman  (1964).  Equations  8.16  and  8.18  were  derived  to  in¬ 
clude  the  angle,  oj,  which  permits  a  vertical  component  of  acceleration  to 
be  included  in  the  analysis.  By  modifying  the  simplified  technique  pro¬ 
posed  by  Goodman  and  Seed,  a  simplified  technique  which  includes  vertical 
as  well  as  horizontal  accelerations  was  developed.  From  those  equations 
the  yield  accelerations  in  Table  8-4  were  predicted. 

Using  a  procedure  described  in  this  chapter,  the  range  of  yield  ac¬ 
celerations  observed  during  the  various  model  tests  were  determined,  and 
were  compared  with  the  predicted  values  in  Table  8-5.  It  was  difficult 
to  draw  definitive  conclusions  concerning  these  comparisons  both  because 
of  the  variation  in  acceleration  amplitude  from  top  to  bottom  of  the 
slopes,  and  because  of  the  relatively  large  increments  by  which  the  table 
accelerations  were  increased. 

It  was  concluded  from  this  study  that,  because  of  these  limitations 
and  the  general  technical  difficulties  of  making  accurate  measurements  of 


accelerations,  it  is  difficult  to  obtain  highly  precise  values  of  observed 
yield  accelerations.  Nevertheless,  it  appears  that  the  yield  acceleration 
may  be  predicted  by  the  method  outlined  in  this  chapter,  at  least  as  ac¬ 
curately  as  it  may  be  directly  measured.  The  results  summarized  in  Table 
8-5  appear  to  support  the  conclusion  that  the  predicted  values  of  yield 
acceleration  are  reasonable  estimates  of  the  true  values. 

It  is  recommended  that  direct  superposition  of  simultaneous  horizon¬ 
tal  and  vertical  accelerations  be  made,  as  outlined  in  this  chapter,  if 
practical  cases  should  arise  where  yield  accelerations  must  be  predicted. 
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Chapter  9 

Summary  and  Conclusions 

Introduction 

The  objective  of  this  research  was  to  study  the  interaction  effects 
of  combined  compression  and  shear  loading  on  the  response  of  sands  to 
dynamic  loading.  Various  experimental  studies  were  performed  to  evaluate 
the  response  of  Monterey  No.  0  sand  to  various  types  of  combined  compres¬ 
sion  and  shear  cyclic  loading.  The  behavior  observed  in  these  tests  was 
then  used  to  evaluate  postulated  theories  of  soil  response  under  these 
combined  loading  conditions.  The  results  of  these  studies  and  the  main 
conclusions  are  summarized  in  this  Chapter. 

Effects  of  Combined  Dynamic  Loading 

During  the  course  of  this  research  significant  effects  were  observed 
during  combined  compression  and  shear  cyclic  loading.  The  primary  ob¬ 
served  effect  of  combined  loading  was  the  more  rapid  degradation  of  modu¬ 
lus  with  strain  than  would  otherwise  occur. 

Because  the  effect  of  primary  practical  interest  is  the  influence  of 
simultaneous  shear  straining  on  the  degradation  of  the  compression  modu¬ 
lus,  it  was  this  effect  that  was  evaluated  in  detail  in  this  study. 


Strain  Ratio  Method 


Two  methods  were  developed  and  presented  for  calculating  the  degra- 


dation  of  compression  modulus  with  strain  under  combined  loading  conditions 
The  first  of  these,  called  the  Strain  Ratio  Method,  requires  the  computa¬ 
tion  of  either  the  instantaneous  or  an  overall  average  ratio  of  shear 
strain  amplitude  to  compression  (normal)  strain  amplitude.  By  entering 
the  sets  of  curves  presented  in  Figure  7-7,  the  additional  degradation 
in  compression  modulus  due  to  the  presence  of  shear  strain  may  be  deter¬ 
mined.  This  method  requires  logarithmic  interpolation  among  the  curves, 
depending  on  the  value  of  normal  compression  strain  amplitude  and  mean 
confining  pressure. 

Octahedral  Shearing  Strain  Method 

A  simpler  method,  called  the  Octahedral  Shearing  Strain  Method,  re¬ 
quires  that  either  the  instantaneous  or  an  average  value  of  octahedral 
shearing  strain  be  calculated  from  Equation  2-42.  If  the  principal 
strains  are  known,  the  simpler  Equation  2-44  may  be  used  for  this  calcu¬ 
lation.  Once  the  octahedral  shearing  strain  is  known,  the  total  degra¬ 
dation  in  compression  modulus  may  be  determined  directly  from  Figure  7-6, 
interpolating  only  for  the  mean  confining  stress.  This  total  degradation 
in  modulus  includes  both  the  component  due  to  normal  compressive  straining 
and  that  due  to  the  presence  of  shear  strains. 

Using  either  of  these  methods  a  reasonable  estimate  of  the  degrada¬ 
tion  of  modulus  with  the  total  straining  may  be  obtained.  This  is  a 
significant  improvement  over  the  current  state  of  engineering  practice. 
However,  analysis  of  Figure  7-7  shows  that  the  degradation  in  modulus  due 
to  interaction  is  of  importance  only  when  the  shear  strain  is  significant 
compared  to  the  normal  strain.  In  other  cases  the  use  of  this  technique 
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provides  no  advantage  over  current  practice. 

Combined  Loading  Effects  on  Slope  Stability 

A  series  of  large-scale  shaking  table  tests  were  conducted  upon  slope 
models  constructed  of  sand.  These  samples  were  subjected  to  purely  hori¬ 
zontal  or  combined  horizontal  and  vertical  cyclic  accelerations  in  an 
effort  to  determine  the  point  at  which  yielding  in  the  slopes  began.  The 
results  of  these  tests  were  summarized  in  Table  8-5. 

Analysis  Method 

A  method  of  analysis  was  proposed  for  use  in  practical  cases  where 
yield  accelerations  for  slopes  under  combined  horizontal  and  vertical  ac¬ 
celerations  are  needed.  This  method  is  an  extension  of  a  method  originally 
developed  by  Seed  and  Goodman  (1964)  for  horizontal  excitation  alone. 

The  analysis  method  requires  the  calculation  of  a  yield  acceleration, 
ky,  from  Equation  8.18  for  the  values  of  plane  strain  friction  angle,  <p; 
slope  angle,  a;  angle  of  acceleration,  <d;  shear  strength  intercept,  ; 
the  length  of  the  slope,  L;  and  the  density  of  the  soil,  y^.  This  value 
of  ky  is  calculated  for  several  values  of  the  depth  of  sliding  surface, 
d;  and  for  the  minimum  value  of  the  internal  force,  P,  for  those  values 
of  d.  The  minimum  ky  is  determined  by  an  iterative  process,  utilizing 
Equations  8-16  and  8-18. 

A  modified,  simplified  approach  was  also  presented.  This  simplified 

method  requires  the  determination  of  the  factor,  4>  r  from  Figure  8-16, 

SL 

and  the  calculation  of  the  horizontal  ky  using  Equation  8-19.  This  value 
is  then  modified  to  account  for  the  acceleration  angle,  u>,  by  use  of 


Equation  8-18.  An  example  of  this  latter  process  was  presented  using  the 


shaking  table  strength  characteristics  and  slope  geometry,  with  the  effect 
of  variation  of  the  angle,  u>,  shown  in  Figure  8-15. 

This  latter,  simplified  method  was  used  for  the  slopes  and  strength 
characteristics  of  this  testing  program,  with  the  results  presented  in 
Tables  8-4  and  8-5. 

It  was  concluded  from  this  study  that,  because  of  the  limitations  of 
the  testing  program  and  the  general  technical  difficulties  of  making  ac¬ 
curate  measurements  of  accelerations,  it  is  difficult  to  obtain  highly 
precise  values  of  experimentally  determined  yield  accelerations.  Never¬ 
theless,  it  appears  that  the  yield  acceleration  may  be  predicted  by  the 
analytical  method  outlined  at  least  as  accurately  as  it  may  be  directly 
measured.  The  results  of  this  testing  program,  summarized  in  Table  8.5, 
support  the  conclusion  that  the  analysis  method  proposed  provides  reason¬ 
able  estimates  of  the  true  yield  acceleration  values. 

It  is  therefore  recommended  that  direct  superposition  of  simultaneous 
horizontal  and  vertical  accelerations  be  made,  as  outlined  in  Chapter  8, 
should  practical  cases  arise  where  yield  accelerations  for  slopes  must  be 
predicted . 
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Appendix  A 
Testing  Apparatus 

The  following  testing  apparatus  were  utilized  in  this  study: 
A-l  Resonant  Column-Torsional  Shear  Apparatus 

A-2  Hollow  Cylinder  Test  Apparatus 

A- 3  Test  Box  for  Slope  Model  Studies 
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Appendix  A-l 

Resonant  Column-Torsional  Shear  Apparatus 

Introduction 

The  Resonant  Column-Torsional  Shear  testing  apparatus  used  in  this 
study  was  developed  by  Professor  Vincent  P.  Drnevich  of  the  University  of 
Kentucky.  The  apparatus  is  capable  of  exciting  cylindrical  specimens 
either  vertically  or  torsionally,  or  in  both  modes  simultaneously.  The 
testing  set-up  and  a  close-up  of  a  sample  prepared  for  testing  are  shown 
in  Figure  Al-1. 

As  seen  in  Figure  Al-l(b),  vertical  excitation  is  obtained  by  ap¬ 
plying  an  A.C.  voltage  across  the  large  vertical  coil  mounted  atop  the 
specimen.  The  application  of  this  voltage  creates  a  sinusoidally  varying 
magnetic  field  which  interacts  with  the  field  of  the  large  permanent  mag¬ 
net  mounted  to  the  counterbalanced  framework,  resulting  in  a  sinusoidally 
varying  vertical  force  being  applied  upon  the  specimen. 

By  varying  the  frequency  and  amplitude  of  the  A.C.  voltage,  the 
cyclic  stress  applied  upon  the  specimen  can  be  controlled.  A  vertically 
oriented  accelerometer  is  mounted  within  the  cap  assembly  of  the  speci¬ 
men,  which  measures  the  acceleration  response  to  the  applied  cyclic  stress. 

Similarly,  torsional  stress  is  applied  to  the  specimen  by  applying 
an  A.C.  voltage  across  the  four  torsionally  oriented  coils  which  are 
wired  in  series  and  mounted  upon  the  counterbalanced  framework  alongside 
the  specimen.  The  resulting  electromagnetic  field  reacts  with  the  mag¬ 
netic  field  of  the  four  torsionally  oriented  permanent  magnets  affixed 


to  the  cap  assembly  of  the  specimen,  creating  a  cyclic  torsional  stress. 


FIGURE  Al-1  RESONANT  COLUMN  TESTING  APPARATUS  AND  SETUP 
CLOSE-UP  OF  SAMPLE  PREPARED  FOR  TESTING  (b) 
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A  torsionally  oriented  accelerometer  is  mounted  in  the  cap  assembly  for 
measuring  the  acceleration  response  to  the  applied  torsional  cyclic  stress. 
Because  the  two  loading  systems  are  mechanically  independent,  it  is  possible 
to  simultaneously  excite  the  specimen  vertically  and  torsionally  without 
mechanical  interaction  of  the  loading. 

Test  Set-up 

The  complete  test  set-up  is  shown  in  Figure  Al-l(a).  Just  below  the 
testing  apparatus  on  the  first  shelf  of  the  cabinet  is  a  1  ft  cube  con¬ 
crete  block.  This  block  acts  as  a  counter  balance  for  the  applied  cyclic 
loads,  and  is  securely  bolted  to  the  steel  framed  cabinet.  The  base  of 
the  testing  apparatus,  including  the  framework  alongside  the  specimen  and 
the  large  vertical  permanent  magnet  are  securely  attached  to  this  counter¬ 
balancing  concrete  block. 

The  testing  apparatus  is  completely  enclosed  within  an  airtight 
lucite  cell.  This  enclosure  allows  for  the  application  of  cell  confining 
pressures  up  to  approximately  7  KSC.  Additionally,  the  specimen  may  be 
confined  with  intercell  vacuum  of  less  than  1  KSC.  Although  the  equip¬ 
ment  is  designed  to  allow  both  sample  saturation  and  the  use  of  water  or 
other  fluid  for  applying  the  external  cell  confining  pressure,  these 
options  were  not  utilized  in  this  testing  series. 

The  A.C.  voltages  used  to  excite  the  specimens  in  both  the  vertical 
and  torsional  directions  were  produced  by  a  low-frequency  sine  wave 
generator  and  power  amplifier  combination.  The  two  generator/amplifier 
combinations  are  shown  in  the  right  side  of  Figure  Al-l(a),  above  the 
utility  cabinet.  The  sine  wave  generators  used  in  this  testing  series 


were  Hewlett-Packard  Model  202C,  and  the  power  amplifiers  were  Hewlett- 


Packard  Model  6824A. 

The  output  of  these  power  amplifiers  is  wired  to  a  control  box  which 
allows  the  operator  to  quickly  switch  the  voltage  to  the  coils  when  desired. 
This  control  box,  which  is  located  on  the  right  hand  side  of  the  second 
shelf  below  the  testing  apparatus  in  Figure  Al-l(a),  also  provides  a 
sampling  point  for  precise  measurement  of  the  amplitude  and  frequency  of 
the  voltages  applied  to  the  coils. 

Located  just  to  the  left  of  the  control  box  in  the  center  of  the 
second  shelf  is  a  Fluke  Model  1900a  frequency  counter  used  for  precise 
frequency  measurements.  The  applied  A.C.  voltage  is  displayed  on  the 
vertical  plates  of  an  oscilloscope;  in  this  study  a  Tectronics  Model  565 
Dual  Beam  Oscilloscope  was  used,  and  is  shown  on  the  left  side  in  Figure 
Al-l(a).  A  dual  beam  oscilloscope  was  helpful  in  this  testing  program 
because  it  allowed  for  simultaneous  monitoring  of  the  vertical  and  tor¬ 
sional  response  during  conbined  loading. 

The  two  accelerometers  mounted  in  the  cap  assembly  are  Colombia 
electrolytic  devices,  Model  200-1-H,  and  require  a  charge  amplifier  or 
cathode  follower  to  condition  their  outputs  before  they  can  be  read  by 
conventional  high  impedance  measuring  equipment.  Two  charge  amplifiers 
were  used  to  supply  an  acceleration  response  signal  to  the  oscilloscope. 
These  charge  amplifiers,  which  are  shown  on  the  left  side  of  the  second 
shelf  in  Figure  Al-l(a),  are  Colombia  Model  4102. 

The  output  of  the  charge  amplifiers  are  applied  to  the  horizontal 
plates  of  the  oscilloscope  to  produce  Lissajous  figures  with  the  voltages 


applied  to  the  driving  coils.  The  acceleration  amplitude  may  also  be 


measured  using  the  oscilloscope  display. 

Cylindrically  shaped  samples  are  formed  with  a  rubber  membrane  and 
conventional  forming  molds.  The  vertical  weight  of  the  cap  assembly  is 
offset  with  a  constant-force  spring  shown  in  the  top  of  Figure  Al-l(b). 

Determination  of  the  strain  amplitudes  and  damping  factors  from  the 
raw  test  data  is  described  in  Chapter  4;  and  the  determination  of  dynamic 
compression  and  shear  moduli  is  discussed  in  Appendix  C-3.  Several  typi¬ 
cal  data  sheets  are  shown  in  Figures  Al-2  through  Al-4,  and  the  test  re¬ 
sults  from  this  raw  data  are  presented  in  Appendices  Bl  through  B3. 


208 


BKIVERSITX  OF  CAlIFC.wA 
RESOIIANT  COL'JIIIJ  -  TCHSIOIJAL  S!EAR 
DATA  SRSST 

Date  £d2±J31 

Operator  f“c _  Test  To.  iifb-Z  Page  _S_  0f  3 


Vfbl  ______  WbO  _  Spill 


FIGURE  A1 -3(c)  DATA  SHEET  FOR  RESONANT  COLUMN 
TEST  NO.  12B-2  (continued) 


I 


210 


Appendix  A- 2 

Hollow  Cylinder  Test  Apparatus 

The  thin-walled  hollow  cylinder  testing  apparatus  used  in  this  testing 
program  was  a  modification  of  a  static  loading  device  originally  developed 
by  Professor  Poul  V.  Lade  of  the  University  of  California  at  Los  Angeles 
and  Professor  J.  M.  Duncan  of  the  University  of  California  at  Berkeley 
(Lade,  1972) .  A  front  view  of  the  testing  apparatus  showing  the  control 
panel  is  shown  in  Figure  A2-1,  and  a  diagram  of  the  main  loading  chain  is 
shown  in  Figure  A2-2.  The  various  numbered  components  from  the  load  chain 
diagram  are  listed  in  Table  A2-1. 

This  testing  apparatus  is  capable  of  cyclic,  stress-controlled  loading 
in  both  the  vertical  and  torsional  directions,  both  separately  and  simul¬ 
taneously.  When  vertical  and  torsional  loading  are  applied  simultaneously, 
both  cyclic  excitations  are  at  the  same  frequency,  but  may  be  out  of  phase. 
The  apparatus  is  capable  of  operating  at  any  frequency  between  approxi¬ 
mately  0.2  Hertz  and  20  Hertz.  All  tests  in  this  testing  series  were  per¬ 
formed  at  a  frequency  of  approximately  0.3  Hertz. 

This  apparatus  is  capable  of  loading  hollow  cylinder  specimens 
cyclically  at  a  peak-to-peak  stress  of  up  to  100  psi  vertically  and  tor- 
sionally,  simultaneously. 


Loading  System 

Simply  stated,  dynamic  load  is  applied  to  the  test  specimen  as  follows 
1.  A  variable  speed  motor,  controlled  on  the  front  panel,  operates 
at  the  desired  frequency  of  loading  (i.e.  at  0.3  Hertz). 


2.  Two  offset  cam  assemblies  which  are  attached  to  the  shaft  of 


FIGURE  A2-2  DETAIL  OF  LOAD  CHAIN  OF  HOLLOW 
CYLINDER  TESTING  APPARATUS 
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TABLE  A2-1 

Components  of  Hollow  Cylinder  Testing  Apparatus 

1.  Proximeter  Target 

2.  Center  Shaft  -  Sample  Top  Clamping  Assembly 

3.  Sample  Top  Piece 

4.  Sample  Cap 

5.  Sample  In  Rubber  Membrane 

6.  Center  Shaft 

7.  Bellofram  Rubber  Cell  Seal 

8.  Sample  Base 

9.  Sample  Bottom  Plate 

10.  Top  of  Cabinet 

11.  Upper  Standard  Framework 

12.  Center  Shaft  Bearing  Holders  (2) 

13.  Torsional  Moment  Arms 

14.  Center  Standard  Plate 

15.  Lower  Standard  Framework 

16.  Vertical  Load  Cell 

17.  Thrust-Tension  Joint 

18.  Double  Ball  Joint 

19.  Spacer  Shaft 

20.  Spacer  Shaft 

21.  Center  Shaft  of  Double  Acting  Cylinder  (Thompson  Ball  Shaft) 

22.  Top  and  Bottom  Plates  of  Double  Acting  Cylinder  (2) 

23.  Top  and  Bottom  Walls  of  Double  Acting  Cylinder  (2) 
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24.  Center  Wall  of  Double  Acting  Cylinder 

25.  Piston  of  Double  Acting  Cylinder 

26.  Piston  Face  Plates  of  Double  Acting  Cylinder  (2) 

27.  Bellofram  Rubber  Membranes  of  Double  Acting  Cylinder  (2) 

28.  Rotolin  Center  Shaft  Bearings  (2) 

29.  Thompson  Ball  Bearings  (2) 


i 

■i 
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the  variable  speed  motor  cause  two  pistons  to  traverse  into  and 
out  of  two  small  pressure  regulators  in  a  cyclic  manner,  at  the 
same  frequency  as  the  motor. 

3.  The  output  of  the  two  small  pressure  regulators  are  fed  to  the 
control  of  two  large  pneumatic  volume  booster  relays  which  are 
fed  with  100  psi  air  pressure  and  1/2-inch  supply  lines. 

4.  The  output  of  these  pneumatic  volume  booster  relays  are  fed  into 
two  large  reservoir  cylinders  which  contain  an  air-oil  interface 
and  a  supply  of  low-viscosity  loading  oil. 

5.  The  output  of  one  of  the  two  large  cylinders  is  then  fed  directly 
to  the  lower  half  of  the  vertical  double  acting  cylinder,  while 
the  output  of  the  other  cylinder  is  fed  directly  to  one  of  two 
torsionally  balanced  torsional  loading  cylinders. 

6.  The  cyclic  pressure  in  the  lower  half  of  the  vertical  double 
acting  cylinder  reacts  with  the  constant  pressure  in  the  upper 
half  of  that  cylinder,  resulting  in  a  sinusoidally  varying  force 
which  is  transmitted  through  the  center  shaft  to  the  sample  top 
piece,  and  thus  to  the  sample. 

7.  The  cyclic  pressure  in  one  of  the  two  torsional  loading  cylinders 
reacts  with  the  constant  pressure  in  the  other  cylinder,  re¬ 
sulting  in  a  torsional,  sinusoidally  varying  stress  which  is 
transmitted  through  the  torsional  moment  arms  and  the  center 
shaft  to  the  sample  top  piece,  and  thus  to  the  sample. 

A  photograph  of  the  left  side  of  the  testing  apparatus  showing  the 
variable  speed  motor,  offset  cams,  small  pressure  regulators,  and  the 
large  pneumatic  volume  booster  relays  is  shown  in  Figure  A2-3.  The  speed 
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control  for  the  variable  speed  motor  is  visible  in  the  lower  right  corner 
of  the  front  control  panel  in  Figure  A2-1. 

Careful  examination  of  Figure  A2-3(a)  shows  that  the  two  pneumatic 
volume  booster  relays,  which  are  mounted  on  the  cross-brace  in  the  center 
of  the  photograph  are  vented  at  the  control  port  at  their  tops  by  muffled, 
pin  type  vent  values.  This  is  done  because  the  small  pressure  regulators 
are  not  very  efficient  at  venting  air  pressure  when  their  regulated  pres¬ 
sure  is  dropped,  and  they  need  additional  externally  provided  venting  to 
keep  their  time  constants  low  during  that  portion  of  the  loading  cycle. 

There  are  several  details  of  interest  in  Figure  A2-3(b).  Close  ex¬ 
amination  of  the  shaft  couplers  used  between  the  center  shaft  of  the 
variable  speed  motor  and  the  offset  cams  reveals  a  small  micro-switch 
mounted  on  the  upper  coupling.  This  microswitch  was  installed  to  provide 
a  simple  means  of  insuring  that  cyclic  loading  begins  and  ends  at  the  same 
place  in  the  loading  cycle  for  each  test.  The  switch  is  normally  adjusted 
so  that  loading  begins  and  ends  at  the  balance  point  at  the  center  of  the 
sine  wave  of  loading. 

The  phase  lag  between  vertical  and  torsional  loading  is  easily  ad¬ 
justed  by  adjusting  the  relative  rotation  of  the  offset  cams  with  the 
shaft  couplings.  The  offset  cams  are  shown  precisely  in  phase  in  the 
photograph. 

The  balance  pressure  at  the  center  of  the  sine  wave  of  loading  is 
adjusted  by  positioning  of  the  shaft  between  the  offset  cam  and  the  piston 
of  the  small  pressure  regulator.  Lengthening  the  shaft  will  result  in  a 
higher  balance  pressure.  Care  must  be  used  when  this  adjustment  is  made 
that  the  selected  balance  pressure  will  allow  sufficient  "range"  of  peak- 
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to-peak  pressure  to  accommodate  the  range  of  load  anticipated  in  each  test. 

A  photograph  showing  the  rear  of  the  testing  apparatus  is  provided  in 
Figure  A2-4(a);  and  a  photograph  showing  the  right  side  and  the  two  tor¬ 
sional  loading  cylinders  is  presented  as  Figure  A2-4 (b) .  A  3/4-inch  pres¬ 
sure  hose  may  be  seen  at  the  rear  of  the  test  cabinet.  This  hose  provides 
air  pressure  for  the  loading  system  in  this  apparatus. 

The  two  large  reservoir  cylinders,  containing  the  air-oil  interface, 
are  in  the  bottom  of  the  cabinet,  out  of  the  view  of  the  photographs. 

Each  cylinder  is  approximately  seven  inches  in  diameter  and  approximately 
fourteen  inches  high.  It  should  be  noted  at  this  time  that  the  use  of 
air-oil  cylinders  and  oil  in  the  loading  cylinders  was  provided  in  this 
testing  apparatus  so  that  it  could  also  be  used  for  strain-controlled 
static  testing.  If  the  apparatus  was  intended  for  use  only  for  stress- 
controlled  testing,  it  would  have  been  quite  possible  (and  desirable)  to 
feed  the  output  of  the  pneumatic  volume  booster  relays  directly  into  the 
loading  cylinders,  using  air  only. 

The  vertical  double  acting  cylinder  is  visible  in  the  center  of 
Figure  A2-4(a)  near  the  bottom.  A  small  differential  pressure  trans¬ 
ducer  is  visible  at  the  top  of  this  cylinder.  An  identical  differential 
pressure  transducer  is  visible  on  the  under-side  of  the  right  torsional 
loading  cylinder  in  Figure  A2-4(b).  It  is  these  two  transducers  which 
are  calibrated  to  provide  a  measure  of  the  vertical  and  torsional  stress 
applied  upon  the  specimens  during  testing. 

Torsional  strain  is  measured  with  a  small  LVDT  mounted  with  long 
hose-clamps  directly  upon  the  sample  (not  visible  in  photographs) .  The 
vertical  strain  is  measured  with  a  proximeter  device,  which  may  be  seen 
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on  the  top  of  the  outer  pressure  cell  in  Figure  A2-4(b). 

External  cell  pressure  is  measured  with  the  pressure  gage  mounted 
directly  behind  the  proximeter  on  the  top  of  the  pressure  cell. 

The  levels  of  constant  pressure  are  maintained  in  the  upper  half  of 
the  vertical  double  acting  cylinder  and  in  one  of  the  torsional  loading 
cylinders  by  means  of  pressure  regulators  controlled  on  the  front  control 
panel.  The  ensure  that  constant  pressure  is  maintained  during  loading, 
two  reservoir  cylinders  are  provided  in  series  with  the  pressure  lines. 
These  cylinders  are  visible  at  the  bottom  of  the  cabinet  in  Figure  A2-4(a) 
on  the  two  sides  of  the  cabinet.  It  should  be  noted  that  two  pneumatic 
volume  booster  relays  would  have  served  equally  well  in  this  position. 

The  vertical  load  cell  shown  in  Figure  A2-2  is  used  primarily  for 
calibration  purposes,  for  adjustment  purposes  and  to  estimate  the  magni¬ 
tude  of  mechanical  friction  developed  in  the  load  chain.  The  thrust- 
tension  joint  and  the  double  ball  joint,  which  are  also  shown  in  that 
figure,  are  designed  to  de-couple  the  vertical  and  torsional  loading  ap¬ 
plied  to  the  center  shaft.  The  double  ball  joint  corrects  any  small  mis¬ 
alignment  between  the  vertical  double  acting  cylinder  and  the  center 
shaft,  and  the  thrust-tension  joint  allows  transfer  of  vertical  normal 
stresses  to  the  center  shaft,  but  rotates  freely,  preventing  any  tor¬ 
sional  stresses  from  transmitting  to  the  vertical  cylinder. 

The  two  torsional  loading  cylinders  transfer  their  load  to  the  tor¬ 
sional  moment  arms  by  means  of  1/4-inch  cables  18-inches  long.  Resultant 
lateral  forces  from  these  cables  are  balanced  out  by  the  two  large  Rotolin 
bearings  through  which  the  center  shaft  passes. 

Recording  of  vertical  stress  and  strain  during  this  testing  program 
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FIGURE  A2-5  TEST  RECORDS  FOR  4  TYPICAL  HOLLOW  CYLINDER  TESTS 


APPENDIX  A3 


Test  Box  for  Slope  Model  Studies 

The  slope  model  studies  were  performed  in  a  wooden  test  box  affixed 
to  the  20  ft  x  20  ft  shaking  table  at  the  Earthquake  Simulator  Laboratory, 
Earthquake  Engineering  Research  Center,  University  of  California,  Berkeley, 
California.  The  shaking  table  is  capable  of  reproducing  with  excellent 
accuracy  dynamic  accelerations,  velocities,  or  displacements  stored  on 
magnetic  tape  or  disk.  The  test  box  was  designed  both  to  contain  the 
sand  test  slopes  and  to  efficiently  transfer  the  dynamic  table  motions 
to  the  slope  models. 

A  schematic  diagram  of  a  typical  slope  specimen  is  shown  in  Figure 
A3-1.  The  inner  dimensions  of  the  test  box,  as  illustrated  in  this  figure, 
are  approximately  84.2  in  long  x  42.5  in  wide  x  21.7  in  high.  A  spacer 
board  has  been  added  as  shown  in  the  figure  to  reduce  the  total  volume  of 
sand  needed  for  the  construction  of  a  test  specimen.  This  spacer  board, 
which  is  positioned  at  an  angle  of  approximately  45°  with  the  horizontal, 
is  approximately  15.8  in  on  a  side.  Also  visible  in  Figure  A3-1  is  a  2  in 
steel  angle  reinforcing  system  on  the  outer  edges  of  the  test  box.  This 
reinforcing  system  provides  rigidity  to  the  outer  walls  of  the  test  box. 

A  view  of  the  test  box  affixed  to  the  shaking  table  is  shown  in 
Figure  A3-2(a).  As  the  box  was  placed  upon  the  table,  a  layer  of  high- 
strength  "hydrostone"  cement  mortar  was  placed  between  the  box  and  the 
table.  The  box  was  then  firmly  bolted  down  to  the  table  by  seven  1-1/2 
in  diameter  steel  bolts.  Three  of  these  bolts  were  tightened  against  a 
1/4  in  thick  steel  plate  which  was  placed  in  the  bottom  of  the  box,  and 
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the  remaining  four  bolts  were  tightened  against  the  reinforcing  steel 
angles  on  the  outer  corners  of  the  test  box.  two  of  these  outside  bolts 
are  visible  on  the  outer  edge  of  the  test  box  in  Figure  A3-2(a),  and  the 
steel  plate  and  bolts  in  the  bottom  of  the  test  box  may  be  seen  in  Figure 
A3-2 (b) . 

Behind  the  spacer  board  a  stiff  bracing  system  has  been  constructed 
to  ensure  that  the  deflection  of  the  board  during  loading  would  be  very 
small  compared  with  the  displacements  of  the  soil  slope.  This  bracing 
system  consisted  of  six  triangular  1-1/2  in  thick  plywood  wedges,  cut  to 
exactly  fit  between  the  board  and  the  inside  corner  of  the  test  box. 

Before  slope  specimens  were  constructed  in  the  test  box,  a  2  in  thick 
layer  of  high-strength  "hydrostone "  cement  mortar  was  placed  in  the  bottom 
of  the  box  on  top  of  the  steel  plate.  The  "hydrostone”  was  carefully 
leveled  and  a  thin  layer  of  the  Monterey  No.  0  test  sand  was  cemented  to 
the  top  of  it.  Also,  a  thin  layer  of  the  test  sand  was  epoxyed  with  a 
high  strength  resin  epoxy  to  the  spacer  board  and  the  back  wall  of  the 
test  box.  This  treatment  insured  a  level  bottom  and  good  soil-box  contact 
between  soil  slope  specimens  and  the  test  box. 

A  photograph  showing  the  specimen  forming  "mold"  in  position  within 
the  test  box  is  shown  in  Figure  A3-3(a).  The  mold  assembly  consists  of 
a  7  ft  x  4  ft  plywood  panel,  1  in  thick,  which  is  reinforced  with  three 
longitudinal  3  in  steel  angles  on  the  back  side,  and  associated  mounting 
hardware.  The  plywood  panel  is  positioned  to  the  desired  specimen  slope 
angle,  and  rigidly  secured  in  place  as  shown  in  Figure  A3-3(a).  After 
the  specimens  are  formed  within  the  mold  by  table  vibration,  the  plywood 
panel  is  carefully  removed  with  an  overhead  crane,  leaving  a  smooth  slope 
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face  on  the  specimen. 

The  mounting  hardware  which  secures  the  plywood  panel  was  designed 
to  provide  rigid  support  without  distorting  the  smooth  "plane"  inner  sur¬ 
face  of  the  plywood.  Two  steel  channels  along  the  sides  of  the  test  box, 
covered  with  strips  of  compressible  foam  rubber,  provide  upward  support 
at  the  edges  for  the  plywood  panel,  while  a  longitudinal  steel  angle  and 
two  2  in  x  4  in  lengths  of  lumber  provide  the  securing  downward  force. 

A  7  ft  long  lumber  2x4  provides  the  toe  footing  support  at  the  bottom  of 
the  plywood  panel. 

Instrumentation  is  mounted  both  within  the  test  slopes  and  on  the 
test  box  in  preparation  for  each  test.  As  the  slope  specimens  are  formed 
two  very  small  accelerometers,  one  vertical  and  one  horizontal,  are  care¬ 
fully  positioned  within  the  sand  near  the  top  face  of  the  slope.  The 
accelerometers  are  mounted  on  1/32  in  thick  perforated  aluminum  wafers, 
approximately  7/8  in  on  a  side.  This  allows  for  relatively  accurate 
orientation  with  respect  to  the  horizontal  and  vertical  directions.  Also 
two  accelerometers  are  mounted  on  the  back  side  of  the  test  box,  behind 
the  top  of  the  slope  on  the  steel  reinforcing  angle.  These  two  accelero¬ 
meters  are  also  oriented  for  the  vertical  and  horizontal  directions,  and 
may  be  observed  in  Figure  A3-1. 

A  cross-brace  is  placed  across  the  center  of  the  test  box  after  the 
test  specimen  is  formed,  and  four  DCDT's  are  affixed  to  it  on  a  mounting 
framework  as  shown  in  Figure  A3-3(b).  The  slope  pictured  in  this  figure 
has  already  failed,  and  characteristic  sloughing  of  the  sand  downslope 
may  be  observed.  Note  that  the  rod  extensions  of  the  DCDT  cores  are  bent 
slightly  as  a  result  of  the  soil  mobility.  This  effect  could  result  in 


inaccurate  readings  at  significant  strains,  as  discussed  in  this  report. 

One  additional  modification  performed  on  the  test  box  during  this 
testing  series  was  the  cutting  of  a  triangular  shaped  "window"  on  one  side 
of  the  box.  This  was  to  permit  the  photographing  of  the  slope  face  with 
a  high  speed  data  camera  during  testing.  The  opening  was  approximately 
10  in  x  10  in  with  a  diagonal  at  45°  with  the  horizontal.  This  modifica¬ 
tion  was  necessary  for  low  angle  photography,  nearly  parallel  with  the 


Example  Test  Results 


The  following  examples  are  included  of  test  results  from  computer 
analysis  of  the  raw  test  data: 

Resonant  Column  Test  Results 
B-l  Test  No.  1A-3 

B-2  Test  No.  12B-2 

B-3  Test  No.  18C-3 

Hollow  Cylinder  Test  Results 

B-4  Test  No.  A26L-2 

B-5  Test  No.  B29-3 

B-6  Test  No.  B31-1 

The  complete  test  results  are  available  in  Griffin  (1980). 
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APPENDIX  B-5 

Hollow  Cylinder  Test  No.  B29-3 
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Appendix  C 
Derivations 


The  following  derivations  were  made  for  use  in  this  study: 
C-l  Strain  Tensor  for  Resonant  Column  Specimen 
C-2  Strain  Tensor  for  Hollow  Cylinder  Specimen 
C-3  Calculation  of  Moduli  for  Resonant  Column  Tests 
C-4  Principal  Stresses  for  Resonant  Column  Specimen 
0-5  Principal  Stresses  for  Hollow  Cylinder  Specimen 
C-6  Calculation  of  Compression  Modulus  for  Hollow  Cylinder  Tests 
C-7  Calculation  of  Volume  of  Hollow  Cylinder  Specimen 
C-8  Calculation  of  Relative  Density 
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Appendix  C-l 

Strain  Tensor  for  Resonant  Column  Specimen 


From  the  expressions  for  the  strains  contained  in  Equations  2.5 
through  2.11  and  from  the  boundary  conditions  presented  in  Equations  2.12 
and  2.13,  the  following  expressions  may  be  written: 


and  u  is  unknown, 


£e 

£R0 

ERZ 

£9i 

but 


v  3v 

because  gg-  *  0 


3u 


+  zA0  -  zA0  ) 
o  o/ 

because  =  0 


1  3u 

2  r30 


(C-l.l) 

(C-l. 2) 

(C-l. 3) 

(C-l. 4) 


u 

From  Hooke's  Law, 

e 

x 

and  e 

x 

Also, 


'  er 


it  is  known  that: 


e  *  -ye 

y  z 


+  e  +  e  -  e  , 

y  z  vol 


(1  -  2y)e 

Z 


(C-l. 5) 


(C-l. 6) 

(C-l. 7) 


e  ,  ■  e_  +  e0  +  e  -  (1  -  2y)e  (c-l. 8) 

vol  R  0  z  z 

or  -2ye  (C-l.  9) 

R  v  Z 

Now,  if  one  considers  the  new  volume,  V..,  and  relates  it  to  the  old 

N 

volume,  V  ,  it  can  be  stated  that: 
o 
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vM  =  v  <1  +  e  . )  =  v  [1  +  (l  -  2y)e  ] 

N  o  vol  o  z 


(C-1.10) 


But  is  also  equal  to: 


V  *  H  (1  +  £  )  •  TIR, 
NO  Z  N 


(C-l.ll) 


where  Ho  is  the  old  height  and  R^  the  new  radius.  We  also  know  that  v 

is  a  function  of  H  and  R  ,  the  old  radius ,  as  follows : 
o  o 


V  *  H  *  TTR 
o  o  o 


(C-1.12) 


Combining  the  last  three  equations,  C-1.10  through  C-1.12,  the 


following  expression  for  R^  may  be  developed: 

Tl  +  (1  -  2y)el1/2 

>‘L  j 


R  *  R 
N  O 


(C-1.13) 


and  it  is  also  known  that: 


N- V1  +  ,k> 


(C-1.14) 


It  is  now  possible  to  write  an  expression  for  eR  as  follows: 


tl  +  (1  -  2y)ez' 
l+V  1 


1/2 


-  1 


(C-1.15) 


Now,  from  Equation  C-1.5: 


u  =  /  eR  dr 


or 


— (p^r-) 


(C-1.5) 


(C-1.16) 


It  is  possible  now  to  write  expressions  for  £g,  £rq,  and  and  as 


follows : 


I 


It  is  also  known  that: 

E  +  EQ  =  -2 y£  (C-1.19) 

R  o  Z 

so  that  one  may  write: 

e  =  eQ  =  -ye  (C-1.20) 

R  o  Z 

The  strain  tensor  may  now  be  written : 

-ye  0  0 

z 

e  =  0  -ye  k  rAe  (2.14) 

z  2  o 


It  should  be  further  noted  that  the  Cgz  component  is  a  function 

of  the  radius  r.  For  practical  application  purposes,  the  average  £gz 

2 

will  be  assumed  to  exist  at  r  =>  y  R,  and  be  considered  a  constant  at 
that  value. 
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Appendix  C-2 

Strain  Tensor  for  Hollow  Cylinder  Specimen 


From  the  expressions  for  the  strain  contained  in  Equations  2.5,  2.8 
through  2.11,  2.15,  and  2.16,  and  from  the  boundary  conditions  expressed 
in  Equations  2.12  and  2.13,  the  following  expressions  may  be  written: 


eR0  =  2  (r30  +  A0o2 


Ae  z)  = 

o  /  2  r30 


(C-2. 1) 


1  3u  .  3w 

eRZ  -  2  SI  '  because  37  =  0 


(C-2. 2) 


_  1  3v  1  _  aA  ,  3w 

=  —  v—  =  —  R  ,  where  =  0 

0z  2  3z  2  avg  o  30 


(C-2. 3) 


and  u  is  unknown,  but 


u  =  /  eR  dr 


(C-2. 4) 


If  it  is  assumed  that  Hooke's  Law  applies,  the  following  expressions 
may  be  written: 


1  U  p 

e  =  -  a  -  =•  a  -  =•  oQ 

R  E  R  E  z  E  0 


(C-2. 5) 


and 


1  u  u 

e  =  —  a  -  —  a  -  —  oa 

Z  E  Z  E  R  E  0 


.1  y  y 

e0  =  E  a0  ~  E  °R  "  E  az 


(C-2. 6) 

(C-2. 7) 


If  Eg  *  0,  Equation  C-2. 7  may  be  rewritten  as  follows: 


a0«y(oR  +  oz) 


(C-2. 9) 


Now  e  and  e  may  be  written: 

Z  R 
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e  =  ^  a  -  ^  a  -  jj--  [y  (a  +  a  )  ] 

R  E  R  E  2  E  R  z 


(C-2.9) 


e„  -  ^  a  -  ^  or  -  £*(u(a_  +  a  >  ] 
Z  £2  er  E  R  2 


(C-2. 10) 


Now,  the  change  in  the  radial  strain,  Ae  ,  during  loading  may  be  written: 

R 


(C-2. 11) 


And  because  Aa  *  0  in  this  testing  program, 

K 


(C-2. 12) 


Similarly, 


*.  •  w-y 


(C-2. 13) 


Combining  the  above  two  equations,  C-2. 12  and  C-2. 13  gives: 


•(rf)' 


Now,  because 


(C-2. 14) 


u  «  /  £  dr 
R 


(C-2. 4) 


for  Ravg  <  r  <  R2, 


"L  '(rf>4s  dr 

avg  H 

■  -  RW 


(C-2. 15) 


■  f^b 


-  W 


The  expressions  for  £Rg  and  may  now  be  written: 


e  »  I  *  o 
eR0  2  r36 


_  1  3u 

£  “  — —  St  0 

Rz  2  3z 


The  strain  tensor  may  now  be  written: 


1  R  A0 

2  avg  o 


4  R  A0 

2  avg  o 


(C-2.17) 

(C-2.18) 


(2.17) 


It  is  also  worth  noting  that  for  the  special  condition  where 
*  0,  Equation  (2.17)  is  transformed  so  that  the  values  of  £R,  Eg, 
and  £z  are  also  the  principal  strain  values,  e^,  and  respectively. 
In  this  special  case  the  strain  tensor  is  that  of  a  plane  strain 


condition. 
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Appendix  03 


Calculation  of  Dynamic  Moduli  for  Resonant  Column  Tests 


Compression  Modulus 


The  basic  expression  for  the  dynamic  compression  modulus,  E,  is  as 


follows : 


p  *  V‘ 


(C-3.1) 


where  p  is  the  mass  density  of  the  specimen  material,  and  where  V  is  the 

P 

"P-wave"  (compression  wave)  velocity  of  propagation.  For  the  solid 
cylindrical  specimen. 


(C-3.2) 


where  W  is  the  specimen  weight,  V  is  the  specimen  volume,  and  g  is  the  ac¬ 
celeration  of  gravity.  The  value  of  may  be  calculated  from  the  fol¬ 
lowing  expression: 


2-ir-f  -L 
v 


(C-3.3) 


where  f^  is  the  resonant  frequency  in  compression,  I.  is  the  length  (or 
height)  of  the  cylindrical  specimen  and  where  is  the  root  of  the  fre¬ 


quency  equation: 


\j>  *  tan  il) 
rv  v 


(C-3.4) 


where  WTQp  is  the  top  cap  weight.  The  value  of  in  this  testing 

series  is  1576.6  gm. 

During  the  resonant  column  test,  the  resonant  frequency  is  influenced 
by  dynamic  coupling  between  the  sample  and  the  cap  system.  If  a  test 
could  be  conducted  under  ideal  conditions,  where  the  top  cap  weight  was 


4 
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negligible  when  compared  with  the  specimen  weight,  there  would  be  no 
coupling  effect  and  f  would  reflect  the  true  resonant  frequency  of  the 
soil  specimen.  In  this  special  case,  *  tt/2,  and  Equation  C-3.3  reduces 
to  the  following: 

V  -  4  •  f  •  L  (C-3.5) 

P  v 

In  the  more  general  case,  the  dynamic  compression  modulus  may  be 
calculated  directly  from  Equations  C-3.1  through  C-3.4.  The  only  dif¬ 
ficulty  in  this  calculation,  albeit  a  minor  one,  is  that  the  value  of  ip^ 
must  be  calculated  from  Equation  C-3.4  by  an  iterative  process. 

These  calculations  are  made  in  Computer  Program  RC,  which  is  in¬ 
cluded  as  Appendix  D-l. 

Shear  Modulus 

The  basic  expression  for  the  dynamic  shear  modulus,  G,  is  as  follows: 

G  -  P  •  V*  (C-3.6) 

where  Vg  is  the  "S-wave"  (shear  wave)  velocity  of  propagation,  and  where 
p  is  the  mass  density  of  the  specimen  material  and  is  defined  in  Equation 
C-3.2.  The  value  of  Vg  may  be  calculated  as  follows: 

2*ir*  f  •  L 

Vs  -  X  ~ -  (C-3.7) 

where  fg  is  the  resonant  frequency  in  torsional  shear,  L  is  the  length 
(or  height)  of  the  cylindrical  specimen,  and  where  ip^  is  the  root  of  the 
frequency  equation: 

Vtan  K  “  J~ 

o 


(C-3.8) 


where  J  is  the  torsional  moment  of  inertia  of  the  specimen,  and  Jq  is  the 
torsional  moment  of  inertia  of  the  top  cap  system.  The  units  of  J  and  Jq 
are  mass-length-time2.  For  a  cylindrical  specimen,  J  may  be  calculated 
as  follows: 


J 


(C-3.9) 


where  D  is  the  diameter  of  the  specimen. 

As  discussed  earlier,  the  factor  <1^  is  necessary  to  correct  for  the 
dynamic  coupling  influence  of  the  top  cap  system  upon  the  measured  resonant 
frequency.  If  the  torsional  moment  of  inertia  of  the  top  cap  system  were 
negligible  when  compared  with  the  torsional  moment  of  inertia  of  the  spe¬ 
cimen,  -  ir/2,  and  Equation  C-3.7  would  reduce  to: 

V  -  4  •  f  •  L  (C-3.10) 

s  s 

In  the  more  general  case,  the  dynamic  shear  modulus  may  be  calculated 
directly  from  Equations  C-3.6  through  C-3.9.  The  value  of  in  this 
testing  series  is: 


Jq  ■  31.45  gra-cm-sec2 


(C-3.11) 


The  value  of  ^must  be  calculated  from  Equation  C-3.8  by  an  iterative 
process . 

These  calculations  are  made  in  Computer  Program  RC,  which  is  in¬ 
cluded  as  Appendix  D-l. 


J 
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Appendix  C-4 

Principal  Stresses  for  Resonant  Column  Specimen 

The  Mohr's  circle  diagram  in  Figure  C4-1  represents  the  state  of 
stress  within  the  resonant  column  specimen,  viewed  from  the  "@z  plane”, 
for  the  condition  of  simultaneous  maximum  vertical  compression  and  tor¬ 
sional  shear  stress.  The  principal  stresses  are  shown  on  this  figure  as 
a. i  0. ,  and  a  .  The  values  Ao  and  Atq  as  recorded  raw  data  were  double 

1  e  w  Z  V  Z 

amplitude  values  and  thus  must  be  divided  by  2  for  plotting  on  Mohr's 
circle . 

The  center  of  the  Mohr's  circle,  a  ,  may  be  calculated  as  follows 

Aa 

a  .  *  a.  +  — ~  •  sin(u)  t)  (3.18) 

Ct  0  4  z 


where  <*)  is  the  rotational  frequency  of  vertical  loading  and  is  constant, 
z 

Note  thats 


lc 


2c 


3c 


(C-4.1) 


for  this  testing  series,  where  CJ^c,  a2c»  and  c3c 
stresses  during  consolidation. 

From  the  figure, 


tan(2@) 


At, 


0Z 

.  —  •  sin(p*o)  *t) 
2  _  Z 


Aaz 

— —  •  sin(u  »t) 
.  4  z 


are  the  principal 


(C-4. 2) 


where  p 


Now: 


6 


r2  •  4,e, 

L  • 


•  sin(p-ci)  -t) 
Z 

sin(w  *t) 
z 


(3.17) 
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Now: 

Ac  At 

a.  *  a_  +  — —  •  sin(wt)  +  — — ■  •  sin(p*uj  *t)  •  tan(B)  (3.14) 
i  R  l  z  2  z 

Note  that  because  of  the  geometry  of  the  Mohr's  circle,  0^  may  also  be 
written : 

Ac 

C,  *  c  +  ——  •  sin(u)  *t) 

1  R  4  z 

+  rsc  atT 

%  t— r— •  •  sin(ai*t)  ]2  +  { — —5-  •  sin(p-w  *t)l  2  (C-4.3) 

4  Z  2  Z 

The  remaining  two  principal  stresses  are  as  follows: 

C2  -  0R  (3.15) 

and 

C3  -  2Cct  -  0L  (3.16) 

It  should  be  noted  that  C^  and  0 3  are  the  major  and  minor  principal 
stresses  at  the  condition  of  simultaneous  maximum  vertical  compression 
and  torsional  shear  stress  as  shown  in  Figure  C4-l(a).  As  time  goes  on, 
however,  the  amplitude  and  direction  of  these  principal  stresses  will 
change,  and  at  some  point  "c^"  will  be  numerically  smaller  than  "C3". 

At  that  point  the  principal  stress  directions  "reverse"  and  o3  becomes 
the  new  c^. 
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Appendix  C-5 

Principal  Stresses  for  Hollow  Cylinder  Specimen 


The  Mohr's  circle  diagram  in  Figure  C5-1  represents  the  state  of 
stress  within  the  hollow  cylinder  specimen,  viewed  from  the  "02  plane", 
for  the  condition  of  simultaneous  maximum  vertical  compression  and  tor¬ 
sional  shear  stress.  The  principal  stresses  are  shown  on  the  figure  as 
a^  and  o The  minor  principal  stress,  0 3,  is  unchanged  during  loading 
and  is: 

a3  "  a3c  "  (C-5*1) 


The  vertical  and  torsional  stresses  shown  in  Figure  C5-l(a)  may  be  written 

3 o  Ac 


sin(X-t) 


(C-5. 2) 


and 


3t 


9z 

3t 


At 


02 


sin  [(X-a)  •  t)  ] 


(C-5. 3) 


where  X  is  the  rotational  frequency  of  loading  and  is  constant. 


The  center  of  the  Mohr’s  circle,  a  ,  at  any  time  may  be  calculated 


as  follows: 


°ct  ~ 


Ac  Aa. 

[alc  +  — 2~  •  sin(X-t)!  +  [a2c  +  ~Y~  ’  sin(X*t)] 


(C-5. 4) 


and  because  a 


2c 


W(Cic  +  a3c> '  and 

Aa 


yAa 


(C5-5.5) 


ct 


.  talc  +  -5-  •  sin(X-t)}  +  [W(ale  +  a3c)  + 


lc  3c' 


•sin(X-t) ] 


which  reduces  to: 


ct 


(H^He  •  3in,X-t>  +  “T2  (3-3,) 


j 


FIGURE  C5-1  STRESS  TIME  HISTORIES  (a)  AND  MOHR'S  CIRCLE  IN 
STRESS  (b)  FOR  SIMULTANEOUS  DYNAMIC  LOADING  OF 
HOLLOW  CYLINDER  SPECIMENS  UNDER  CONDITION  OF 
MAXIMUM  INSTANTANEOUS  VERTICAL  AND  TORSIONAL  STRESS 
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Now,  from  Figure  C5-1,  the  following  expression  may  be  written: 


At 


9z 


tan(2B) 


•  sin[(A-a)  •  t] 


Act 

|alc  +  ~T  *  sin(X't)  *  °ct 


(C-5.6) 


or 


j  tan-1 


ATqz  •  sin[(A-a)  *  t] 


Ao 

2[*ic  +  -T 


sin(A't)  -  oct] 


(3.36) 


The  principal  stresses  may  now  be  written: 


Ao  Atq 

^  =  <Jlc  +  — 2~  ’  sin  (A- 1)  +  -  — -  •  sin[(A-ci)  •  t]  •  tan(B)  (3.33) 


=  2a 


ct 


a.  = 


-  2(a1  - 


act> 


(3.34) 


and  a3  «  03c  =  0R  (3.35) 

It  should  be  noted  that  with  time,  the  numerical  values  of  the  major 

A  0 

and  minor  principal  stresses  will  change.  If  — y  exceeds  0.85 -0.^  at  any 
time,  then  the  principal  stress  directions  will  "reverse"  for  a  period 
during  each  cycle  of  loading,  and  03  will  become  the  new  0^  for  that 


period. 
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Appendix  C-6 

Calculation  of  Compression  Modulus  for  Hollow  Cylinder  Tests 


From  Hooke’s  Law  the  strains  may  be  written  as  follows: 


-K 

+  V 

(C-6.1) 

+  V 

(C-6. 2) 

‘  I(CTZ 

+  v 

(C-6. 3) 

Because  eQ  =  0  in  this  testing  series,  Equation  C-6. 3  may  be  rewritten 
as  follows: 


CTe  *  W(a2  +  crR) 


Now  Equation  C-6.1  may  be  written: 


(C-6. 4) 


(C-6. 5) 


During  dynamic  loading,  the  change  in  vertical  strain  may  be 
written  as  follows: 


(C-6. 6) 


but  because  Ao  =  0, 

K 


,  Acr 

E  =  (1  - 


(C-6. 7) 

(4.13) 
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Appendix  C-7 

Calculation  of  Volume  of  Hollow  Cylinder  Specimen 

For  most  tests  in  this  testing  program,  the  specimen  height,  H, 
and  thickness,  T,  were  measured  immediately  following  the  initial 
isotropic  increment  of  consolidation.  For  some  earlier  tests,  only 
the  sample  height  was  accurately  measured.  Assuming  the  specimen  in 
Figure  C7-l(a)  is  typical,  then: 


T  - 


(C-7. 1) 


R  is  approximately  10.0  cm  (3.937  in)  for  all  samples, 
avg 

The  volume  of  a  specimen,  V,  is: 

V  -  H  *  7T  •  <R22  -  Rl2) 

T  T 

but  because  R_  *  R  +  —  ,  and  R,  *  R  -  —  , 

2  avg  2  1  avg  2 


Expanding  this  expression  gives: 


(C-7. 2) 


(C-7. 3) 


(C-7. 4) 


or  V  =■  H  •  tt  •  2  •  T  •  R  (C-7. 5) 

avg 

To  calculate  the  volume  in  cubic  inches,  since  T  and  H  are  measured 
in  inches,  Equation  C-7. 5  may  be  written: 


V  -  H  •  T  •  24.73695 


(C-7. 6) 


w  (Z  DIRECTION) 


FIGURE  C7-1  FREE  BODY  DIAGRAM  AND  STATE  OF  STRAINS 
FOR  HOLLOW  CYLINDER  SPECIMENS 


285 


If  T  was  not  recorded,  an  average  value  of  T  *  0.779  in.  may  be 
used  with  good  accuracy,  and  Equation  C-7.5  may  be  written: 

V  =  H  •  19.270084  (C-7.7) 

Equations  C-7.6  and  C-7.7  are  used  in  Comouter  Program  HC, 
included  as  Appendix  D-2,  to  calculate  the  volume  of  hollow  cylinder 
specimens . 


Calculation  of  Relative  Density 

The  soil  used  in  this  study  was  Monterey  No.  0  sand,  a  uniformly 

graded,  fine  grained  quartz  sand  processed  from  beach  sand.  A 

gradation  analysis  of  this  sand  is  shown  in  Figure  C8-1. 

In  order  to  calculate  relative  densities  for  the  various  specimens 

constructed  through  this  study,  laboratory  tests  were  performed  to 

determine  the  maximum  and  minimum  densities,  Y,  and  v.  .  ,  which 

d  max  d  min 

could  be  achieved.  These  two  limiting  densities  were  found  to  be  as 
follows : 


and 


Y.  =  1.70709  gm/cc 
d  max  * 

=  1.42532  gm/cc 

u  min 


(C-8.1) 

(C-8.2) 


The  relative  density  of  a  sand,  DR,  is  defined  as  follows: 

e  -  e 

D  =  - — -  x  100%  (C-8.3) 

R  e  -  e  . 
max  min 

where  e  is  the  void  ratio  and  is  defined  as : 

V 

e  =  ^  (C-8.4) 

s 

where  Vv  is  the  volume  of  voids,  and  Vg  the  volume  of  solids,  which  may 
also  be  written: 


W 


w 


In  this  equation,  W  is  the  weight  of 


negligible) ,  Gg  is  the  unit  weight  of 
unit  weight  of  water,  y  ■ 


(C-8.5) 

soil  (weight  of  voids  is  assumed 

the  solid  materials,  and  y  is  the 

w 


1  gm/cc.  The  volume,  V,  may  be  written: 
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Appendix  D 


Computer  Programs 


The  following  computer  programs  were  developed  for  use  in  this 
study: 

D-l  Program  RC — A  program  to  Reduce  Data  from  the  Triaxial  Resonant 
Column  Testing  Series. 

D-2  Program  HC--A  Program  to  Reduce  Data  from  the  Thin-Walled  Hollow 


Cylinder  Testing  Series. 


i o  n o on nonn 


APPENDIX  D1 

Program  RC— A  Program  to  Reduce  Data  from  the 
Triaxial  Resonant  Column  Testing  Series 


PROGRAM  BC  1  INPUT  .OUTPUT  .PU"CH>  1" 

r  !• 

c  A  COMPUTES  PRjGSam  Ts.  PEOuCE  oat*  ercn  the  TP  I  ax  I  al  RESONANT  CCLjHN  J- 

r  TeatinG  'ERIE').  ** 

f  s: 

f  JHfMO  °v  PAT  GOt**!“  —  h-th^cr  r'»  l«TO  A* 


c  this  program  is  »as'o  on  ilassical  elastic  t“c'c(  at.  hoc<e.s  l»*. 
c  the  input  parameters  include  the  sample  physical  dimensions,  the 
C  SYSTEM  DYNAMIC  PROPERTIES,  calibration  factors  AND  The  RA„  test  sc- 

c  DATA,  INCLUDING  THE  DRIVISS  voltage.  AfCELEROMeT-S  RESPONSE  VDlTaqe,  SSC 

c  AND  THE  RESONANT  EBEDUErCY  In  VERTICAL  and  THRSIHMAL  L,'*0I"C.  The  ’  »C 

r  program  COMPUTES  The  oomoressidn  ano  shear  yhh.j|_|,  strain  AMPLITUDES  13: 

c  AND  DAMPING  RATIOS.  THE  MAGNITUDE  AN'  DIRECTION  Dc  ThE  PRINCIPAL  !*- 

C  STRESSES  AND  Thc  INCRevcn  AL  PRINCIPAL  STRAINS.  AND  T“f  D-TAHEDRaL  1  S  1 

C  STRAIN  and  the  OCTRahEDRal  shearing  strain.  «“ERE  ThE’E  I A  A  1‘- 

C  REOjNDANCT  or  DATA,  such  AS  with  THE  COYS  I  NED  LOADING  TESTS.  The  IV 

r  PRDGRay  ALSO  Compute,  i  POISSON'S  RATIO  A.'O  the  Ratio  oc  1  S' 

C  HOROZONTAL  Shear  strain  to  VERTICAL  normal  strain.  is- 


r 

r 

7*r  rocLOWfNC-  f$  *  ne  Tnr  VA»!lal.?  N*-vrS  JS”"  1  •)  T-f? 

r 

*V8 

•9  *»“ 

r 

*• 

*»4" 

r 

TORS'DNAL  acceleration  Resodnae  re 4 o I ng 

* 

c 

4Cf?V 

vertical  acceleration  RESPONSE  reading 

26- 

c 

ACFT 

torsional  ACCELERaOICN  calibration  eactcr 

27: 

c 

*crv 

VERTICAL  ACCELERATION  C»Ll“E*TIDN  'ACTD® 

26* 

c 

9FT4 

ANGLe  OE  RE-ORIENTATION  OE  The  ooinCI  =  »L  STRESS  «'! 

29* 

/■ 

^0»l  T 

TORS'DNAL  hRtv'Ng  vDLTagb  saam«!'. 

c 

con  v 

VFRTtCtL  DRIVING  VOLTAGE  reading 

11* 

c 

04VOT 

TOPSIRaial  OAmo;«-g 

ti¬ 

c 

DAMPV 

vertical  OAmoInG 

ll' 

r 

or  ft 

TORSIONAL  damping  GALIRBATTON  eacTdr 

14* 

c 

ocrv 

vertical  OAMoing  -ALIRRATIO'  ei-T-5 

'1*. 

c 

Deo<V 

VERTICAL  ‘TRAIN 

u* 

c 

0G**'T 

torsional  «t»ai" 

1  >•' 

z 

OU* 

Ot*“ETER  «E  5AMOI.E 

u* 

05I3V 

VERTICAL  'TRESS 

c 

OTijT 

TORSIONAL  stress 

4-- 

c 

?4T4 

ANGLE  OE  RE-ORIENTATION  OE  The  in-REvenTal  PRINCI=»L  STRAINS 

4l- 

c 

•PSCT 

CENTER  OF  EPS I  -  EPS2  MOHR'S  CIRCLE 

42* 

c 

E°SOC 

OCTAHEDRAL  STRAIN 

43  r 

?©$1 

major  INCREMENTAL  RRInCIPal  STRAIN 

44 ; 

* 

?PS2 

Intermediate  incremental  0RINCt°*L  STRAIN 

46: 

c 

W51 

MINOR  INCREMENTAL  PRINCIRAL  strain 

4f>"' 

c 

FO^T 

TORSIONAL  RESONANT  EReO'JenCY 

4*’~ 

pp 

PPOV 

vertical  resonant  creojeaj'y 

49' 

c 

G4«9C 

octahedral  shearing  strain 

49- 

c 

H5TGMT 

height 

6 : : 

c 

MOOT 

SHEAR  MODULUS 

51* 

r 

MOOV 

YOUNG'S  MODULUS 

52: 

ROIS'ON'S  RATIO 

*  ^  - 

* 

ViJ^ 

C»LC'JL»TEO  POI  S  SD*t .  <  0  A  T  I  — 

*  4  * 

c 

*P 

number  0*  DATA  POINTS  for  a  sample  (vaxImiu  50  J 

56: 

r 

MS 

number  H»  samples 

r 

PM| 

R0Lc»  mass  OE  INERTIA  OE  sample 

4  7- 
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MASS  DENSITY  0 1  SA«Pl* 

CFMT*B  0*  SIS!  -  i!0>  -DNS'S  T!»'LC 
majOB  PBINCI°»l.  aiBEsa 
j«4T«»NEDI  ‘  IE  ®91NCIP*I.  aiBPaa 
ul^g  ooIvOIOAt  Sr**?S 
CONFINING  °»E$SU°? 

eyeT?v  ui<s  OATI" 

SYSTEM  POLAR  NASS  3F  ! NTS  t  I  A 
Systfv  »Ot.»*  -ASS  0»  tsruriA  BaTIO 
BAT  ID  Of  AVG*  DGAmi  io  DEBa* 
mo  CAB  WEIGHT 

unit  weight  or  samou 

VO(.U'<r  ho  SA-Blr 
VEIOMT  or  SAMB|.r 


DIMENSION  COILV(SOl.ACCEVISO) tFRCiVl S3 > .C3ILT ( 1 3 I .ACCET ( SI i ,BB3T ( SC 

! ) 

BEAU  V'JtV:jC«vOOV*MODT 
COOVAT I  I « I 

I  *  I 

mo-AT  ( on  j.A  i 
fh9«aT(  1  h*  •/////! 

rCSM»T(l3X.'Tr«T  pceuLTS  BOB  5A-PL*  NO.1.!*) 

Fhp«ATi/////I*X.151S*G  »’  01*''  HEIGHT  NP</ 

;iax.<xSC  GBAvS  C«  CM*) 

eOB-AT(/.lGX.triG.2.I4l 

rOBM»T(///HX»  *TCw  »crv  DG*Y  »CrT  5C  T 

txp/IAXi  •<■,»»»«•  I 

e'y»'*ATC/»!Ox.s*!:.?i  .  . 

FORMAT  t // / l AX • • VOL  yxT  S»HI  SMB  SB  I  / 

ITAXt'CC  G/CC'I 

C0BVATu:x.4Pl*.?l  ...  -  .  *  r 

pOBvaT(///I*x..c0H.V  * C'*M  *BOV  D  '  ItTvu  „*  •  uesT,, 

IBOT./!AX.*«V-»o  »V-»°  hybt:  mv-bb  mv-bb  -pb  ’• 

’pOr/ATIi.x .  WAO..t««  -  «!«!  IP  N-T  Y.T.Y.'Y  ST«1  **••! 

ph».,»TI/////IOX..CAL-JL»tx-  aT»’r  aM'  ru»«Tl-  C3N«  ANTS,  *A 

it*  seT  NO.>.n///nx»'«oov  ded'v  o*wt,'/  OT 

,T  TavflT  -UC  STB«T*/1!«.  ........ 

;,0SI  PERCENT  PERCENT  RSI  PERCENT  B.bc.NT  /■  I 

poo-aTi pbbob  in  VERTICAL  -GOJLJS  ' ALE JLt I  ION  •»•«) 
POBMAT(I«x,*p»»  EBB3B  IN  T3R5I3NAL  vOOJLJS  CALCJIAIION  ••••) 

fORvaTI SX.C1 3»?»510.4*?ri3,2»*l''»*«SrlI«A* 

C^5v|  T  f •**  *  T  **•  ’«ff  *  c*  '  •  9  f 

I ?*X *r 1  - .? • 51  - •4»rl 

f 

POBv  AT  I//4X»'MJ  SIG1  BIG?  I!  COC' 

,  cpc;  Pp3)  p  4  T  A  EBSCC  C»VCC«/14*. 

2.BSI  PSI  »<!  3E0BCC  5  bpb:pnT  Pebe-nt  BEBCPN 

,t  oe09£CS  ar»C**r</i 

CrtOMjT  •  9  *  ,r *  •  '*  , 


;?itT  c 

D  VV-»©  wtOT’t 


■jc;:  IMO-JT  }AT* 


?5A}<»«*nS 
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APPENDIX  D2 

Program  HC--A  Program  to  Reduce  Data  From  the 
Thin-Walled  Hollow  Cylinder  Testing  Series 
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